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Abstract 


A  lot  of  work  is  devoted  to  the  study  of  Ion  Assisted  jieposition  (LAD)  with  different 
materials:  metals,  Al,  Ag,  and  dielectrics,  nitrides,  fluorides  and  oxides. 

Concerning  oxides,  a  comparison  between  IAD  and  I6n  Hating  techniques  is  in  progress. 

For  the  characterization,  we  describe  the  two  techniques  used  for  determination  of  optical 
losses  which  are  photothermal  deflection  and  measurements  of  angular  scattering  distribution. 

Concerning  photothermal  deflection,  the  results  of  characterization  obtained  on  the  OSA 
samples  of  1102  are  given;  instabilities  versus  time  are  observed  when  the  samples  are  illuminated. 

The  theory  and  experiments  performed  on  scattering  are  applied  to  coatings  of  special 
designs  produced  by  Balzers,  OCLI,  OJAI  and  Spectra  Physics  to  determine  the  correlation 
between  interfaces  of  the  multilayers. 

These  results  .which  were  presented  at  the  OSA  1989  meeting^in  Florida  are  summarized  . 
here.  This  work  leads  indirectly  to  the  measurement  of  the  grain  size  of  the  materials  in  thin  film 
form.  A  completely  different  approach  is  given  by  growth  modeling,  but  it  remains  to  link  the  two 
techniques  for  a  better  understanding  of  the  process  of  growth. 
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INTRODUCTION 

For  several  years,  and  due  to  US  Air  Force  grants,  the  two  laboratories  of  Prof 
H.A.  Macleod  at  the  Optical  Sciences  Center  -  Tucson-  and  E.  Pelletier  at  the  Ecole 
Nationale  Sup&ieure  de  Physique  -  Marseilles-  have  coordinated  their  work  of 
research,  playing  in  this  way  complementary  roles. 

The  distance  between  Tucson  and  Marseilles  brings  obviously  some  difficulties 
because  the  number  of  meetings  between  the  members  of  the  two  groups  is 
necessarily  limited.  This  year  we  should  have  waited  for  the  end  of  the  annual  OSA 
meeting  in  Orlando  (oct  89)  to  take  stock  of  the  situation  for  the  works  as  part  of  the 
cooperative  program  and  to  draw  the  future  prospects. 

The  present  report  shows  first  the  evolution  of  the  work  in  each  of  the  two 

labs: 

-  the  Tucson  group  is  interested  in  the  development  of  Ion  Assisted  Deposition 
(IAD)  techniques,  and  more  particularly  for  metals  (Al,  Ag)  and  dielectrics  (nitrides, 
lanthanides  trifluorides)  interesting  for  their  potential  application  to  infrared 
optical  components.  The  results  obtained  on  IAD  induced  the  group  to  deal  with  film 
growth  modeling  to  understand  in  detail  the  growth  mechanisms  effectively 
involved. 

-  with  a  complementary  but  totally  different  approach,  the  group  of 
Marseilles  is  also  interested  in  this  same  problem. 

In  comparison  with  the  Tucson  group,  we  bring  few  results  on  the  IAD 
technique  because  our  activity  is  limited  on  the  study  of  very  classical  materials 
such  as  Ti02.  Si02  and  Ta20s.  However,  in  Marseilles,  we  began  a  comparative  study  of 
different  deposition  techniques  since  we  have  at  our  disposal  a  Balzers  Ion  Plating 
plant.  To  this  we  can  add  important  attempts  in  the  development  of  optical 
characterization  techniques.  Here  we  will  deal  particularly  with  two  subjects  that 
seem  to  be  more  prominent  for  this  year  and  that  are  respectively 

-  measurement  of  absorption  losses  by  photothermal  technique 

-  and  mainly,  measurement  of  the  scattering  curve  to  look  into  the  question  of 
correlation  between  successive  interfaces  of  a  layer  stack.  On  this  last  question,  we 
will  made  a  rapid  evaluation  of  the  progresses  made  and  we  will  especially  present 
the  results  obtained  with  the  active  cooperation  of  the  OSA  Technical  Group  of  Thin 
Films  led  by  J.  Apfel  (O.C.L.I.).  We  meet  the  concern  of  Tucson  on  the  study  of  layer 
growth  mechanism,  where  we  will  have  to  combine  our  efforts  to  progress. 
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ION  ASSISTED  DEPOSITION  AND  MODELING  AT  THE  UNIVERSITY  OF 

Arizona 

Relevant  work  at  the  University  of  Arizona  during  the  period  of  this  contract 
has  been  supported  by  the  USAF  under  the  URI  program.  This  has  placed  emphasis 
on  the  investigation  of  the  use  of  ion-assisted  deposition  for  modification  of  film 
structure  and  determining  its  utility  for  improving  the  performance  in  optical  thin 
films.  The  materials  studied  by  ion-assisted  deposition  has  ranged  over  metals, 
especially  silver  and  aluminum,  and  dielectrics,  especially  aluminum  oxide,  nitride 
and  oxynitride  and  the  fluorides  of  the  lanthanide  series.  In  modeling,  some 
progress  has  been  made  with  the  Monte  Carlo  model  of  film  growth.  Much  of  the 
Arizona  work  reported  here  is  also  the  subject  of  reports  to  the  USAF  from  the  URI- 
supported  Center  for  Thin  Film  studies. 

Ion-assisted  deposition  of  metals 

It  was  of  interest  to  us  to  determine  if  the  great  benefits  of  io-assisted 
deposition  in  the  fabrication  of  dielectric  coatings  would  also  be  realised  in  the 
deposition  of  metals  of  high  optical  performance,  such  as  silver  and  aluminum.  The 
results  were  disappointing  in  that  the  optical  properties  were  found  to  be  degraded 
in  the  sense  that  the  real  part  of  refractive  index  rose  with  bombardment  while  the 
extinction  coefficient  fell.  Typical  results  for  silver  are  shown  in  figures  1  and  2. 
Note  that  the  changes  in  the  optical  constants  accompany  increases  in  the  implanted 
argon  of  the  films.  Broadly  similar  results  were  obtained  for  aluminum  although  the 
changes  for  very  small  bombardment  were  not  as  pronounced,  figures  3  and  4, 
making  it  possible  that  very  light  bombardment  might  be  acceptable  for  aluminum. 
Although  the  changes  are  not  large,  nevertheless  they  are  in  a  negative  sense. 
Similar  changes  were  found  in  other  aspects  of  film  performance,  the  DC  resistivity 
for  example  (which  rose  with  bombardment).  Implanted  argon  also  correlated  well 
with  the  other  changes.  The  conclusion  was  that  straightforward  ion-assisted 
deposition  was  not  as  favorable  a  process  for  metal  film  deposition  as  it  is  for 
dielectrics. 

The  experiments  up  to  this  stage  were  all  involving  ion  energies  in  the  range 
200  to  500  eV.  Before  leaving  the  question  of  metal  films  we  thought  that 
experiments  using  lower  bombardment  energies  should  be  carried  out.  Furthermore, 
although  the  optical  property  changes  appear  disappointingly  negative,  there  are 
aspects  of  ion-assisted  deposition  of  metals  that  are  usefiri,  particularly  the 
improvement  in  adhesion  that  affects  both  metals  and  dielectrics.  The  bombardment 
that  takes  place  before  and  during  the  early  stages  of  film  growth  is  particularly 
effective  in  this  respect.  A  series  of  experiments  were  therefore  undertaken  to 
investigate  the  use  of  lower  energies  in  bombardment,  20  to  50  eV,  and  also  the  effect 
of  limiting  the  bombardment  to  the  growth  of  the  first  part  of  the  films  only. 


Momentum 


Figure  1 

Variation  of  n  and  k  at  632.8  nra  of  Ag  films  measured  by  the  SPR  as  a  function 
of  P.  where  P  =  y  V~E  and  y  is  the  ratio  of  arrival  rates  of  ions  and  evaporant  atoms. 
(Hwangbo) 
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Since  we  found,  in  the  higher  energy  series  of  experiments,  that  the  DC  resistivity  of 
the  films  was  a  good  indicator  of  film  quality,  this  was  the  measured  parameter  in 
this  new  series. 

Results  are  shown  in  figure  5.  All  films  had  the  same  total  thickness  of  55 
nm.  The  vertical  scale  is  the  measured  resistivity  of  the  films  after  deposition.  The 
horizontal  scale  is  the  thickness  of  film  that  was  subjected  to  bombardment  after 
which  the  bombardment  ceased  and  the  remainder  of  the  film  was  deposited  without 
interruption  by  conventional  evaporation  without  bombardment.  The  results  were 
not  very  consistent,  although  the  error  bars  calculated  for  the  experiment  (shown 
against  only  one  set  of  results  ior  clarity)  were  very  much  less  than  the  variations 
exhibited.  Zero  bombardment  gave  consistently  the  value  of  resistivity  that  had  been 
established  before.  Bombardment  of  the  entire  film  thickness  led,  as  might  be 
expected  to  resistivities  that  increase  both  with  the  energy  and  current  of 
bombarding  ion  beam.  For  bombardment  of  the  initial  stages  of  film  growth  only, 
however,  there  is  a  tendency  for  the  resistivity  to  show  a  slight  decrease  compared 
with  the  unbombarded  film.  The  decrease  is  unfortunately  not  well  marked,  and 
beyond  a  5  nm  thickness  of  bombarded  material  the  scatter  in  the  results  increases 
considerably,  bombardment  at  50  eV  and  low  ion  current  actually  showing  a 
substantial  decrease.  This  particular  result  might  be  worth  further  investigation.  At 
this  stage  we  can  conclude  from  the  results  of  ion-assisted  deposition  of  the  early 
stages  of  a  metal  film  appears  acceptable  and  that  this  should  make  it  possible  to 
obtain  some  of  the  benefits  of  the  ion-assisted  process  without  suffering  the 
degradation  that  can  also  accompany  it. 

Ion-assisted  deposition  of  nitrides 

The  usefulness  of  aluminum  oxide-nitride  materials  was  illustrated  by  the 
construction  of  a  simple  coating  of  the  rugate  type  where  the  index  contrast  was 
approximately  1.83  to  1.66.  The  control  of  the  process  was  only  semiautomatic  and 
therefore  instead  of  a  continuous  variation  of  refractive  index  a  10-step 
approximation,  figure  6,  was  used.  The  resulting  filter  characteristic  curve  is 
shown  in  figure  7.  In  the  construction  of  this  filter  the  only  variable  was  the 
background  pressure  of  oxygen  in  the  coating  plant.  The  evaporant  material  was 
aluminum,  evaporated  at  a  constant  rate.  The  bombarding  ions  were  nitrogen.  Since 
uxygen  is  much  more  reactive  than  nitrogen,  the  composition  of  the  resultant  film 
could  be  varied  right  up  to  virtually  pure  aluminum  oxide  simply  by  varying  the 
background  oxygen  pressure. 

More  recently  we  have  been  investigating  the  nature  of  the  aluminum- 
oxygen-  nitrogen  material  with  the  primary  objective  of  determining  whether  the 
material  produced  by  our  ion-assited  deposition  is  simply  a  mixture  of  aluminum 
oxide  and  aluminum  nitride  or  a  true  aluminum  oxinitride  compound.  If  it  were  a 


Resistivity  0*0- 
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Figure  5: 

Resistivity  of  low  energy  IAD  silver  films  as  a  function  of  the  thickness  of  the 
bombarded  portion.  All  films  were  55  nm  in  total  thickness.  (Hwangbo) 


Physical  Thickness  (nm) 


Figure  6: 

A  step-index  profile  of  A10*Ny  films  from  Fig.  1  for  a  rugate  filter.  One  cycle 
consists  of  10  homogeneous  step-index  layers  and  has  a  168.6  nm  physical  index.  A 
sinusoidal  curve  is  fitted  to  the  step-index  profile.  (Hwangbo) 
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Wavelength  (nm) 


Figure  7: 

Reflectance  (R)  and  transmittance  (T)  of  the  real  and  simulated  rugate  filters. 


mixture  then  we  might  expect  the  x-ray  photoelectron  spectra  to  show  evidence  of 
this.  The  aluminum  2p  peaks  of  the  compound  film  show  no  evidence  of  two  separate 
compounds  involving  aluminum.  The  same  appears  true  of  the  oxygen  Is  peaks  but 
the  nitrogen  spectrum,  figure  8,  is  more  complicated.  At  low  nitrogen  content, 
there  is  an  extra  nitrogen  peak  at  a  somewhat  higher  energy  than  the  main  peak. 
This  peak  has  been  seen  by  other  workers  when  aluminum  nitride  was  baked  at 
elevated  temperatures  and  attributed  to  aluminum  oxynitride.  Since  we  see  the  peak 
only  when  the  nitrogen  content  is  low  and  the  oxygen  content  high,  we  believe  it 
may  be  rather  an  interaction  between  nitrogen  and  oxygen  directly.  The  aluminum 
and  oxygen  peaks  at  this  composition  tend  to  show  the  same  energy  as  AI2O3.  At 
higher  nitrogen  concentrations  the  small  second  peak  disappears  and  we  appear  to 
have  true  aluminum  oxynitride.  A  fuller  discussion  of  these  interesting  results, 

including  the  complete  set  of  spectra,  is  to  be  found  in  Hwangbo’s  dissertation. 

Lanthanide  trifluorides 

The  lanthanide  trifluorides  have  been  the  subject  of  investigation  by  the  Thin 
Film  Group  at  the  University  of  Arizona  for  some  time.  They  are  materials  with  great 
promise  for  optical  coating  applications  because  of  their  wide  range  of  transparency 
from  the  extreme  UV  through  to  the  far  IR.  They  are  of  low  refractive  index  in  the 
infrared  and  this  increase  still  further  their  potential  usefulness.  Unfortunately 
they  suffer  from  exceedingly  high  intrinsic  tensile  stress  when  conventionally 

deposited,  stress  that  is  sufficiently  high  to  lead  to  tensile  failure  when  layers  are 
thick  enough  to  be  useful  in  the  infrared,  or  when  they  are  used  in  multilayers.  Ion- 
assisted  deposition  is  the  only  process  that  promises  complete  stress  control  in 
coatings  and  therefore  it  was  of  interest  to  investigate  the  use  of  ion-assisted 
deposition  in  studying  the  lanthanide  fluorides.  Briefly,  some  of  the  series  are 
polymorphic  and  in  their  deposition  the  predominant  crystalline  form  which  we 
obtained  at  low  substrate  temperatures  was  the  high-temperature  structure,  while  at 
high  substrate  temperatures  we  found  a  predominance  of  the  low  temperature  form. 
This  result  is  a  consequence  of  the  nonequilibrium  nature  of  the  process  and 

especially  the  rapid  quenching.  Particular  attention  was  paid  in  the  ion-assisted 
experiments  to  samarium  fluoride.  Figure  9  shows  the  infrared  transmittance  of 
both  conventionally  deposited  and  ion  assisted  samarium  fluoride.  Conventional 

deposition  of  samarium  trifluoride  at  room  temperature  gives  the  high  temperature 
hexagonal  phase  with  the  (111)  orientation  stronger  than  the  (110).  Increasing  the 
substrate  temperature  to  200°C  and  above  led  to  the  appearance  of  the  low 
temperature  orthorhombic  phase  which  tended  to  show  a  (111)  orientation.  Ion- 
assisted  deposition  at  low  temperature  should,  if  the  mechanism  is  the  same  as  that  of 
substrate  heating,  tend  to  produce  the  orthorombic  phase.  No  such  tendency  was 
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XPS  measurements  of  the  main  Nls  peaks  of  A10*NV 
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Figure  9: 

The  infrared  transmittance  of  films  of  samarium  fluoride  (Lingg) 
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observed.  Low  temperature  ion-assisted  deposition  gave  the  hexagonal  phase  but  the 
strong  orientation  was  (110)  with  some  evidence  of  (111).  Increasing  the  degree  of 
bombardment  still  did  not  lead  to  the  orthorhombic  phase  but  eventually  resulted  in 
the  appearance  of  x-ray  diffraction  peaks  characteristic  of  a  number  of  samarium- 
fluorine  compounds  with  F  to  Sm  ratios  less  than  three  together  with  a  further 
reduction  in  the  hexagonal  (111)  intensity.  Thus  ion-assisted  deposition  leads  to  a 
change  in  the  preferred  orientation  of  the  crystalline  structure  rather  than  a 
complete  change  in  structure  and  the  mechanism  does  not  appear  to  be  simply 
thermal.  In  fact  in  our  previous  reports  we  have  described  the  evidence  that  leads  us 
to  believe  that  a  principal  mechanism  in  the  compaction  of  the  growing  film  is 
transfer  of  momentum.  We  are  not,  however,  in  a  position  to  state,  unequivocally  that 
it  is  this  same  momentum  transfer  that  is  responsible  for  the  changes  in  crystalline 
orientation,  but  just  that  the  mechanism  is  not  a  straightforward  thermal  one.  The 
investigation  of  the  lanthanide  series  is  described  in  Lingg's  dissertation. 

Film  growth  modeling 

The  statistical  properties  of  surfaces  of  films  grown  using  a  Monte  Carlo  model 
where  the  adatoms  have  a  certain  mobility  on  the  surface  after  arrival.  This  is  a  good 
approximation  to  the  more  realistic  case  of  thermally  activated  hopping.  The  atoms 
that  take  part  in  thermally  activated  hopping  are  confined  largely  to  the  surface 
and  arc  inhibited  from  further  hopping  when  they  are  buried  under  freshly 
arrived  material.  Thus  the  mobility  in  the  Monte  Carlo  model  is  to  be  considered  to  be 

the  permitted  amount  of  movement  of  the  adatom  before  it  becomes  pinned  under 

further  material  and  is  an  increasing  function  of  deposition  temperature  and 
decreasing  function  of  deposition  rate.  Larger  scale  deposition  models  become 
possible  with  the  scheme  than  could  be  achieved  by  the  thermally  activated  hopping 
model.  It  then  becomes  interesting  to  examine  the  surface  statistics  and  eventually  to 
compare  them  with  measurements  on  real  films.  Figure  10  onwards  illustrate  some 
of  the  results  and  a  summary  is  given  in  Table  1.  The  conclusions  that  we  can  draw 
from  this  study  so  far  is  that  over  several  orders  of  magnitude  the  simple  ballistic 
deposition  model  predicts  that  thin  film  surface  have  a  fractal  nature,  that  the  low- 
frequency  cut-off  of  the  fractal  behaviour  appears  to  decrease  as  film  thickness 
increases  and  that  increasing  adatom  mobility  appears  to  decrease  the  fractal 
dimension.  Just  as  an  information,  figures  11  to  16  give  an  idea  of  the  results  given 
by  our  model  of  film  growth.  We  will  find  the  problems  of  surface  statistics  at  the 
scale  of  optical  defects  in  the  section  devoted  to  scattering  studies  in  Marseilles.  More 

information  on  this  topic  will  eventually  be  published  in  Sargent's  PhD  dissertation, 

in  preparation  at  present. 
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GR0H89B  Cprogran  date  2  June  19891 

RUN  ID:  MIMSURF-000  DATE:  9  July  1989 

Surface  Mobility:  Tn/Ts  =  0  (Henderson-type  Simulation). 

1000x500 


Figure  10: 

Simulation  of  film  grown  under  conditions  of  zero  mobility.  (Sargent) 
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Properties  as  a  function  of  adatom  mobility 


Mobility 


rm/rd 

chi 

n 

K 

D 

.dens 

avg 

rms 

0 

0.59 

1.88 

.68 

1.56 

.63 

463 

12.19 

1 

1.10 

1.93 

.48 

1.54 

.72 

463 

11.50 

2 

1.19 

1.90 

.67 

1.55 

.74 

465 

11.45 

4 

1.44 

1.99 

.40 

1.51 

.77 

465 

12.73 

8 

1.92 

2.03 

.18 

1.49 

.82 

479 

7.29 

16 

2.41 

2.08 

.14 

1.46 

.84 

481 

6.07 

32 

2.77 

2.10 

.15 

1.45 

.87 

482 

6.63 

64 

2.98 

2.05 

.20 

1.48 

.90 

480 

8.67 

chi: 

average  surface  diffusion 
adatom  diameters. 

length 

in 

n: 

value  of  power  in  inverse 

power 

law 

K: 

spectrum. 

value  of  constant  in  inverse 

law 

power 

D: 

spectrum. 

Hausdorff  fractal  dimension. 

dens: 

density  of  structure  relative 
close-packed. 

to  2D 

hep 

avg: 

average  height  of  surface  (in  cells). 

rms: 

rms  deviation  of  surface  (in  cells). 

Table  l. 

Results  of  a  number  of  film  growth  simulations. 
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Figure:  Surface  power  spectra  of  run  minisurf-000  for  nominal  thicknesses  of  100,  200, 
300,  400,  and  500  cells. 


Figure:  Surface  power  spectrum  of  run  minisurf-000  plotted  on  a  log-log  scale  for  a 
nominal  thickness  of  500  ceils. 


Figure  II: 

Statistics  of  the  film  surface  in  figure  10.  (Sargent) 


Surface  Power  Spectrum 


Figure:  Surface  power  spectra  of  run  minisurf -008  for  nominal  thicknesses  of  100.  200. 
300,  400,  and  500  cells. 


Figure:  Surface  power  spectrum  of  run  minisurf-008  plotted  on  a  log-log  scale  for  a 
nominal  thickness  of  500  cells. 


Figure  13: 

Statistics  of  film  surface  in  figure  12.  (Sargent) 


2 
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Surfaca  Power  Spectrum 


Figure:  Surface  power  spectra  of  run  minisurf -064  for  nominal  thicknesses  of  100,  200, 
300,  400.  and  500  cells. 


Figure:  Surface  power  spectrum  of  run  minisurf -064  plotted  on  a  log-log  scale  for  a 
nominal  thickness  of  500  cells. 


Figure  15: 

Statistics  of  the  surface  in  figure  14.  (Sargent) 


Surface  Power  Spectrum 


Figure:  Surface  power  spectra  of  run  surf-002  for  nominal  thicknesses  of  100,  200,  300, 
400,  and  500  cells. 


Figure:  Surface  power  spectra  of  run  surf-002  plotted  on  a  log-log  scale  for  nominal 
thicknesses  of  100,  200,  300,  400,  and  500  cells. 


Figure  16: 

Statistics  of  a  very  large  simulation  showing  the  variation  of  the  statistical 
parameters  with  film  thickness.  (Sargent) 
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Comparative  study  of  ion-assisted  Deposition  IAD)  and  reactive 
ion  plating  (RIP)  techniques  for  production  of  T1O2,  ta20s,  S1O2 

At  the  present  time  IAD  and  RIP  techniques  are  currently  used  in  many 
optical  laboratories;  we  give  here  only  a  short  review  concerning  these  two 
techniques  before  the  analysis  of  the  optical  characterization  concerning  the  layers 
produced  in  Marseilles  Laboratory. 

Ion  Assisted  Deposition  (IAD) 

As  in  the  case  of  classical  reactive  evaporation,  the  deposition  is  made  by 
evaporation  of  materials  with  an  electron  gun.  The  relative  position  of  the  electron 
gun  with  respect  to  the  substrates  has  been  adjusted  to  obtain  the  best  deposition 
uniformity.  The  control  and  the  regulation  of  the  deposition  rate  is  obtained  with  a 
quartz  microbalance.  For  the  deposition  of  metallic  oxides  as  TiC>2  or  Si02  it  is 
necessary  to  operate  in  oxygen  reactive  pressure  servo  controlled  on  the  O2  partial 
pressure  (  with  a  quadripolar  mass  spectrometer)  or  on  the  total  pressure.  In  the 
case  of  a  ion  beam  assistance,  it  is  not  necessary  to  heat  the  substrates;  the  layer  in 
formation  is  bombarded  with  a  beam  of  Argon  and  Oxygen  ions  with  energies  of 

some  hundreds  of  eV. 

The  hot  cathode  ion  gun  we  use  includes  a  8  cm  in  diameter  output  grid  and  is 
40  cm  apart  from  the  monitoring  substrate  center  (Fig  17). 

We  can  consider  that  the  ions  exchange  their  energy  with  the  condensing 

material  by  transfer  of  momentum.  The  effect  obtained  depends  mainly  on  the  ion 
mass,  the  ion  energy  and  on  the  ion  current  density  on  the  coated  surface. 

Electrons  are  coupled  to  the  ion  beam  with  a  filament  in  such  a  way  to  neutralize 
the  total  space  charge.  We  use  defocused  grids  to  obtain  a  slightly  open  ion  beam. 

Reactive  Ion  Plating  (RIP) 

The  Reactive  Ion  Plating  processus  is  a  combination  of  evaporation  and 
pulverization  (Fig  18).  The  coating  apparatus  is  a  BAK  800  IP  from  Balzers. 

The  setting  of  the  different  elements  in  the  vacuum  chamber  is  practically 
the  same  as  for  evaporation  techniques,  but  the  substrate  holder  is  electrically 
isolated.  The  electron  guns  are  isolated  from  ground  and  a  potential  of  about  40  V  is 
established  between  a  heated  Filament,  emitting  electrons,  and  the  evaporation 

crucible  in  process.  The  argon  gas  is  introduced  close  to  the  Filament;  the  oxygen 
gas  is  introduced  close  to  the  crucible  and  the  discharge  is  established  between 
Filament  and  crucible. 


Figure  17: 

Ion  Assisted  Deposition 


Insulated 

rotory 


Figure  18: 

Reactive  Ion  Plating 


electron  guns 


Plasma 
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The  evaporation  is  performed  in  the  presence  of  the  gas  discharge.  In 
collisions  and  electron  impact  reactions,  coating  material  ions  are  formed  and 
accelerated  in  the  electric  field  so  that  condensation  and  film  formation  take  place 
under  the  influence  of  ion  bombardment.  This  complex  process  is  not  well  defined 
but  has  numerous  advantages:  in  addition  with  acceleration  of  coating  material  in 
the  plasma  collisions,  charge  effects  make  the  reaction  between  materials  easier. 
This  process  is  especially  efficient  for  the  formation  of  oxides  and  nitrides.  The 
relatively  high  gas  pressure  in  the  evaporation  chamber  (10‘3  Torr)  and  the  large 
size  of  the  crucibles  give  uniform  coatings  on  a  large  area. 


Characterization  of  optical  performances 
The  main  interest  of  the  IAD  and  RIP  techniques  in  comparison  with  classical 
reactive  evaporation  is  that  the  layers  produced  are  not  sensitive  to  hygrometry 
when  the  parameters  of  deposition  have  been  well  adjusted.  One  can  demonstrate 

that  the  index  shift  during  air  admission  is  completely  negligible. 

For  these  layers  we  will  present  our  results  concerning  characterization  and 

obtained  in  the  study  of  three  widely  used  materials:  Ti02,  Ta20s  and  Si02- 

For  each  deposition  technique,  we  tried  systematically  to  adjust  the  parameters 
to  obtain  a  good  repeatability  in  layer  production.  Concerning  IAD,  we  show  that 

the  quality  of  vacuum  is  a  preponderant  parameter,  by  comparing  results  obtained 
on  layers  realized  with  cryogenic  pumping  and  with  classical  diffusion  pump.  The 

results  are  given  in  Table  II. 

-  the  refractive  index  is  represented  as  the  development 
B  C 

n  =  A  +  +  ^4"  where  X  is  expressed  in  pm 


For  an  easy  comparison,  we  give  also  the  value  of  n  obtained  for  X  =  0,63  pm . 

Inhomogeneity  defects  are  very  low.  With  our  notations,  An  represents  the 
global  variation  in  the  volume  of  the  layer  and  the  ratio  A n/n  is  also  given. 

Concerning  the  extinction  coefficients  measured  with  photothermal  methods, 
we  can  observe  that  IAD  leads  to  values  lower  than  Ion  Plating.  At  the  present  time 
we  make  some  attempts  to  see  in  what  extent  it  is  possible  to  reduce  absorption  losses 
with  this  Ion  Plating  technique. 
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The  grain  size  measured  on  I  AD  layers  is  interesting  because  it  is  much  lower 
than  in  Ti02  layers  obtained  by  classical  evaporation.  On  this  subject  we  have  not 
yet  got  any  result  with  Ion  Plating. 

To  conclude,  for  Ti02  layers,  we  have  shown  that  both  IAD  and  Reactive  Ion 
Plating  exhibit  numerous  advantages  compared  with  classical  evaporation:  first, 
spontaneous  moisture  adsorption  is  complety  eliminated.  Furthermore,  the 
repeatability  of  the  index  value  is  easily  obtained  with  the  additional  advantage  of  a 
negligible  defect  of  homogeneity. 

A  detailed  comparison  of  IAD  and  RIP  methods  must  also  concern  the  problems 
of  uniformity  of  deposition.  As  soon  as  great  dimensioned  surfaces  are  involved 
(several  tens  of  centimeters  in  diameter)  Ion  Plating  seems  to  be  best  fitted.  With  IAD 
we  are  limited  by  the  dimensions  of  the  ion  beam.  It  can  be  deduced  from  this  study 
that,  to  obtain  minimal  losses,  the  IAD  technique  with  a  cryogenic  pumping  can  lead 
to  the  best  results.  Obviously,  we  cannot  be  sure  that  the  performances  given  here 
for  RIP  techniques  cannot  be  improved. 

It  would  be  interesting  to  develop  all  these  studies  to  other  materials  to 
produce  optical  layers,  but  it  is  impossible  to  perform  all  this  work  without  an 
extensive  cooperation  between  several  laboratories.  The  exchange  of  results 
between  Tucson  and  Marseille  is  very  useful  to  help  this  study  of  the  techniques  of 
deposition. 

MEASUREMENTS  OF  ABSORPTION  LOSSES  BY  COLLINEAR  PHOTOTHERMAL 
DEFLECTION  TECHNIQUE 

Our  technique  of  measurement  of  absorption  losses  in  a  layer  by  photothermal 
deflection  spectroscopy  has  been  improved  in  order  to  have  a  good  calibration  of  the 
measurements.  From  them  by  calculation  we  can  deduce  the  value  of  the  extinction 
coefficient  k,  knowing  the  value  of  the  real  part  of  the  refractive  index  and 
thickness  (determined  both  by  spectrophotometric  measurements). 

In  this  report  we  will  not  recall  the  detail  of  the  principle  of  the  collinear 
photothermal  apparatus  we  used;  it  has  been  described  elsewhere  in  the  literature. 

Extinction  coefficients  of  titanium  dioxide  films  prepared  using 

different  techniques. 

We  have  applied  this  collinear  photothermal  deflection  method  to  the 
determination  of  extinction  coefficients  of  single  layer  Ti02  films.  The  results 

presented  concern  the  samples  prepared  for  the  1986  OSA  annual  meeting.  All  these 

Ti0  2  layers  have  been  deposited  by  different  techniques  on  identical  silica 
substrates.  Table  III  gives  the  list  of  the  19  samples  we  have  measured,  with  the 
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Number 

Technique 

Supplier 

■EI+EBKHi 

E.B. 

Spectra  Physics,  CA 

E.B. 

Spectra  Physics,  CA 

E.B. 

Leybold  Heraeus,  Hanau 

138 

E.B. 

Optical  Sciences  Center,  AZ 

052 

E.  B. 

Indian  Inst,  of  Science 

094 

E.  B. 

Martin  Marietta  Aerospace,  NM 

036 

E.B. 

Opto  Mechanik,  FL 

024 

I.A.D. 

Leybold  Heraeus,  Hanau 

137 

I.A.D. 

Optical  Sciences  Center,  AZ 

143 

I.A.D. 

Optical  Sciences  Center,  AZ 

013 

I.A.D. 

Univ.  of  New  Mexico 

Oil 

I.A.D. 

Univ.  of  New  Mexico 

135 

I.B.S.D. 

Rockwell,  CA 

078 

I.B.S.D. 

Colorado  State,  Univ.  CO 

120 

I.B.S.D. 

O.C.L.I. 

057 

A.R.E. 

Indian  Inst,  of  Science 

093 

R.F.S. 

Battelle  -  N.W.L.,  WA 

165 

I.P. 

Balzers 

114 

I.P. 

O.C.L.I. 

026 

bare  substrate 

124 

bare  substrate 

125 

bare  substrate 

E.B. 

Electrom  Beam  Evaporation 

I.A.D. 

Ion  Assisted  Deposition 

I.B.S.D. 

Ion  Beam  Sputter  Deposition 

AJLE. 

Activated  Reactive  Evaporation 

R.F.S. 

RF  Diode  Sputtering 

l.P. 

Ion  (or  Plasma)  Plating 

Table  III: 

Suppliers  of  Titania  films 

(OSA  meeting  1989  -  published  in:  "Comparison  of  the  properties  of  Titanium  dioxide  films 
prepared  using  various  techniques,"  J.M.  Bennett,  E.  Pelletier,  G.  Albrand,  J.P.  Borgogno,  B. 
Lazaridfcs,  C.K.  Camiglia,  R.A  Schmell;  T.H.  Allen,  T.  Tuttle-Hart,  K.H.  Guenther,  A.  Saxer, 
Applied  Optics,  28,  3303-3317  (1989) 


preparation  methods  used  and  suppliers  name.  Table  IV  gives,  for  each  sample,  the 
values  of  refractive  index  a  for  X  *  600  nm  and  of  thickness  e,  values  determined  at 
ENSPM  by  spectrophotometry,  together  with  the  absorption  losses  A  measured  by  the 
pbotothermal  method.  Finally,  the  values  of  the  extinction  coefficient  k  of  the 
material  constituting  the  layer  calculated  from  those  of  A,  n  and  e  are  given.  We  can 
state  that  the  k  values  determined  by  the  pbotothermal  technique  are  consistent 
with  those  determined  at  ENSPM  by  spectrophotometry  (not  given  here). 
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Sample 

N° 

Technique 

n* 

(X*600nm) 

e(nm)* 

A(t-0)x10s" 
(l-30°-Pol  S) 

k  xIO5” 

extinction 

coef. 

AA 

A 

m 

EB 

2,292 

308 

22 

4,1 

-0,14 

2,279 

315 

53 

1  0 

-0,05 

WMM 

2,338 

345 

122 

23 

-0,04 

138 

2,280 

310 

60 

1  1 

+.0,48 

052 

2,181 

386 

30 

4,4 

-0,12 

094 

2,222 

333 

43 

8 

-0,05 

036 

2,330 

299 

83 

1  5 

-0,17 

024 

1AD 

2,472 

300 

2 

0,4 

0 

137 

2,401 

283 

30 

5,5 

+  0,55 

143 

2,393 

267 

200 

33 

+  0,09 

013 

2,318 

238 

52 

1  2 

-0,15 

01 1 

2,438 

314 

70 

1  6 

-0,22 

135 

IBSD 

2,446 

285 

24 

4,8 

0 

078 

2,386 

209 

120 

39 

-0,04 

120 

2,508 

266 

98 

1  8 

0 

057 

ARE 

2,172 

374 

32 

5,3 

-0,03 

093 

RFS 

2,500 

264 

28 

5,1 

0 

165 

IP 

2,483 

272 

20 

3,6 

0 

114 

2,575 

261 

21 

4.1 

0 

Table  IV: 

OS  A  samples.  Measurements  of  n,  e,  A  performed  in  Mareilles  by. 

♦Spectrophototnetric  measurements  (6  -  7) 

**  Photothermal  deflection  measurement. 


Absorption  has  been  measured  at  different  points  of  the  sample  surface.  The 
uniformity  is  rather  good  on  a  diameter  of  5  mm  since  the  relative  deviations  are 
smaller  than  20%;  the  only  exception  is  sample  078  (1BSD)  which  possesses  a  very 

poor  uniformity  (deviations  of  more  than  40%). 

The  first  question  we  can  ask  is:  does  the  value  of  k  depend  on  the  deposition 

technique  used?  The  results  concerning  the  statistics  of  the  values  of  the  extinction 

coefficient  are  given  in  the  form  of  4  histograms  (fig.  19).  The  values  of  the 
extinction  coefficient  are  given  on  the  abscissa  in  logarithmic  scale.  Each  of  the 
histograms  a,  b,  c,  d,  is  related  to  a  deposition  technique:  respectively  Electron  Beam 
samples  (2-a),  Ion  Assisted  Deposition  samples  (2-b),  Ion  Beam  Sputter  Deposition 
samples  (2-c)  and  samples  from  other  techniques  (2-d)  are  pictured  as  striped 
columns.  In  each  histogram  the  k  values  for  all  the  other  samples  (dotted  columns) 

are  also  given. 
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Figure  19: 

Histograms  of  extinction  coefficients.  Each  figure  is  related  to  a  deposition 
technique  -  a:  EB  samples;  b;  IAD  samples;  c;  [BSD  samples;  d;  other  samples  (ARE,  RFS, 
IP) 
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If  we  take  a  mean  value  of  the  extinction  coefficient  on  the  whole  set  of 
samples,  we  obtain  1.2  10‘4. 

-  we  find  the  same  mean  value  for  the  EB  samples  (1.1  10'4)  and  for  the  IAD 
samples  (1.3  10‘4). 

-  this  mean  value  is  higher  for  the  IBSD  samples  (2.1  10'4).  However  this  set 
consists  of  only  three  samples.  Among  them  the  n°  078  which  is  non  uniform  and 
highly  absorbing  (k  =  4  xlO*4).  This  leads  us  to  ignore  it 

-  the  mean  value  of  k  is  lower  (0.4  10'4)  for  the  four  samples  of  the  last  group, 
corresponding  to  other  different  techniques  (  ARE,  RFS,  IP). 

We  must  note  that,  in  each  set  of  samples,  the  standard  deviations  are  very 
high  except  for  the  last  one  (  ARE,  RFS,  IP).  These  very  high  values  seem  to  be  a 
consequence  of  the  fact  that  the  samples  have  been  produced  by  different 
laboratories,  with  different  deposition  conditions. 

For  all  these  19  samples  measured  there  is  no  obvious  correlation  between  the 
value  of  the  extinction  coefficient  of  the  deposited  material  and  the  technique  of 
film  deposition  .  EB,  IAD  or  IBSD  give  high  as  well  as  low  absorbing  layers.  RFS  and 
IP  samples  have  somewhat  lower  values  of  k  but  we  don't  have  enough  samples  of 
these  techniques  to  reach  any  conclusion. 

The  refractive  index  n  is  plotted  versus  k  in  figure  20a  (  for  all  samples 
(Table  III).  We  already  know  that  the  n  values  are  related  to  the  deposition 
technique:  EB  samples  have  the  lowest  n  values,  IAD  samples  have  higher  n  values 
and  IBSD  samples  or  IP  samples  have  the  highest  n  values.  Furthermore,  from  Figure 
20a  we  can  see  that  the  five  layers  with  both  high  refractive  index  and  low 
extinction  coefficient  (upper  left  part  of  the  figure  3)  arc  produced  by  IP  (2 

samples/2),  RFS  (1  sample/1),  IBSD  (1  sample/3)  and  even  IAD  (1  sample/5).  If  we 

consider  only  the  mean  values  of  n  and  k,  for  each  technique,  we  find  that  IP  and 
RFS  produce  titania  films  with  both  the  highest  n  value  and  the  lowest  k  value. 
Further  experiments  with  more  samples  prepared  by  each  technique  seem  necessary 

to  confirm  these  results. 

Photo-induced  absorption  change  In  titania  layers 
Evolution  of  illuminated  layers 

The  values  of  absorption  coefficient  given  in  Table  II  are  initial  values:  the 

first  photothermal  measurement  is  performed  about  1  minute  after  beginning 
illumination  of  the  sample  (the  measurement  itself  takes  a  few  seconds).  For  each 
sample,  the  photothermal  signal  is  recorded  over  at  least  half  an  hour. 

When  the  samples  are  illuminated  by  the  dye  laser  for  several  hours,  an 
evolution  in  the  photothermal  signal  is  often  observed.  It  is  the  case  for  13  of  the  19 
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Figure  20: 

Refractive  index  versus  extinction  coefficient.  Each  sample  is  pictured  as  a 
square.  In  figure  (a),  the  deposition  technique  is  indicated  inside  the  square.  In 
figure  (b),  the  A  A/A  values,  expressed  as  a  percentage,  are  given  inside  the  squares. 
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Ti02  samples  from  OSA.  It  is  strange  to  observe  that  for  ten  of  them  absorption 
decreases  with  time  while  it  increases  for  the  three  other  samples.  The  six 
remaining  samples  were  perfectly  stable  during  laser  exposure,  even  after  several 
hours. 

It  is  interesting  to  study  more  carefully  this  phenomenon  and  to  display  the 
law  of  absorption  variation  with  time  t.  In  fig.  21  are  plotted  the  curves  giving 
A(t)/A(t=0)  for  some  of  the  OSA  samples  compared  with  samples  from  our  laboratory, 
a  -  Ti02  n°  11  -  OSA  (Univ.  of  New  Mexico  -  IAD) 
b  -  Ti02  (Marseille  -  IAD) 
c  -  Ti02  n°  137  -  OSA  (OSC  -  IAD) 
d  -  Ti02  n°  120  OSA  (OCLI  -  IBSD) 

If  we  consider  the  curve  21a,  we  observe  a  rapid  decrease  of  absorption  during 
the  first  ten  minutes  of  illumination,  then  the  change  of  absorption  becomes  nearly 
constant  with  time.  Some  films  manufactured  in  Marseille  present  the  same 
evolution  (curve  21b). 

Curve  21c  gives  the  other  kind  of  evolution  of  absorption;  it  corresponds  to 
one  of  the  three  OSA  samples  in  which  the  absorption  increases  with  exposure  time 
(Ti02  n°  137,  n°  138  and  n°  143,  all  of  them  produced  at  OSC).  It  is  interesting  to  note 
the  effect  of  annealing.  These  three  samples  were  heated  for  48  hours  at  200°C.  After 
this,  the  evolution  of  A(t)  is  completely  different,  since  we  now  observe  a  decrease 
versus  time.  Thus  after  annealing,  the  3  samples  have  reached  a  "standard"  behavior 
that  is  similar  to  that  of  most  of  T1O2  layers. 

Curve  21d  concerns  a  Ti02  sample  whose  stability  is  excellent. 

We  can  verify  that  absorption  has  changed  only  in  the  illuminated  area  of  the 
surface;  if  the  sample  is  moved,  the  same  phenomenon  is  observed  at  this  other 
point.  We  have  also  noticed  that  the  relative  variation  of  absorption  does  not  depend 
on  the  power  of  the  exciting  beam,  at  least  in  the  range  30  mW  -  300  mW 

Let  us  now  give  a  brief  analysis  of  the  reasons  of  this  evolution.  The  power 
density  received  by  the  sample  being  rather  high  (about  400  W/cm2),  we  could  think 
that  it  gives  rise  to  a  local  overheating.  This  is  why  we  have  observed  the  surface 
with  a  high  sensitivity  infrared  camera.  No  temperature  increase  greater  than  0.2°C 
(which  is  the  thermal  resolution  of  the  camera)  can  be  detected  for  samples  with 
absorption  coefficient  lower  than  5  10'2,  The  phenomenon  of  absorption  change 
observed  is  thus  not  connected  with  an  increase  in  temperature,  but  is  a 
photoinduced  effect.  The  photothermal  deflection  method  permits  the  study,  in  real 
time,  of  the  problems  of  instability  of  illuminated  materials  in  thin  film  form. 

We  can  wonder  to  what  extent  there  is  a  correlation  between  temporal 
instability  of  Ti02  layers  and  other  characteristic  parameters;  refractive  index  n, 
extinction  coefficient  k,  and  obviously,  the  deposition  technique.  For  that  we  have 
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Figure  21: 

Relative  variations  of  absorption  versus  time  for  4  illuminated  single  layer 
titania  films,  a:  n°  11  (OSA);  b:  TiOj  (ENSPM);  c:  n°  137  (OSA);  d:  n°  120  (OS A). 


given  in  the  last  column  of  table  IV  the  algebraic  relative  variation  after  ten 

minutes  of  illumination: 

AA  A(t=10mn)  -  A(t=0) 

A  "  A(l=0) 

Figure  20b  gives  a  cartography  of  the  samples  in  the  n,  k  plane.  The  figures, 
inside  the  square  representing  each  sample,  correspond  to  expressed  as  a 

percentage.  Tbe  6  stable  samples  are  marked  with  a  0.  It  clearly  appears  that  the 
only  stable  samples  are  those  whose  index  is  greater  than  2.45.  Concerning  the 
extinction  coefficient,  the  situation  is  more  complex;  we  can  find  both  unstable 
samples  with  low  k  and  very  stable  samples  with  high  extinction  coefficients. 
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However,  in  figure  3  it  is  obvious  that  five  of  the  six  stable  samples  are  grouped 
together  in  the  upper  left  part  of  the  figure,  that  is  they  have  both  high 
refractive  index  and  low  extinction  coefficient.  On  the  other  hand,  it  is 
interesting  to  note  that  these  six  stable  samples  correspond  to  the  following 

techniques:  2  IP  (out  of  2),  1  RFS  (out  of  1),  2  IBS  (out  of  3),  1  IAD  (out  of  5).  None  of 

the  7  samples  obtained  by  EB  is  stable  when  illuminated. 

Experiments  carried  out  in  our  laboratory  show  that  it  is  very  difficult  to 
deposit  Ti02  films  with  high  stability,  even  by  IAD.  We  have  observed  that  IAD  films 
have  generally  a  better  stability  than  EB  films.  However,  IAD  technique  allows  a 
systematic  fabrication  of  stable  Ta205  layers. 

Intermittent  illumination 

As  indicated  above,  for  most  of  single  Ti02  layers,  absorption  under 

illumination  decreases.  When  illumination  is  stopped,  we  can  observe  a  partial 
reversibility  of  the  phenomenon.  An  example  (sample  11)  is  given  in  figure  22. 
Curve  a  corresponds  to  a  first  illumination  of  20  min.  Between  curves  a  and  b.  the 
sample  was  unilluminated  for  20  hours,  in  the  same  conditions  of  temperature  and 
moisture  (T  =  20°C;  60%  moisture).  During  this  time,  absorption  has  partially 
recovered.  The  same  experiment  can  be  repeated  several  times  with  similar 
evolution  under  the  laser  beam,  followed  by  a  partial  reversibility.  The 
unillumination  times  are:  24  hours  between  curves  b  and  c,  then  96  hours  between  c 
and  d  and  one  hour  between  d  and  e.  The  regenerating  times  are  thus  rather  large 
with  respect  to  the  illumination  times. 

Figure  22  shows  that: 

-  the  mean  evolution  rate  is  slower  out  of  the  beam  than  under  the  beam; 

-  these  two  rates  are  high  at  the  beginning  and  then  rapidly  stabilize. 

On  the  other  hand,  the  reversible  part  of  the  absorption  change  depends  in  a 
great  part  on  the  sample,  but  we  have  not  found,  for  this  set  of  samples,  an  obvious 
relation  between  instability  under  illumination  and  reversibility. 

Influence  of  the  hygrometric  degree  of  atmosphere 

It  is  well  known  that  the  moisture  content  of  air  modifies  optical  properties  of 
thin  layers.  We  will  show  now  that  it  modifies  also  their  stability  under  and  out  of 
the  flux.  For  that  a  single  layer  TiOj  film  has  been  submitted  to  the  following  cycle 
(figure  23): 

moisture  degree  60  % _ ^ _ (point  A)  95  % _ _ 60  % _ 

temperature  _  20°C  _ 
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Figure  23  shows  first  that  an  increase  in  moisture  leads  to  a  decrease  in  absorption 
and  that  this  phenomenon  is  reversible.  Then  it  appears  that  water  saturation 
prevents  the  regeneration  (point  A  indicated  in  the  figure  6).  We  have  observed 
these  two  phenomena  on  other  samples  of  a  TiC>2  layer. 

Conclusion 

We  have  used  the  very  sensitive  collinear  photothennal  deflection  method  to 
measure  the  extinction  coefficients  k(X)  of  materials  deposited  in  thin  layers  on 
transparent  substrates.  The  limit  of  detectivity  of  our  apparatus  corresponds,  for  a 
Ti02  layer  on  fused  silica  substrate,  to  an  extinction  coefficient  k  of  some  10-7. 

Measurements  made  on  the  19  Ti02  films  manufactured  by  different 
techniques  (EB,  IAD,  IBS,  RFS,  IP)  do  not  display  an  obvious  correlation  between  the 
value  of  k  (extreme  values:  4  xl0‘®  and  390  x  10‘®  for  X  =  600  nm)  and  the  deposition 
technique.  However,  if  we  consider  both  n  and  k,  we  find  that  high  values  of  n  and 
low  losses  can  be  obtained  with  the  RFS,  IP,  IBSD,  and  even  IAD  techniques. 

In  addition  we  demonstrate  that  the  PDS  method  is  a  good  means  to  observe  a 
photoinduced  absorption  change  in  some  of  the  layers  placed  in  a  laser  beam.  Most 
of  the  TiC>2  layers  have  decreasing  absorption  during  half  an  hour  of  exposure  to  a 
beam  of  irradiance  of  400  W/cm2. 

The  relative  change  of  absorption  depends  in  a  great  part  on  the  sample;  it 
can  be  very  high  (55%  after  ten  minutes  of  illumination),  but  it  is  null  for  the  six 
higher  refractive  index  samples.  The  IP,  RFS  and  IBSD  techniques  seem  to  be  the 
methods  of  layer  production  that  lead  to  the  highest  index  values  and  to  the  best 
stability  under  flux. 
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Measurement  of  scattering  curves:  study  of  the  correlation 

STATE  BETWEEN  INTERFACES  IN  A  LAYER  STACK 

Studies  of  light  scattering  from  optical  coatings  are  a  good  mean  to 
characterize  materials  in  thin  film  form.  We  have  at  our  disposal  in  Marseilles 
experiments  and  theoretical  methods  that  are  both  very  powerful.  Our  main  attempt 

was  to  show  how  we  can  obtain  results  on  layer  microstructure  by  using 
simultaneously  theory  and  experiment.  It  is  of  interest  to  review  here  the  main  steps 
of  the  reasoning  leading  to  these  results.  Calculations  are  not  given  here.  More 
information  about  this  can  be  found  in  the  published  papers  or  thesis  (in  French) 

from  C.  Amra  and  P.  Roche. 

We  know  that,  in  multilayer  coatings,  the  relation  giving  the  energy  balance 
is:  1=R  +  T  +  A  +  D 

where  A  and  D  represent  respectively  absorption  and  scattering  losses. 

These  losses  A  and  D  are  very  often  extremely  low  and  do  not  exceed  some 
thousandths.  Nevertheless  these  two  values  must  be  perfectly  controlled  for  high 
performance  optical  coatings.  For  this,  we  must  go  further  in  the  details  of  this 

equation  and  take  account  of  external  scattering  directly  accessible  by 
measurement  and  of  internal  scattering  which  is  practically  added  to  intrinsic 
absorption  of  materials  in  thin  films.  With  our  notations,  D  represents  the  integral 
in  all  the  space  of  the  amount  of  light  scattered  in  each  direction:  it  is  the  external 
global  scattering. 

Many  work  has  been  devoted  to  the  study  of  scattering  and  particularly  on  the 

the  law  of  dependence  of  D  with  wavelength.  A  simple  example  is  given  here  in 

Figure  24.  It  concerns  a  single  layer  of  optical  thickness  Xo/2  and  of  index  nc  =  2.3. 
In  this  calculation;  the  roughnesses  of  interfaces  subtrate-layer  and  layer-air  are 
assumed  to  be  identical  and  a  represents  the  correlation  between  interfaces.  We 

have  considered  here  only  the  two  extreme  cases:  a  =  0  and  a  =  1  which  correspond 
respectively  to  a  perfect  correlation  and  to  a  complete  uncorrelation.  Compared  with 

scattering  of  the  bare  substrate,  it  can  be  seen  that  the  law  of  variation  of  D  versus  X 

depends  considerably  on  o. 

Numerous  studies  deal  with  scattering  in  multilayer  coatings  but,  due  to  the 
high  number  of  parameters  involved,  the  conclusions  are  very  difficult  since  there 
are  always  many  solutions  to  explain  the  experimental  results. 

Even  in  the  case  of  a  single  layer,  if  we  have  acces  to  the  measurement  of 
angular  distribution  of  scattering  between  0  and  90°,  one  can  explain  the  value  of 
scattering  by  reflection  by  considering  one  or  the  other  of  the  two  extreme 
hypothesis  of  correlated  roughness  (a  =  1)  or  uncorrelated  roughness  (a  =  0). 
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Figure  24: 

External  global  scattering  versus  X.  The  results  of  the  calculation  are  given 
for  a  X/2  layer  of  tic  =  2.30.  The  roughnesses  are  identical  for  the  two  cases: 
correlated  interfaces  a  =  1  and  uncorrelated  interfaces  a  =  0. 


It  is  obvious  that  this  ambiguity  will  be  vety  difficult  to  solve  in  the  case  of  a  stack  of 
many  layers  since  we  can  vary  independently  the  roughnesses  of  different 
interfaces  and  the  states  of  correlation.  We  will  see  here  how  the  use  of  a  causal 
model  together  with  the  measurements  of  anisotropy  variation  of  scattering  permits 
an  easy  investigation  of  microroughness  and  leads  to  obvious  conclusions  to  solve 
this  problem. 

We  want  to  give  a  response  to  this  question:  what  are  the  key  parameters  of 
scattering  and  what  information  do  they  give  us  on  material  microstructure  at  the 
scale  of  optical  frequencies? 

The  apparatus  intended  for  the  plotting  of  scattering  curve  for  the 
wavelength  X  =  0.6328  pm  has  been  decribed  elsewhere.  In  the  same  way.  the  theory 
has  been  given  in  the  literature  (Elson,  Amra,  Bousquet).. 

With  this  vector  theory,  we  can  write  that  the  amount  of  light  scattered  in  the 
direction  (0,  0)  is  given  by  the  expression: 

P 

1(0.  B)  =  I  Cij  (0,0)  Oij  (9.0)  fjj  (0,0)  ( 1 ) 

i.j 


where 
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-  Cij  are  called  ideal  coefficients  since  they  depend  only  on  the  perfect  stack 
supposed  free  of  roughness:  they  depend  on  the  coating  design  and  on  the 

illumination  and  observation  conditions. 

-  Yjj  is  the  roughness  spectrum  of  the  interface  j  and  contains  all  the 

information  relative  to  the  surface  defects,  it  is  the  Fourier  Transform  of  the 
autocorrelation  function  of  the  interface  profile. 

-  ay  is  the  correlation  between  the  interfaces  i  and  j.  It  can  be  considered  as 
the  degree  of  coherence  between  the  scattering  sources  at  the  interfaces  i  and  j  or  as 
the  degree  of  similarity  of  the  surface  profiles  i  and  j.  Usually,  the  roughness 
spectrum  is  assumed  to  be  the  sum  of  the  Hankel  Transform  of  an  exponential  and  a 
gaussian  function;  so  4  parameters  are  used  for  each  spectrum.  For  a  multilayer,  this 
leads  to  many  parameters:  5p+4  parameters  for  a  p-layer  stack. 

Certainly  the  key  parameter  of  scattering  is  the  correlation  between 
interfaces.  Indeed,  if  we  have  uncorrelated  surfaces,  scattering  is  always  increasing 
with  an  increasing  number  of  layers  and  decreasing  in  the  case  of  identical 
surfaces.  We  will  have  interferences  between  the  scattering  waves,  and  providing 
that  these  interferences  are  destructive,  we  can  reduce  the  scattering. 

This  phenomenon  has  been  observed  on  a  single  layer  and  described  as 
"antiscattering  effect"  and  it  is  also  the  case  of  Fabry-Perot  filters  since  we  know 
that  scattering  can  vary  by  a  factor  100  according  to  the  state  of  correlation  between 
interfaces. 

Usually,  the  correlation  is  assumed  to  be  a  constant  value  between  zero  and 
unity.  In  fact,  we  must  consider  a  more  realistic  model  ;  for  instance,  we  must  admit 
that  the  surface  defect  at  interface  j  is  due  to  the  total  or  partial  reproduction  of  the 
defects  of  interface  i  added  with  residual  roughness  due  to  the  grain  size  that 
characterizes  the  microstructure  of  the  material.  First  we  will  consider  here  the 
model  where  the  Fourier  Transform  of  residual  roughness  is  negligible  in  the  range 
of  spatial  frequencies  responsible  for  optical  scattering,  since  we  know  that 
columnar  microstructure  of  materials  is  not  responsible  for  optical  scattering. 
Anyway,  we  will  first  consider  coatings  produced  on  very  rough  substrates  to  be 
sure  that  the  residual  roughness  will  be  negligible  compared  with  the  reproduction 
of  substrate  defects.  Therefore  the  consequence  is  that  the  material  acts  as  a  linear 
filter  on  the  interface  and  the  correlation  is  the  Fourier  Transform  of  a  filtering 
function  characterizing  the  action  of  the  material.  The  advantage  of  this  model  is 
that  we  have  only  few  parameters  since  cross-correlation  laws  and  roughness 
spectra  are  connected  by  the  following  relationship:  Yii  =  Yjj  I  <*ijl 2 

The  scattered  intensity  is:  I  (9)  =  X  Cy  ay  Yjj 
So,  instead  of  (5p+4),  we  now  need  only  two  parameters  to  predict  scattering 
introduced  by  the  multilayer.  The  roughness  of  the  substrate  is  given  by  4 
parameters  and  can  be  measured  before  coating. 
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C  ij  is  given  by  theory  (connected  with  the  design) 

Tjj  is  the  interface  roughness  and  is  deduced  from  the  roughness  spectrum  of 
the  substrate  if  the  correlation  law  is  given:  Yjj  =  f  (Ys) 

Ys  must  be  measured  before  coating  and  the  correlation  between  i  and  j  is  given  as  a 
function  of  two  filtering  functions  that  characterize  the  two  materials  H  and  L  of  the 
stack:  a*j  =  f  (ay+i,  ....  Oj-ij) 

With  two  alternate  materials,  we  write: 

<*2k,2k+l  =  OH(o) 

and 

«2k+l,2k+2  =  <*L(®) 

vihich  are  the  two  functions  to  be  determined  for  prediction  of  scattering. 

Moreover,  if  we  consider  the  material  microstructure,  it  is  obvious  that  large 
period  defects  are  perfectly  reproduced  (o  (a)  =  1)  while  small  period  defects  are 
only  partially  reproduced.  Practically,  we  can  consider  that  a(o)  will  behave  as  a 
low-pass  filter. 

Briefly,  provided  that  the  substrate  was  measured  before  coating,  prediction  of 
scattering  for  a  p-layer  stack  requires  only  the  knowledge  of  two  cut-off 
frequencies  oh  and  ol-  Figure  25  concerns  a  layer  of  index  n  =  2.3:  it  shows  an 
example  calculated  of  the  scattering  curve  in  the  case  of  normal  incidence 
illumination,  for  wavelength  0.6  pm  and  for  a  series  of  cut-off  frequencies  from  oc  = 
lpm'1  to  oc  =  “• 


Figure  25: 

Layer  of  index  2.3,  optical  thickness  X/2  for  X  =  0.6  pm.  Scattering  curves  calculated 
with  different  values  of  cut  off  frequency.  oc  *  «>  corresponds  to  a  perfect 
correlation  between  the  two  interfaces  of  the  layer. 
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To  go  further  in  our  investigation,  we  can  use  the  maps  of  anisotropy  before 
and  after  coating.  It  is  possible  to  work  with  an  angular  autocorrelation  function 

and  to  demonstrate  how  the  isotropy  degree  depends  on  spatial  frequency. 

The  process  of  investigation  of  optical  microroughness  with  the  fit  between 
prediction  and  experiment  is  summarized  in  figure  26. 

In  this  presentation,  we  will  show  some  results  concerning  the  study  of  the 
mean  plane  section  of  angular  scattering  on  a  particular  design. 

With  the  Thin  Film  technical  Group  of  OSA  (meeting  on  October  16,  1989)  we 
have  presented  the  results  of  an  experiment  that  gathers  works  done  by  several 
laboratories. 

According  to  the  theory  of  scattering  in  multilayers  (Claude  Amra  - 

Marseille),  one  can  determine  the  correlation  of  interface  roughness  in  multilayer 
by  a  simple  test.  The  test  consists  in  measuring  the  scattering  of  a  design  composed 
of  alternating  high  and  low  index  layers  at  the  wavelength  where  the  layers  are  one 
halfwave  thick.  Four  laboratories  (Bazers,  OCLI,  Spectra  Physics,  OJAI)  have 

prepared  coatings  of  design:  Glass  2H  2L  2H  for  633  nm  and  have  sent  them  to 
Marseille  for  examination.  Each  laboratory  used  a  different  deposition  process  or 

material  combination. 

The  curves  (figure  27)  show  how  the  scattering  curve  is  theoretically  a 

function  of  the  value  of  the  correlation  coefficient  a  betwen  the  interfaces  of  this 
stack.  If  we  compare  these  curves  with  the  angular  scattering  of  the  substrate,  we 
can  also  demonstrate  that  the  angular  scattering  of  the  substrate  before  coating  is 

very  close  from  the  angular  scattering  after  coating  if  the  surfaces  are  perfectly 
correlated  (only  a  slight  difference  at  large  angles  9). 

The  ratio  rj  of  global  scattering  after  coating  to  global  scattering  before 
coating  is,  for  a  roughness  of  0.5  nm. 

n  =  ]'?  =  1.06  for  a  =  1  t\  28.2  for  a  =  0 

6.6  ppm  6.6  ppm 

Table  V  gives  the  values  D  of  scattering  losses  measured  for  all  coatings 
together  with  the  ratio  h  =  ^  (bef  ore  Soaring )'  The  asterisk  (*)  indicates  that  the 

corresponding  sample  was  kept  uncoated.  It  appears  from  this  table  that  all 
techniques  lead  to  good  correlation:  t\  is  much  lower  than  the  value  28  calculated  for 
uncorrelated  surfaces. 

The  measured  angular  scattering  for  the  4  samples  are  given  before  and  after 
coating  on  figures  28,  29,  30,  31.  We  can  see  that  the  correlation  is  not  completely 
perfect  for  the  coatings  since  the  two  curves  before  and  after  coating  are  not 
superimposed. 
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Investigation  of  optical  microroiighnessj 


-  Mean  plane  section  of  angular  scattering 
•  Angular  autocorrelation 

-  Variation  of  isotropy  with  spatial  frequency 


and  complementary  information  with  A1  layer  technique 


in 


Figure  26: 

Optimum  approach  for  progress 


study  of  optical  microroughness. 


Figure  27: 

Glass  2H  2L  2H  for  633  nm.  Scattering  diagrams  calculated  for  different  values 
of  the  correlation  coefficient  a  between  interfaces. 


Technique 

Materials 

re-cleaning 

Scattering  D(ppm) 
from  bare  subst. 

N*  1st  clean  2nd  clean 

scattering  D 
from  2H2L2H 
(ppm) 

■ 

30 

7.0 

10.2 

1.46 

Spectra-  tp 

React-evap. 

Ti02/Si02 

Y 

35 

6.0 

9.5 

1.58 

37* 

6.6 

6.0 

13* 

Bl 

Q 

OCU 

Plasma  PI. 

Zr0i/Si02 

Y 

1  5 

Ex 

35.5 

4.86 

1  7 

6.1 

25.7 

4.21 

6.5 

5.5 

Balzers 

Ion  pi. 

Ta205/Si02 

Y 

03 

6.2 

8.4 

1.35 

? 

9.8 

21 

6.9 

12.6 

1.83 

OJAI 

IBS 

TiC>2/Si02 

N 

22 

6.9 

11.0 

1.59 

Table  V: 

Results  obtained  on  the  OS  A  samples.  D  is  given  before  and  after  coating. 
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Figure  28: 

Measured  scattering  diagrams  before  and  after  coating  2H  2L  2H 
Coating  manufactured  by  Spectra  Physics  -  Global  scattering  losses  are  given  on  the 
curves 


Figure  29: 

Measured  scattering  diagrams  before  and  after  coating  2H  2L  2H 
Coating  manufactured  by  O.C.L.t.  -  Global  scattering  losses  are  given  on  the  curves 


Figure  30: 

Measured  scattering  diagrams  before  and  after  coating  2H  2L  2H 
Coating  manufactured  by  Balzers  -  Global  scattering  losses  are  given  on  the  curves 


Figure  31: 

Measured  scattering  diagrams  before  and  after  coating  2H  2L  2H 
Coating  manufactured  by  OlAl  -  Global  scattering  losses  are  given  on  the  curves 
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For  the  coatings  of  OSA  experiments,  we  can  explain  the  better  values  of  r\  by  a 
correlation  of  0.990  (for  Balzers  sample  003),  of  0.986  (for  Spectra  Physics  sample 
030),  of  0.981  (for  OJAI  sample  022)  and  of  0.885  (for  OCLI  sample  017). 

Some  complementary  experiments  have  been  performed  in  our  Laboratory  for 
a  better  understanding  of  this  problem.  We  used-  Ti02-Si02  layers  produced  by  Ion 
Assisted  Deposition.  The  different  figures  32  show  that  the  agrement  between  the 
angular  curves  before  and  after  coating  depends  strongly  of  the  roughness  of  the 
bare  substrate:  it  is  obvious  that  the  correlation  is  easy  to  obtain  with  rough 
substrates.  Moreover,  whatever  the  substrate  roughness,  we  observe  a  perfect 
correlation  at  small  angles  and  the  disagreement  is  obvious  at  large  angles. 

To  give  a  more  complete  conclusion  on  this  work,  we  have  been  induced  to 
study  the  consequences  of  thickness  errors  on  the  layers  of  the  stack:  2H  2L  2H.  We 
will  present  later  the  measurements  of  reflectance  versus  wavelength  performed  on 

the  OSA  coatings;  the  centering  on  Xq  =  0.633  |im  is  not  always  perfect.  We  try  to 
apply  synthesis  methods  for  determining  the  effective  designs  of  the  multilayers 
and  we  will  then  have  the  possibility  to  study  theoretically  the  consequences  of 
these  thickness  errors  on  angular  scattering  curve. 

This  whole  work  must  be  the  subject  of  a  publication  in  common  with 
Marseille  and  OCLI  (J.  Apfel)  and  the  coating  manufacturers.  The  main  question 
which  will  be  interesting  for  us  in  the  future  is  the  law  of  dependence  of  scattering 
with  wavelength.  This  law  depends  strongly  on  the  correlation  between  interfaces 
and  it  should  be  a  further  verification  of  our  conclusion. 

Publications  and  participation  in  conferences 

The  grant  AFOSR-88-0139  is  of  considerable  help  to  link  the  laboratory  of 

Marseilles  with  Macleod's  group  and  allowed  the  participation  to  conferences  in  the 
United  States. 

C.  Amra,  M.  Commands  and  E.  Pelletier  attended  the  OSA  Conference  in 

November  88.  2  papers  were  presented. 

C.  Amra  and  E.  Pelletier  participated  in  the  Boulder  Damage  Symposium  -  Nov. 
88  and  Nov.  89  -  2  papers  were  presented. 

C.  Amra  and  E.  Pelletier  participated  in  the  OSA  Conference  -  Nov.  89  in 

Orlando.  One  paper  and  Technical  Meeting  group. 

One  post-doc.  (C.  Hickey)  and  one  student  (R.  Sprague)  from  Tucson  spent  two 
years  in  Marseille  (1988  -  1990).  R.  Sprague’s  thesis  will  be  a  work  common  to  Tucson 
and  Marseille  (April  90). 
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Further  information 

Ion  Assisted  Deposition  and  Modeling 

Many  of  the  above  topics  will  be  found  in  the  following  dissertations: 

C.K.  Hwangbo,  "Optical  thin  films  prepared  by  ion-assisted  and  ultrasound-assisted 
deposition." 

L.  Lingg,  "Study  of  thin  films  of  the  lanthanide  trifluorides  with  emphasis  on  the 

optical  properties."  (Final  text  in  preparation). 

John  Lehan,  "Microstructural  analysis  of  thin  films  by  Rutherford  Backscattering 
Spectrometry,  electron  diffraction  and  spectrophotometry.”  (In  preparation). 
Rober  Sargent  ....  "The  modeling  of  thin  film  growth  (tentative  title)."  (In 
preparation). 

Recent  publications 

C.K.  Hwangbo,  L.J.  Lingg,  J.P.  Lehan,  H.A.  Macleod,  J.L.  Makous  and  S.Y.  Kim,  "Ion- 
assisted  deposition  of  thermally  evaporated  Ag  and  A1  films,"  Applied  Optics,  28, 
2769-2778,  (1989) 

C.K.  Hwangbo,  L.J.  Lingg,  J.P.  Lehan,  H.A.  Macleod  and  F.  Suits,  "Reactive  ion-assisted 
deposition  of  aluminum  oxynitride  thin  films,"  Applied  Optics,  28,  2779-2784 
(1989) 

J.D.  Targove,  L.J.  Lingg,  J.P.  Lehan  and  H.A.  Macleod,  "Effect  of  oxygen  incorporation 
on  the  structure  of  ion-beam-assisted  LaF3  films,"  Applied  Optics,  27,  213-215 
(1988) 

B.G.  Bovard  and  H.A.  Macleod,  "Thermal  model  of  a  switching  Fabry-Perot  etalon," 
Proc.  SPIE,  821,  187-197  (1988) 

R.B.  Sargent  and  H.A.  Macleod,  "Computer  simulation  of  the  growth  of  two-element 
films,"  Proc.  SPIE,  821,  113-119  (1988) 

R.B.  Sargent.  D.Y.  Song  and  H.A.  Macleod,  "Computer  simulation  of  substrate  defect 
propagation  in  thin  films,"  Proc.  SPIE,  821,  133-140  (1988) 

J.D.  Targove,  B.G.  Bovard,  L.J.  Lingg  and  H.A.  Macleod,  "Densification  of  aluminum 
fluoride  thin  films  by  ion-assisted  deposition,"  Thin  Solid  Films.  159.  L57-L59 
(1988) 

B. G.  Bovard  and  H.A.  Macleod,  "Nonlinear  behavior  of  optical  coatings  subjected  to 

intense  laser  irradiation,"  J.  Mod.  Opt.,  35,  1151-1168  (1988) 

J.D.  Targove  and  H.A.  Macleod,  "Verification  of  momentum  transfer  as  the  dominant 
densifying  mechanism  in  ion-assisted  deposition,”  Applied  Optics,  27,  3779-3781 
(1988) 

C. K.  Hwangbo,  "Optical  thin  films  prepared  by  ion-assisted  and  ultrasound-assisted 

deposition."  PhD,  University  of  Arizona,  1988 
P.J.  Martin,  R.P.  Netterfield  and  W.G.  Sanly,  J.  Appl.  Phys.,  55,  235  (1984) 

H.K.  Pulker,  Coatings  on  Glass,  Elsevier,  Amsterdam  (1984) 

J.P.  Borgogno.  B.  Lazaridis,  E.  Pelletier,  "Automatic  determination  of  the  optical 
constants  of  inhomogeneous  thin  films,"  Applied  Optics,  21,  4020  (1982) 

M.  Commandrd,  E.  Pelletier,  "Measurements  of  absorpttion  losses  in  Ti02  films  by 

collinear  photothermal  deflection  technique,"  Conf.  presented  at  "OSA  Annual 
Meeting",  Santa  Clara,  California  -  Oct.  31  -  Nov.  1,  1988  -  Submitted  to  Applied 
Optics. 

Photothermal  Deflection  Spectroscopy 

A.C.  Boccara  and  D.  Fournier,  W.  Jackson  and  N.M.  Amer,  "Sensitive  photothermal 
deflection  technique  for  measuring  absorption  in  optically  thin  media,"  Optics 
letters,  5.  377  (1980) 
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Annexe  I 


Some  recent  publications 


The  Journal  of  Applied  Optics  (15  July  1989  issue)  gives  9  papers  of  the 
laboratories  of  Tucson  and  Marseilles.  We  give  here  two  of  the 
articles  of  Macleod's  group. 


Ion  assisted  deposition  of  thermally  evaporated  Ag  and  A1  films 


Chang  Kwon  Hwangbo,  Linda  J.  Lingg,  John  P.  Lehan,  H.  Angus  Macleod,  John  L.  Makous,  and  Sang 
Yeol  Kim 


Optical,  electrical,  and  microstructural  effect a  of  Ar-ion  bombardment  and  Ar  incorporation  on  thermally 
evaporated  Ag  and  A1  thin  films  are  investigated.  The  results  show  that  as  the  momentum  supplied  to  the 
growing  films  by  the  bombarding  ions  per  arriving  metal  atom  increases,  the  refractive  index  at  632.8  nm 
increases  and  the  extinction  coefficient  decreases,  lattice  spacing  expands,  grain  size  decreases,  electrical 
resistivity  increases,  and  trapped  Ar  increases  slightly.  In  Ag  films,  stress  reverses  from  tensile  to  compres¬ 
sive  and  in  A1  films  compressive  stress  increases.  In  the  A1  films  the  change  in  optical  constants  can  be 
explained  by  the  variation  in  void  volume.  The  reversal  of  stress  from  tensile  to  compressive  in  Ag  films 
requires  a  threshold  level  of  momentum.  The  increase  in  electrical  resistivity  is  related  to  the  decrease  in 
grain  size  and  increase  in  trapped  Ar  in  both  types  of  film.  Many  of  these  properties  correlate  well  with  the 
momentum  transferred,  suggesting  that  the  momentum  is  an  important  physical  parameter  in  describing  the 
influence  of  ion  beam  on  growing  thin  films  and  determining  the  characteristics  of  thin  metal  films  prepared 


by  ion  assisted  deposition. 


I.  Introduction 

Ion  assisted  deposition  (IAD)  is  known  to  produce 
thin  films  with  characteristics  that  conventionally  de¬ 
posited  thin  films  do  not  have.1"1  For  example,  the 
IAD  of  dielectric  thin  films  has  shown  that  it  increases 
the  packing  density  and  refractive  index,  reduces  the 
moisture  adsorption,  and  improves  the  adhesion  in 
certain  material  dependent  conditions.1"'1  IAD  can 
also  alter  the  stoichiometry  and  crystallinity  of  dielec¬ 
tric  thin  films.1011  For  metal  thin  films,  IAD  en¬ 
hanced  adhesion  in  Au,  Ag,  and  A1  films  on  glass  sub¬ 
strates.12-15  IAD  of  Cr  thin  films  resulted  in  an 
increase  in  reflectance  and  reversal  of  stress  from  ten¬ 
sile  to  compressive  when  the  ion  flux  exceeded  a  criti¬ 
cal  value. 16  In  addition,  the  nonbulklike  optical  densi¬ 
ty,  reduction  of  optical  scatter,  and  control  of 
microstructure  in  Cu  films  under  ion  bombardment 
were  reported.17-19 


John  Makous  was  with  the  University  of  Arizona,  Physics  Depart¬ 
ment,  Tucson,  Arizona  85721  when  this  work  was  done  and  is  now 
with  U.S-  Naval  Research  Laboratory,  Washington,  DC  20375-5000. 
When  this  work  was  done  Sang  Yeol  Kim  was  with  Pennsylvania 
Stale  University,  Materials  Research  Laboratory,  University  Park, 
Pennsylvania  16802  and  is  now  with  Aju  University,  Physics  Depart¬ 
ment,  Suwon,  Korea.  The  other  authors  are  with  University  of 
Arizona,  Optical  Sciences  Center,  Tucson,  Arizona  85721. 
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The  correlation  of  optical  properties  and  micro- 
structure  of  IAD  thin  films  with  ion  beam  parameters 
is  very  important  because  it  provides  insights  for  un¬ 
derstanding  and  producing  these  films.  It  is  known 
that  the  ion  energy  and  current  density  are  two  ion 
beam  parameters  which  alter  the  microstructure  of 
naturally  growing  thin  films  and  their  optical  proper¬ 
ties.  It  is  not  clear,  however,  what  physical  parameter 
(energy,  momentum,  or  some  function  of  the  two)  of 
the  ion  beam  affects  the  growth  of  thin  films. 1  Parmi- 
giani  et  a/.20  related  the  normalized  energy,  defined  as 
the  energy  delivered  by  the  bombarding  ions  per  arriv¬ 
ing  metal  atom,  to  the  optical,  electrical,  and  micro- 
structural  changes  of  metal  films  in  dual  ion  beam 
sputtering.21  On  the  other  hand,  Windischmann22  de¬ 
rived  a  model,  based  on  the  linear  cascade  theory  of 
forward  scattering,  so  that  the  intrinsic  stress  of  ion 
beam  sputtered  films  depends  on  the  imparted  mo¬ 
mentum  rather  than  energy. 

In  this  paper,  we  investigate  the  optical,  electrical, 
and  microstructural  effects  of  Ar-ion  bombardment 
and  Ar  incorporation  on  thermally  evaporated  Ag  and 
A1  thin  films.  These  effects  and  incorporated  Ar  are 
presented  in  terms  of  the  momentum  supplied  to  the 
growing  thin  film  by  the  bombarding  ions  per  arriving 
metal  atom.  We  define  the  momentum  P  normalized 
by  the  Ar  mass  as  y\E,  where  y  is  the  (Ar  ion)-to- 
(meta!  atom)  arrival  rate  at  the  substrate  and  E  is  the 
Ar-ion  energy. 

II.  Deposition  Conditions 

Ag  was  evaporated  from  a  Mo  boat,  and  Al  was 
electron  beam  evaporated  from  an  intermetallic  cruci- 
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Fig.  1.  Experimental  configuration  of  a  Ba.2* rs  760  coaling  plant. 


ble  in  a  Balzers  760  box  coater,  as  shown  in  Fig.  1.  The 
base  pressure  was  1  x  10~°  mbar.  The  chamber  tem¬ 
perature  was  ambient,  and  its  variation  during  a  single 
layer  deposition  was  <5°C.  The  deposition  rates  of 
Ag  and  A1  films  were  0.5  and  1.0  nm/s,  respectively,  as 
measured  by  a  quartz  crystal  oscillator.  Prisms,  glass 
slides,  and  graphite  slabs  were  used  as  substrates  for 
the  various  postdeposition  measurements.  Just  prior 
to  deposition,  substrates  were  sputter  cleaned  in  vacua 
by  the  Ar-ion  beam  for  5-60  s,  depending  on  the  type  of 
substrate  used. 

We  employed  a  3-cm  aperture  Kaufman  hot-cath¬ 
ode  ion  source  with  Ar  gas  as  the  bombarding  species. 
Ar  gas  was  admitted  through  the  ion  gun  until  a  pres¬ 
sure  of  8  X  10~5  mbar  in  the  chamber  was  obtained. 
The  Ar-ion  energies  investigated  were  200,  300,  and 
500  eV,  while  the  Ar-ion  current  density  was  varied 
from  0  to  20  jxA/cm2.  These  parameters  resulted  in 
the  mass  normalized  momentum  P  being  varied  from  0 
to  0.48  veV  for  Ag  films  and  0  to  0.46  \ev  for  A1  films, 
respectively.  The  non-IAD  (conventional)  thin  film 
has  zero  P  by  definition. 


III.  Optical  Measurements 

A.  Surface  Plasmon  Resonance 

Optical  constants  at  632.8-nm  waveleng'h  and 
thickness  of  Ag  and  A1  films  were  obtained  using  the 
surface  plasmon  resonance  (SPR)  technique.23  In  the 
Kretschmann  configuration  of  the  SPR,  a  thin  metal 
film  is  deposited  on  the  hypotenuse  of  a  right  angle 
prism.24  When  the  angle  of  incidence  of  a  TM  wave 
inside  the  prism  is  larger  than  the  critical  angle,  the 
incident  wave  is  totally  reflected  and  an  evanescent 
wave  is  absorbed  by  the  metal  film.  At  a  certain 
incident  angle  and  film  thickness,  the  component  of 
the  wave  vector  of  the  electric  field  parallel  to  the  film 


plane  equals  that  of  surface  plasmon  propagating 
along  the  metal-air  interface  and  is  thus  coupled  to 
this  surface  plasmon  mode.  When  this  occurs,  the 
energy  of  the  incident  wave  is  absorbed  by  the  surface 
plasmon,  and  no  wave  is  reflected  by  the  metal  film. 
This  results  in  the  reflectance  showing  a  sharp  reso¬ 
nant  minimum  as  the  incident  angle  varies  above  the 
critical  angle.  The  resonant  angle  and  thickness 
which  yields  zero  reflectance  depend  on  the  optical 
constants  of  the  metal  film.25  Optimization  software 
employing  admittance  matching  between  the  incident 
medium  and  metal  film  was  used  to  calculate  the  reso¬ 
nant  angle  and  thickness.20  Optimum  thicknesses  of 
~50  and  15  nm  were  chosen  for  Ag  and  A1  films,  respec¬ 
tively.  From  the  measured  reflectance  with  respect  to 
the  incident  angle,  the  optical  constants  and  thickness 
were  calculated  using  a  nonlinear  least-squares  meth¬ 
od.21  Reflectance  calculated  from  the  optical  con¬ 
stants  and  thickness  found  by  the  program  was  in  good 
agreement  with  the  experimental  reflectance  profile. 
Measurements  were  performed  within  5  h  after  the 
films  were  taken  out  of  the  vacuum  chamber  because 
the  SPR  technique  is  sensitive  to  the  oxidation  of  the 
metal-air  interface.28 

Refractive  index  n  and  extinction  coefficient  k  of  Ag 
films  at  632.8  nm  are  plotted  in  Fig.  2  in  terms  of  P.  As 
P  increases,  n  increases  and  k  decreases.  Reflectance 
of  Ag  films  of  100-nm  thickness  at  normal  incidence  is 
calculated  from  these  optical  constants  and  listed  in 
Table  I.  At  the  highest  P  (0.48  \eV),  reflectance  drops 
by  2.5%  compared  to  the  non-IAD  film.  In  other 
words,  IAD  Ag  films  become  lossy  when  bombarded 
with  Ar  ions. 

The  n  and  k  of  A1  films  at  632.8  nm  shown  in  Fig.  3 
are  scattered  from  a  line.  Reflectance  of  A1  films  of 
100-nm  thickness  at  normal  incidence  is  also  listed  in 
Table  I.  When  P  is  small,  reflectance  does  not  change 
much  from  that  of  non-IAD  film.  However,  when  P 
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exceeds  0.23  the  reflectance  of  IAD  films  de¬ 
creases. 


0.03 


0.1 


0.2  0J  0.4  0.3 

Momentum  (•vu*) 


Fig.  2.  Refractive  index  and  extinction  coefficient  at  632.8  nm  of 
Ag  films  measured  by  the  SPR  plotted  as  a  function  of  P. 


Fig.  3.  Refractive  index  and  extinction  coefficient  at  632.8  nm  of  A1 
films  measured  by  the  SPR  plotted  as  a  function  of  P. 


Table  I.  Reflectance  of  Ag  and  Al  Films  of  100-nm  Thickness  at  632.8 
nm  at  Normal  Incidence  Is  Caiculaled  from  (he  Optical  Constants  Obtained 
_  from  the  SPR  Measurements 


Ay  films 

P(.5V) 

/mi 

Al  films 

P  (vFV) 

/?«> 

0 

98.3 

0 

88.8 

0.18 

97.1 

0.06 

83.7 

0.22 

97.0 

0.23 

88.6 

0.29 

96.7 

0.36 

87.2 

0.20 

96.7 

0.16 

84.6 

0-27 

96.5 

0.18 

95.8 

B.  Spectroscopic  Eilipsometry 
Recently,  spectroscopic  eilipsometry  (SE)  has  been 
used  as  a  nondestructive  tool  to  depth  profile  multilay¬ 
er  structures  and  void  fractions  in  thin  films  in  addi¬ 
tion  to  its  traditional  iole,  measuring  the  optical  con- 
stants,  dispersion,  and  thickness  of  thin  films.  The 
void  fraction  in  Au  and  amorphous  IAD  Ge  thin  films 
and  depth  profiling  Si,  ZnS,  MgO,  and  ThF.,  thin  films 
have  recently  been  studied.29'33  In  this  work  SE  was 
used  to  measure  the  optical  constants,  dispersion,  and 
thickness  of  metal  thin  films  on  glass  substrates  at  an 
angle  of  incidence  of  70°  in  the  1.5-4.5-eV  spectral 
range  in  steps  of  0.05  eV.  The  void  fraction  in  Al  films 
and  its  oxide  layer  thickness  were  calculated  by  apply¬ 
ing  the  Bruggeman  effective  medium  approximation 
and  linear  regressional  analysis  to  the  SE  data.34-35 
Pseudorefractive  index  np  and  pseudoextinction  co¬ 
efficient  kp,  the  effective  refractive  index  and  effective 
extinction  coefficient  of.  the  film-substrate  assembly, 
for  some  Ag  films  are  plotted  in  Figs.  4(a)  and  (b)  in 
terms  of  photon  energy.  np  of  the  IAD  Ag  films  is 
larger  than  that  of  the  non-IAD  Ag  film  below  the 
plasma  frequency  near  3.8  eV  (at  which  np  equals  k 
But  in  the  2.5-3.0-eV  spectral  region,  n„  of  the  IAD 
films  is  almost  the  same  as  the  non-IAD  film.  In  the 
interband  transitions  above  the  plasma  frequency,  n(, 
of  the  IAD  films  is  lower  than  that  of  the  non-IAD  film. 
kp  of  the  IAD  Ag  films  is  smaller  than  that  of  the  non- 
IAD  film  except  in  the  vicinity  of  the  plasma  frequen¬ 
cy.  Also,  the  high  P  film  shows  larger  nn  and  smaller 
ft,,  than  the  low  P  film,  which  is  consistent  with  the 
SPR  results.  Reflectance  of  Ag  films  at  normal  inci¬ 
dence  is  calculated  from  np  and  kp  and  plotted  in  Fig.  5. 
Reflectance  of  these  IAD  films  drops  by  ~1%  in  the 
intraband  transitions,  while  it  drops  by  ~3%  in  the 
interband  transitions.  Around  the  plasma  frequency, 
surprisingly,  the  reflectance  (interband  absorption) 
increases  (decreases)  by  ~3%  for  IAD  films. 

np  and  kp  of  Al  films  are  plotted  in  Figs.  6(a)  and  (b). 
Near  the  1.5-eV  interband  absorption,  np  of  the  I  AD  Al 
films  is  smaller  than  that  of  the  non-IAD  Al  film.  In 
the  rest  of  the  measured  spectral  region  (from  ~1.S  to 
4.5  eV),  the  high  P  film  shows  larger  np  than  the  non- 
IAD  film,  while  the  small  P  film  behaves  like  the  non- 
IAD  film.  kp  of  the  IAD  films  is  smaller  than  that  of 
the  non-IAD  film  in  the  whole  measured  spectral  re¬ 
gion.  Reflectance  calculated  from  np  and  kp  at  normal 
incidence  for  100  nm  thick  films  is  plotted  in  Fig.  7. 
Reflectance  of  the  IAD  films  drops  by  ~2%  except  at 
the  1.5-eV  interband  absorption. 

Optical  properties  of  metal  thin  films  in  the  spectral 
range  of  intraband  transitions  can  be  described  phe¬ 
nomenologically  by  the  Drude  free-electron  model.17 
At  near  IR  frequencies  W  »  l/r),  "’  the  real  and  imagi¬ 
nary  parts  of  complex  dielectric  constants  can  be  writ¬ 
ten  as  follows: 

t,  a  (1) 


15  July  1989  /  Vol.  28.  No.  14  /  APPLIED  OPTICS 


2771 


Fig.  4.  (a)  Pseudorefractive  index  and  (b)  pseudoextinction  coeffi¬ 
cient  of  Ag  films  measured  by  the  SE  plotted  as  a  function  of  incident 
photon  energy. 


<2  a  uJ/(u3t),  (2) 

where  <o  represents  core  interband  transitions,  is 

the  plasma  angular  frequency  for  intraband  transi¬ 

tions,  r  is  the  intraband  relaxation  time,  and  w  is  the 
angular  frequency  of  the  light.  Ep  and  E  are  the 
energies  of  the  plasma  frequency  and  incident  wave, 
respectively.  The  intraband  plasma  frequency  of  the 
conduction  electrons  wp  is  given  by 

w  »  v4*/V»7m,  (3) 


Fig.  5.  Reflectance  of  Ag  films  calculated  from  the  pseudorefrac- 
tive  index  and  pseudoextinction  coefficient  plotted  as  a  function  of 
incident  photon  energy. 

where  N  is  the  effective  number  density  of  electrons 
participating  in  the  intraband  transitions,  and  m  and  e 
are  the  free-electron  mass  and  charge,  respectively. 
Also,  the  dc  conductivity  a  in  the  Drude  free-electron 
model  is  expressed  by 

«  “  u’r/4ir.  (4) 

Using  Eq.  ( 1 ),  we  can  determine  oip  of  a  metal  film  from 
the  slope  of  a  plot  of  vs  (l/£-).  From  up  and  Eq. 
(3),  A/  can  be  obtained. 

The  plots  of  -fi  vs  (1  /EP)  calculated  from  the  pseu¬ 
dodielectric  constants  for  Ag  and  A1  films  in  the  intra¬ 
band  transitions  are  shown  in  Figs.  8(a)  and  (b),  re¬ 
spectively.  Ep  of  the  IAD  Ag  film  at  P  =  0.48  \  e V  is  7.6 
eV,  while  that  of  the  non-IAD  film  is  7.7  eV.  For  the  A1 
films,  Ep  drops  from  10.7  eV  of  the  non-IAD  film  to 
10.1  eV  of  the  IAD  film  at  P  =  0-46yeV.  The  slope  for 
the  IAD  films  decreases  as  P  increases.  The  decrease 
of  Ep  means  a  reduction  in  N,  which  has  to  do  with 
changes  in  the  microstructure  of  the  metal  film.30-29 
The  fraction  of  voids  in  the  A1  films  is  calculated  by 
using  the  effective  medium  approximation  and  linear 
regressional  analysis,3'''35  and  the  results  are  shown  in 
Table  II.  Because  the  A1  films  had  been  exposed  to 
the  atmosphere  for  a  few  months  before  the  SE  mea¬ 
surement  was  performed,  a  thin  AI2O3  layer  on  top  of 
the  A1  film  was  taken  into  account.38  Since  the  100- 
nm  thickness  of  A1  film  is  much  larger  than  the  skin 
depth  of  A1  (~10  nm),  the  incident  wave  does  not  see 
the  glass  substrate.  Hence  bulk  A1  plus  voids  is  as¬ 
sumed  as  the  substrate.  In  Fig.  9,  the  measured  dielec¬ 
tric  constants  of  the  non-IAD  Al  film  are  in  good  agree¬ 
ment  with  those  calculated  from  the  data  modeled  in 
Table  II.  The  fraction  of  voids  increases  as  P  in¬ 
creases,  while  the  thickness  of  the  oxide  layer  does  not 
seem  to  correlate  with  P. 
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Fig.  6.  (a)  Pseudorefractive  index  and  (b)  pseudoextinction  coeffi¬ 
cient  of  A1  films  measured  by  the  SE  plotted  as  a  function  of  incident 
photon  energy. 

Analogous  to  the  reduction  of  interband  absorption 
in  the  IAD  Ag  films,  the  1.5-eV  interband  absorption  of 
the  IAD  A1  films  decreases.  Interestingly,  the  1.5-eV 
interband  absorption  in  the  A!  films  tends  to  disappear 
as  P  increases.  Similar  observations  for  A1  films  on 
cold  substrates  were  reported  by  Bernland.39  In  his 
work,  as  the  substrate  temperature  decreased  from  198 
to  25  K,  the  1.5-eV  interband  absorption  decreased 
and  finally  disappeared  at  25  K.  Possible  explana¬ 
tions  were  a  lack  of  [200]  translational  symmetry  at  25 
K  and  an  intermediate  state  between  a  crystalline  and 


Fig.  7.  Reflectance  of  Al  films  calculated  from  the  pseudorefractive 
index  and  pseudoextinction  coefficient  plotted  as  a  function  of  inci¬ 
dent  photon  energy. 


an  amorphous  state  in  the  Al  film  on  the  cold  sub¬ 
strate.3940  When  we  investigated  our  Al  films  with  x- 
ray  diffraction  (XRD),  we  did  not  observe  a  [200]  ori¬ 
entation  for  both  non-IAD  and  IAD  Al  films.  Thus 
these  explanations  do  not  seem  applicable  to  our  films. 
If  we  assume  the  reflectance  drop  for  the  intraband 
transitions  due  to  the  fraction  of  voids  is  extended  to 
the  interband  transitions,  the  reduction  of  the  1.5-eV 
interband  absorption  in  the  IAD  Al  films  should  also 
be  due  to  these  voids.  Generally,  the  presence  of  voids 
in  a  thin  film  means  that  the  characteristics  of  the  film 
deviate  from  the  bulk,29-34  Since  the  1.5-eV  interband 
absorption  is  a  characteristic  of  the  bulk,  the  presence 
of  voids  may  influence  it  so  that  the  1.5-eV  interband 
absorption  decreases  as  the  void  fraction  increases. 

Given  that  the  fraction  of  voids  correlates  well  with 
P,  we  believe  that  the  presence  of  voids  in  the  film 
induced  by  Ar  ions  is  an  important  factor  which  affects 
the  effective  number  density  of  electrons  in  the  intra¬ 
band  transitions  and  thus  the  optical  constants. 

IV.  X-Ray  Diffraction  Analysis 

A  Bragg-Brentano  x-ray  diffractometer  (XRD)  was 
used  to  investigate  the  microstructure  of  polycrystal¬ 
line  Ag  and  Al  films  such  as  the  strongest  orientation, 
lattice  spacing,  and  grain  size.  In  addition,  the  strain 
and  average  stress  in  the  film  plane  were  calculated 
from  the  change  in  lattice  spacing. 

For  Ag  films,  the  (Ill)  diffraction  profile  showed  the 
strongest  orientation.  The  intensity  of  (200)  orienta¬ 
tion  is  about  one-seventh  of  the  intensity  of  the  (111) 
orientation.  Thus  we  consider  only  the  (111)  orienta¬ 
tion. 

The  situation  for  A!  films  was  more  complicated. 
Thin  Al  films  (15  nm,  the  same  thickness  at  which  the 
optical  constants  of  Al  films  were  measured  using  the 
SPR  technique)  showed  no  diffraction  peaks.  Since 
these  thin  Al  films  had  been  exposed  to  the  atmo 
sphere  for  a  while  before  the  XRD  measurements,  an 
oxide  layer  comprising  an  average  of  20-30  at.  %  of  the 
film  (as  measured  by  the  Rutherford  backscattering 
spectrometry)  formed  at  the  air-film  interface.  It 
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Fig.  8.  -<i  plotted  as  a  function  of  (1/E2)  for  (a)  Ag  films  and  (b)  A1 
films. 


would  be  possible  that  the  thin  oxide  layer  disturbed 
the  structure  of  very  thin  A1  films  and  resulted  in  no 
diffraction  profiles.  For  100-nm  thick  A1  films,  the 
oxide  layer  comprised  only  an  average  of  2-3  at.  %  of 
the  film,  and  thus  a  diffraction  profile  was  observed 
but  only  the  (111)  orientation. 


Table  II.  Fraction  t  of  Voids  in  Al  Films  and  (he  Thickness  /  of  an 
Aluminum  Oxide  Layer  are  Simulated;  the  Thickness  of  Film  is  100  nm 


P(.JV) 

/(t!) 

t  (nm) 

0 

12.4 

3.3 

0.12 

18.2 

3.5 

0.46 

20.7 

2.4 

Fig  9.  Measured  pseudorefractive  index  and  pseudoextinction  co¬ 
efficient  compared  with  those  calculated  from  the  model  in  Table  11. 


Lattice  spacing  dm  perpendicular  to  the  film  sur¬ 
face  was  calculated  by  the  Bragg  diffraction  formula. 
The  variation  of  dm  for  Ag  films  is  shown  in  Table  III 
in  terms  of  P.  When  P  is  small,  dm  does  not  change 
significantly  from  that  of  the  non-IAD  Ag  film.  When 
P  exceeds  0.2  \eV,  however,  dm  increases  and  levels 
off  at  ?  360  A,  which  is  slightly  larger  than  the  value  for 
Ag  powder  given  by  JCPDS"  (2.339  A).  The  results 
for  Al  films  are  shown  in  Table  IV.  In  contrast  to  Ag 
films,  dni  of  IAD  Al  films  keeps  expanding  slightly  as 
P increases. 


Table  III.  Ag  Films;  the  Lattice  Spacing  Strain  («,)  r-erpendicular 

to  the  Film  Surface,  and  Strain  t,  in  the  Film  Plane  are  Listed;  (he 
Thickness  ol  Ag  Films  Is  SO  nm 


r  uiV) 

dm  (A) 

<;  (10”) 

I,  (10-’l 

0 

2.354 

-2.1200 

2.40 

0.18 

2.354 

-2.1200 

2.40 

0.22 

2.358 

-0.4240 

0.48 

0.29 

2.360 

0.4239 

-4.80 

0.30 

2.360 

0.4239 

-4.80 

0.37 

2.360 

0.4239 

-4.80 

0.48 

2.360 

0.4239 

-4.80 

Powder 

2.359 

Table  IV.  Al  Films;  the  Lattice  Spacing  d111(  Strain  «,  Perpendicular  to 
the  Film  Surface,  and  Strain  <(  In  the  Film  Plane  are  Listed;  the  Thickness 
of  Al  Films  is  100  nm 

PUeV) 

dm  (A) 

UO-’t 

<.  (10’ ') 

0 

2.343 

2.14 

2.01 

0.06 

2.344 

2.57 

2.42 

0.12 

2.344 

2.57 

2.42 

0.29 

2.345 

2.99 

2.82 

0,46 

Powder 

2.346 

2.338 

3.42 

3.22 
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The  change  in  lattice  spacing  is  related  to  the  strain 
t,  perpendicular  to  the  film  plane  as  follows4-: 

<,-(<*-  da)ld„ 

where  d  and  d0  are  the  strained  and  unstrained  lattice 
spacings,  respectively. 

If  we  assume  that  the  film  plane  is  elastically  isotro¬ 
pic  and  follow  the  developments  of  Vook  and  Wilt43 
and  Huang  eta/.,14  the  average  stress  <r.„  and  the  strain 
e,  in  the  film  plane  can  be  written  by 

<,  *  <r-v(S,,  + 

where  Sn,  S12,  and  S:!),  are  the  elastic  constants  in  the 
[111]  orientation  calculated  by  rotating  the  elastic  con¬ 
stants  for  the  [001]  orientation.1’1  The  resulting  t,  and 
«,  for  Ag  and  A1  films  are  listed  in  Tables  III  and  IV, 
respectively.  The  average  stress  in  the  film  plane  is 
plotted  in  Fig.  10  as  a  function  of  P.  Non-IAD  and 
lightly  bombarded  IAD  Ag  films  exhibit  tensile  stress, 
while  IAD  Ag  films  above  P  =  0.2  VeV  show  a  constant 
compressive  stress  value  independent  of  P.  Similar 
stress  reversal  has  been  reported  in  IAD  Nb  films,46 
dual  ion  beam  sputtered  Ag  films,44  and  other  sput¬ 
tered  metal  films.47,43  In  contrast  to  Ag  films,  the 
stress  of  IAD  A1  films  is  compressive  throughout  the 
entire  P  range  and  increases  slightly  with  P. 

It  appears  that  the  compressive  stress  of  IAD  A1 
films  correlates  well  with  P  while  the  stress  reversal 
from  tensile  to  compressive  for  IAD  Ag  films  requires  a 
threshold  value  of  P  Recently,  Windischmann--  pro¬ 
posed  a  simple  scaling  law  which  relates  the  compres¬ 
sive  stress  for  polvcrystalline  thin  films  prepared  by 
ion  beam  sputtering  with  the  momentum.  Also,  Hoff¬ 
man  and  Gaerltner11’ demonstrated  that  the  important 


parameter  for  reversal  of  stress  from  tensile  to  com¬ 
pressive  is  related  to  the  momentum  transfer  rather 
than  energy  by  investigating  evaporated  Cr  films  bom¬ 
barded  with  inert  gases.  Our  results  for  stress  analy¬ 
ses  of  IAD  Ag  and  A1  films  appear  to  support  Win- 
dischmann’s  scaling  law  and  be  in  good  agreement  with 
Hoffman  and  Gaerttner’  results. 

In  addition  to  strain  and  stress,  XRD  allows  us  to 
examine  the  grain  size  in  our  films.  Grain  size  Lui  is 
given  by  the  Scherrer  equation,4-  Lm  =  k\/0  cos 0, 
where  k  is  a  constant  (taken  to  be  11 ,  X  is  1 .542  k  (CuKo 
radiation),  0  is  the  Bragg  diffraction  angle,  and  3  is  the 
angular  width  of  the  pure  diffraction  profile  on  the  20 
scale.  We  corrected  for  instrumental  broadening  by 
assuming  Gaussian  profiles.4-  As  shown  in  Fig.  11, 
L hi  of  IAD  Ag  films  is  smaller  than  that  of  the  non- 
IAD  film,  and  it  does  not  change  significantly  within 
the  experimental  error  as  P  increases.  On  the  other 
hand,  Lm  of  IAD  A1  films  decreases  monotonically  as 
P  increases  and  at  P  =  0.46  NeV  Lm  of  the  IAD  A1  film 
was  only  ~67%  of  that  of  the  non-IAD  film. 

V.  Electrical  Resistivity  Measurement 

The  electrical  resistivity  p  of  Ag  and  A1  films  at 
ambient  temperature  was  measured  by  the  van  der 
Pauw  four-point  method.49  With  the  appropriate  cor¬ 
rection  factors,  p  of  a  flat  metal  film  of  arbitrary  shape 
can  be  measured  if  the  four  contacts  are  at  the  circum¬ 
ference  of  the  homogeneous  film.’10  Resistivity  for  Ag 
and  Al  films  is  shown  in  Fig.  12  as  a  function  of  P.  pof 
IAD  films  increases  as  P  increases. 

It  is  known  that  the  electrical  resistivity  of  a  thin 
metal  film  is  larger  than  that  of  the  bulk  and  decreases 
as  the  thickness  of  the  film  increases.  Fuchs  devel¬ 
oped  a  free-electron  model  to  explain  this  phenome¬ 
non  by  considering  the  size  effects  of  a  thin  film.  ’1  In 
his  theory,  the  thin  film  is  assumed  to  be  a  homoge¬ 
neous  slab,  and  electrons  are  scattered  not  only  in  the 
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Fig.  12.  Resistivity  of  Ag  and  A1  films  as  a  function  of  P. 

bulk  but  also  by  the  smooth  film  surfaces.  Real  poly- 
crystalline  films,  however,  have  different  structures 
from  the  corresponding  bulk.  Modified  models,  based 
on  Fuchs’s  theory  and  taking  into  account  surface 
roughness,52  grain  boundaries,53  and  charge  carriers 
localization,54  were  later  proposed.  These  also  are  not 
complete  because  the  polycrystalline  films  studied 
have  defects,  dislocations,  and  impurities.  With  this 
in  mind,  Sambles55  pointed  out  that  extreme  care  must 
be  taken  in  analyzing  and  interpreting  the  data  of 
resistivity  measurement. 

Since  Ar-ion  bombardment  on  the  growing  Ag  and 
A1  films  introduces  microstructural  changes,  the  pure 
Fuchs  theory  is  not  valid  for  our  films.  The  rms 
surface  roughness  of  our  films  was  ~1  nm  as  measured 
by  a  noncontact  surface  profiler  (WYKO  TOPO  2D) 
and  no  significant  difference  in  the  surface  roughness 
between  non-IAD  and  IAD  films  was  observed.  We 
can  thus  exclude  surface  roughness  as  a  factor  in  the 
differences  in  p  between  IAD  and  non-IAD  films. 
Also,  surface  scattering  due  to  the  film  thickness  ef¬ 
fects  is  relatively  negligible  because  the  grain  size  of  Ag 
and  A1  films  is  smaller  than  the  thicknesses  ( — 50  and 
100  nm  for  Ag  and  A1  films,  respectively),  and  so  vol¬ 
ume  scattering  dominates. 

When  Parmigiani  et  al.20  applied  Mayadas’s53  grain 
boundary  scattering  model  to  ion  bombarded  Ag  films 
prepared  by  dual  ion  beam  sputtering,  it  was  found 
that  the  reflection  coefficient  Rt  at  the  grain  boundary 
depended  on  the  deposition  conditions.  Since  the 
Mayadas’s  model  assumes  a  constant  R,  for  films  of  the 
same  metal,  the  different  Rf  for  each  film  was  attribut¬ 
ed  to  the  different  potential  barrier  at  the  grain  bound¬ 
aries  due  to  the  structural  changes. 

In  this  paper,  we  used  Vancea’s  model54:  that  the 
reduced  effective  number  density  of  electrons  N  par¬ 
ticipating  in  the  dc  conductivity  is  a  cause  for  the 
increase  in  resistivity  of  metal  thin  films.  Vancea's 


model  seems  to  be  an  important  factor  for  the  increase 
in  p  of  our  films  because  we  observed  the  reduction  of 
N  in  the  analysis  of  the  real  part  of  complex  dielectric 
■constants,  which  is  related  to  the  decrease  in  the  dc 
conductivity  of  the  free  electron  model  by  Eq.  (4). 
Also,  in  Vancea’s  model,  each  grain  acts  as  a  potential 
well  with  the  boundaries  as  barriers.  Free  electrons 
which  pass  through  the  potential  barrier  contribute  to 
the  dc  conductivity.  As  the  grain  size  decreases,  the 
number  of  barriers  increases,  more  electrons  are 
trapped  in  grains,  and  the  number  of  effective  elec¬ 
trons  participating  in  the  dc  conductivity  decreases. 
This  model  thus  relates  the  reduction  of  N  and  the 
decrease  in  grain  size  to  the  increase  in  p.  Since  the 
reduction  of  N  implies  an  increase  in  void  fraction  as 
discussed  in  the  SE  section,  the  presence  of  a  larger 
fraction  of  voids  in  IAD  films  seems  to  contribute  to 
the  increase  in  electric  resistivity  via  Vancea’s  model. 

VI.  Rutherford  Backscatterlng  Spectrometry  Analysis 

The  stoichiometry  of  our  films  was  measured  by 
Rutherford  backscattering  spectrometry  (RBS).  In¬ 
cident  singly  ionized  helium  (4He+)  ions  with  1.892- 
MeV  energy  are  backscattered  by  a  thin  film  sample. 
The  number  of  scattered  ions  per  unit  solid  angle  per 
incident  ion  is  measured  and  compared  with  the 
known  Rutherford  scattering  cross  section.  A  simple 
calculation  then  yields  the  absolute  number  per  unit 
area  of  each  atomic  species  in  the  film.56  The  film 
stoichiometry  thus  can  be  measured  very  accurately. 
Graphite  substrates  were  used  because  the  Si  and  O 
peaks  in  a  glass  substrate  would  overlap  the  Al  and  O 
peaks  in  the  films  and  thus  reduce  the  accuracy. 

The  at.  %of  Ar  in  Ag  and  Al  films  is  plotted  in  Fig.  13 
in  terms  of  P.  Even  though  the  amount  of  trapped  Ar 
is  <2  at.  %,  it  responds  well  to  P  in  both  Ag  and  Al 
films.  The  average  at.  %  of  oxygen  in  Ag  films  is  <0.1% 
and,  therefore,  negligible.  Al  films,  on  the  other  hand, 
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have  most  of  their  oxygen  at  the  surface  with  the 
variation  of  average  at.  %  of  oxygen  from  run  to  run 
being  <0.2%.  Therefore,  it  seems  that  Ar  ion  bom¬ 
bardment  on  A1  films  does  not  influence  the  natural 
oxidation  of  the  A1  surface  as  discussed  in  modeling  of 
A1  films  in  the  SE  section. 

It  is  known  that  the  trapped  Ar  atoms  in  films  under 
ion  bombardment  reside  at  grain  boundaries.20  Ac¬ 
cording  to  Vancea’s  model,  Ar  atoms  trapped  at  grain 
boundaries  increase  the  height  of  potential  barriers, 
localizing  more  electrons  in  grains.  Alternatively,  Ar 
atoms  at  grain  boundaries  act  as  electron  traps,  in¬ 
creasing  the  electrical  resistivity.  Also,  Ar  atoms  at 
grain  boundaries  prevent  diffusion  of  metal  atoms 
across  grain  boundaries  and  so  limit  grain  growth. 
However,  we  can  speculate  that  since  Ar  ions  in  films 
have  energy  greater  than  the  adatoms  in  the  film  and 
are  moving  with  a  preferred  direction  into  the  film, 
they  do  not  necessarily  all  reside  at  grain  boundaries 
and  some  may  reside  within  grains.  If  this  hypothesis 
is  correct,  the  small  amount  of  Ar  in  grains  cannot  be 
neglected  as  a  factor  in  the  change  of  electrical  resistiv¬ 
ity  and  stress.  Therefore,  in  either  case,  the  inclusion 
of  Ar  in  films  seems  to  affect  electrical  and  structural 
properties. 

VII.  Conclusions 

The  following  effects  are  observed  in  IAD  Ag  and  Al 
films  as  the  momentum  transferred  to  the  growing  film 
by  the  Ar-ion  beam  per  arriving  metal  atom  increases: 

(1)  In  the  intraband  transitions,  the  refractive  in¬ 
dex  increases  and  extinction  coefficient  decreases,  re¬ 
sulting  in  a  decrease  in  reflectance.  The  interband 
absorption,  however,  decreases  and  reflectance  in¬ 
creases.  The  changes  in  optical  constants  are  consis¬ 
tent  with  the  increase  in  the  void  fraction  of  the  films 
and  the  reduction  of  the  effective  number  density  of 
free  electrons  in  the  intraband  transitions. 

(2)  The  lattice  spacing  perpendicular  to  the  film 
surface  expands.  The  stress  changes  from  tensile  to 
compressive  for  Ag  films.  The  compressive  stress  of 
Al  films  increases  slightly.  The  grain  size  decreases 
for  both  types  of  film.  The  dependence  of  stress  on  P 
supports  Windischmann’s  scaling  law,  which  relates 
stress  to  the  momentum. 

(3)  The  electrical  resistivity  increases.  This  is  at¬ 
tributable  to  the  reduced  effective  number  density  of 
free  electrons.  The  decrease  in  grain  size  increases  the 
number  of  potential  barriers,  which  results  in  the  de¬ 
crease  in  the  effective  number  density  of  free  electrons. 
This  is  in  good  agreement  with  the  observed  reduction 
of  the  effective  number  density  of  free  electrons  in  the 
intraband  transitions  resulting  from  the  increase  in 
the  fraction  of  voids. 

(4)  The  amount  of  trapped  Ar  increases  slightly. 

From  the  results  of  this  study,  IAD  is  likely  to  make 

Ag  and  Al  metal  thin  films  less  bulklike  by  changing 
their  microstructure.  Also,  many  of  their  optical,  elec¬ 
trical,  and  microstructural  properties  correlate  well 
with  the  momentum  transferred  to  the  growing  film  by 
the  Ar  ions  per  arriving  metal  atom,  suggesting  that 


the  momentum  is  an  important  physical  parameter  in 
describing  the  influence  of  ion  beam  on  growing  thin 
films  and  determining  the  characteristics  of  IAD  thin 
metal  films. 
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Reactive  ion  assisted  deposition  of  aluminum  oxynitride 
thin  films 

Chang  Kwon  Hwangbo,  Linda  J.  Lingg,  John  P.  Lehan,  H.  Angus  Macleod,  and  F.  Suits 


Optical  properties,  stoichiometry,  chemical  bonding  states,  and  crystal  structure  of  aluminum  oxynitride 
( A10,N,.)  thin  films  prepared  by  reactive  ion  assisted  deposition  were  investigated.  The  results  show  that  by 
controlling  the  amount  of  reactive  gases  the  refractive  index  of  aluminum  oxynitride  films  at  550  nm  is  able  to 
be  varied  from  1.S5  to  1.83  with  a  very  small  extinction  coefficient.  Variations  of  optical  constants  and 
chemical  bonding  states  of  aluminum  oxynitride  films  are  related  to  the  stoichiometry.  From  an  x-ray 
photoelectron  spectroscopy  analysis  it  is  observed  that  our  aluminum  oxynitride  film  is  not  simply  a  mixture 
of  aluminum  oxide  and  aluminum  nitride  but  a  continuously  variable  compound.  The  aluminum  oxynitride 
films  are  amorphous  from  an  x-ray  diffraction  analysis.  A  rugate  filter  using  a  step  index  profile  of  aluminum 
oxynitride  films  was  fabricated  by  nitrogen  ion  beam  bombardment  of  a  growing  A1  film  with  backfill  oxygen 
pressure  as  the  sole  variation.  This  filter  shows  a  high  resistivity  to  atmospheric  moisture  adsorption, 
suggesting  that  the  packing  density  of  aluminum  oxynitride  films  is  close  to  unity  and  the  energetic  ion 
bombardment  densities  the  Him  as  well  as  forming  the  compound. 


I.  Introduction 

Aluminum  oxynitride  (A10,NV)  can  be  used  as  an 
optical  thin  film  with  variable  optical  constants  be¬ 
tween  those  of  aluminum  oxide  and  aluminum  nitride. 
Rugate  filters,  which  are  designed  with  a  sinusoidal 
refractive  index  profile,  and  other  types  of  inhomogen¬ 
eous  films  can  be  made  by  varying  the  refractive  index 
of  an  aluminum  oxynitride  film  continuously.  A10.N, 
can  also  be  used  as  a  tunable  medium  refractive  index 
material  in  the  near  UV  and  visible,  depending  on  the 
stoichiometry  of  A10xNj.  and  the  application.  Alumi¬ 
num  oxynitride  has  not  been  studied  in  detail  when 
compared  to  aluminum  oxide  and  aluminum  nitride, 
and  relatively  few  papers  have  been  published  on  it. 
Stoichiometry  and  optical  properties  of  aluminum  ni¬ 
tride  and  oxynitride  films  prepared  by  ion  assisted 
deposition  (IAD)  have  been  studied  in  this  laboratory 
in  the  past.1  Composite  aluminum  oxynitride 
[(Al)03)i.c-(AIN)r]  films  prepared  by  laser  assisted 
CVD  showed  a  refractive  index  variation  from  1.67  for 
AliO  i  to  2.08  for  AIN  at  300  nm  by  changing  the  sub¬ 
strate  temperature.2  Also  ion  beam  sputtered  aiumi- 
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num  oxynitride  films  have  been  used  to  encapsulate 
GaAs  substrates  at  high  temperture/' 

In  this  study  we  made  aluminum  oxynitride  films  by 
bombarding  a  growing  A1  film  with  a  nitrogen  ion  beam 
ora  mixture  of  nitrogen-oxygen  ion  beam  in  an  oxygen 
backfill.  We  investigated  the  optical  properties,  stoi¬ 
chiometry,  chemical  bonding  states,  and  crystal  struc¬ 
ture  of  the  films.  To  illustrate  an  application  of  alumi¬ 
num  oxynitride  films,  we  made  a  rugate  filter  using  a 
step  index  profile  and  compared  its  performance  to  the 
predicted  value  of  a  simulated  one. 

II.  Deposition  Conditions 

Aluminum  oxynitride  films  in  our  experiments  were 
deposited  by  reactive  IAD  in  a  Balzers  760  box  coater. 
which  has  been  described  elsewhere. The  base 
pressure  was  ~7  X  10*"  mbar.  The  chamber  tempera¬ 
ture  was  100°C,  and  its  variation  during  a  single  layer 
deposition  was  <10°C.  A1  was  e-beam  evaporated 
from  an  intermetallic  crucible.  The  condensation  rate 
was  3  A/s  as  measured  by  a  quartz  crystal  oscillator. 
Fused  silica,  graphite  slabs,  and  Si  wafers  were  used  as 
substrates  for  the  various  postdeposition  measure¬ 
ments.  Just  prior  to  deposition,  substrates  were  sput¬ 
ter  cleaned  in  vacuo  by  the  ion  beam  for  20  s. 

We  employed  a  3-cm  aperture  Kaufman  hot  cathode 
ion  source  for  two  cases  with  nitrogen  or  a  mixture  of 
nitrogen  and  oxygen  gases  as  the  bombarding  species. 
In  Case  1,  nitrogen  gas  was  admitted  through  the  ion 
source  until  a  pressure  of  8  X  10'  ’  mbar  was  obtained 
in  the  chamber.  Then  oxygen  gas  was  bled  through 
the  side  of  the  chamber  by  an  automatic  gas  pressure 
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Fig.  1.  Case  l:  variations  of  refractive  index  and  extinction  coeffi¬ 
cient  of  AlOrN,  thin  films  at  350  and  550nmas  the  backfilled  oxygen 
pressure  varies  from  0.5  to  1.4  X  10~4  mbar.  The  nitrogen  pressure 
through  the  ion  gun  was  kept  at  8  x  10"r,mbar.  Ion  energy  was  1250 
eV,  and  current  density  was  40  uA/cm2. 


Fig.  2.  Case  2:  variations  of  refractive  index  and  extinction  cocifi 
cient  of  A10,NV  thin  films  at  350  and  550  nm  as  the  pressure  oi  a 
mixture  of  nitrogen/oxvgen  gases  through  the  ion  gun  varies  from  |(^ 
8)  to  (8/0)  X  10"*’  mbar.  The  backfilled  oxygen  pressure  was  kept  at 
4  X  10“f'  mbar.  Ion  energy  was  1250  t  V,  and  current  density  was  -10 
*iA/cm2 


control  valve,  and  its  pressure  was  varied  from  0.4  to 
1.4  X  10"-1  mbar  to  control  the  stoichiometry  of  A10,NV 
films  and  thus  the  optical  constants.  In  Case  2,  a 
mixture  of  nitrogen  and  oxygen  gases  was  admitted 
through  the  ion  gun.  The  total  pressure  of  nitrogen 
and  oxygen  in  the  chamber  was  kept  at  8  X  10"5  mbar, 
while  the  ratio  of  two  gases  (nitrogen/oxygen)  was 
varied  from  8/0  to  0/8.  Also  backfill  oxygen  gas  at  a 
pressure  of  4  x  10~3  mbar  was  let  into  the  chamber.  In 
both  cases,  the  ion  beam  energy  was  1250  eV,  while  the 
current  density  as  measured  by  a  Faraday  cup  was  held 
at  40  nA/cm2. 

III.  Experiments  and  Results 

A.  Optical  Constants  Analysis 

Refractive  index  n  and  extinction  coefficient  k  of 
aluminum  oxynitride  films  were  obtained  by  measur¬ 
ing  the  transmittance  in  a  Cary  2415  spectrophotome¬ 
ter  (Varian)  and  using  an  envelope  method  which  is 
applicable  to  weakly  absorbing  thin  films.6 

In  addition  to  the  basic  envelope  method,  we  cor¬ 
rected  the  rear  surface  reflections  of  the  substrate  and 
assumed  a  Cauchy’s  dispersion  relation  of  the  form 

T  =■  A  +  B/x-  +  C/x* 

for  the  envelopes.  We  then  used  the  Cauchy’s  relation 
to  calculate  the  optical  constants  at  350  and  550  nm. 

In  Case  1  when  a  nitrogen  ion  beam  was  used  to 
bombard  a  growing  Al  film  in  an  oxygen  backfill,  varia¬ 
tions  of  n  and  k  at  350  and  550  nm  vs  backfilled  oxygen 
pressure  are  shown  in  Fig.  1.  As  the  backfilled  oxygen 
pressure  increases  from  0.5  to  1.4  X  lO"1  mbar,  n  at  550 
nm  decreases  from  1.81  to  1.68  and  k  at  550  nm  de¬ 


creases  from  8.2  to  2.4  X  10*'  As  the  amount  of 
reactive  oxygen  gas  increases,  n  and  k  of  AlO.N,  ap¬ 
proach  those  of  Al.'O  i  films. 

In  Case  2  when  the  mixture  of  nitrogen/oxvgen  ion 
beam  was  used  to  bombard  a  growing  Al  film  in  an 
oxygen  backfill,  variations  of  n  and  k  at  350  and  550  nm 
vs  different  nitrogen/oxvgen  mixture  ratios  through 
the  ion  gun  are  shown  in  Fig.  2.  As  the  pressure  of 
nitrogen/oxygen  gases  changes  from  1 0/8)  to  'S/0)  x 
10_r>  mbar,  n  at  550  nm  increases  from  1.65  to  1  S3  and 
k  at  550  nm  increases  1. 6-9.0  x  lO-1.  As  the  ratio  of 
nitrogen/oxygen  gases  increases,  n  and  k  of  A10,NV 
films  increase  toward  those  of  AIN. 

All  aluminum  oxynitride  films  showed  good  trans¬ 
mittance  in  the  near  UV  (>300  nm)  and  visible.  The 
cutoff  wavelength  was  — 190  nm. 

B.  Rutherford  Backscattering  Spectrometry  Analysis 

Stoichiometry  (x  and  y)  of  A10,NV  films  was  mea¬ 
sured  by  Rutherford  backscattering  spectrometrv 
(RBS).7 

In  Case  1,  variations  of  x  and  y  with  an  oxygen 
backfill  are  listed  in  Table  I.  As  the  backfilled  oxygen 
pressure  increases  from  0.5  to  1.4  x  10~J  mbar.  x  in¬ 
creases  from  0.85  to  1.28,  while  y  does  not  vary  appre¬ 
ciably.  Corresponding  changes  in  n  and  h  were  dis¬ 
cussed  in  the  previous  section  and  summarized  in  Fig. 
1. 

In  Case  2,  variations  of  x  and  y  are  listed  in  Table  II. 
As  the  ratio  (nitrogen/oxvgen)  is  varied  from  8/0  to  0/8 
while  keeping  the  sum  of  two  partial  pressures  at  8  X 
10'3  mbar,  x  decreases  from  0.84  to  1.55,  while  y  de¬ 
creases  from  0.60  to  0.  Corresponding  changes  in  n 
and  k  are  shown  in  Fig.  2. 
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C.  X-Ray  Photoelectron  Spectroscopy  Analysis 

While  RBS  provides  quite  accurate  information 
about  the  film  stoichiometry,  it  does  not  provide  infor¬ 
mation  about  the  chemical  bonding  states  of  the  ele¬ 
ments  in  the  film.  A  Perkin-Elmer  PHI  5100  ESCA 
system  for  x-ray  photoelectron  spectroscopy  (XPS) 
was  employed  to  investigate  the  chemical  bonding 
states  of  AlOjNv  films.  Mg  K«  x-rays  ( 1253.6  eV)  were 
used,  and  all  samples  were  sputter  etched  at  70°  for  10 
min  by  Ar  ions  at  4  keV  to  remove  surface  contami¬ 
nants  prior  to  making  measurements.  The  Al>p,  Oi„ 
and  Nls  peaks  were  investigated. 

We  were  primarily  interested  in  whether  aluminum 
oxynitride  is  simply  a  mixture  of  aluminum  oxide  and 
aluminum  nitride  such  as  (AI2O3)  i_c(A1N)c  or  an  alumi¬ 
num  oxynitride  compound  (A10XN;).  If  it  were  a  mix¬ 
ture  of  two  materials,  we  would  expect  the  Alip  peak  to 
be  a  doublet.  Thus  as  x  increases  and  y  decreases,  the 
intensity  of  one  peak  bonded  to  oxygen  increases,  and 
that  of  the  other  peak  bonded  to  nitrogen  decreases 
without  shifting  the  peak  position.  On  the  other  hand, 
if  it  were  a  compound,  we  would  see  only  a  single  Al2P 
peak. 

AlOjNy  films  investigated  by  XPS  were  prepared  by 
varying  the  ratio  of  nitrogen/oxygen  gas  mixture 
through  the  ion  gun  (Case  2  films  as  are  listed  in  Table 
II).  As  references,  we  made  aluminum  oxide  and  alu¬ 
minum  nitride  films  by  bombarding  an  oxygen  ion 
beam  on  a  growing  A1  film  in  an  oxygen  backfill  and  a 
nitrogen  ion  beam  on  a  growing  A1  film  in  a  nitrogen 


Table  I.  Variation  of  Stoichiometry  { x  and  y)  of  AIO.N,  Thin  Films  as  th* 
Pressure  of  Oxygen  Backfill  Changes 


Ion  source 
nitrogen 
(mbar) 

Backfill 

oxygen 

(mbar) 

(xi0-‘) 

x(±0.05) 

AJOfN, 

0  y(±0.05)“ 

8 

0.5 

0.85 

0.46 

8 

0.6 

0.97 

0.49 

8 

0.8 

1.07 

0. 16 

8 

1.0 

1.20 

0.11 

8 

1.2 

1.27 

0.43 

8 

1.4 

1.23 

0.45 

°  This  represents  the  maximum  excursion  from  film  to  film  given 
the  same  nominal  deposition  conditions. 

Table  II.  Variation  of  Slolchlomelry  (jr  and  y)  ol  AIO/fr  Thin  Films  as  the 
Ratio  of  a  Mixture  of  Nitrogen/Oxygen  Gases  Through  the  Ion  Gun  Changes 

Ion  source 

Nitrogen  Oxygen 

(mbar)  (mbar) 

(X10~s)  (Xio-'O 

Backfill 

oxygen 

(mbar) 

(xiO-'-l 

AIO.N, 

i(±0.051"  >(±0.05)" 

0 

8 

4 

1.52 

0 

2 

6 

4 

1.42 

0.18 

■1 

4 

4 

1.33 

0.27 

6 

2 

4 

1.21 

0.30 

7 

1 

■1 

1.00 

O.-IG 

8 

0 

4 

0.81 

0.60 

'*  This  represents  the  maximum  excursion  from  film  to  film  given 
the  same  nominal  deposition  conditions. 

Fig.  3.  Shift  of  the  Mip  peaks  of  A10,Ny  from  that  of  AlOi  52  to  that 
of  AIN1.23:  (a)  AlOjjtfi  (b)  AIOmzNo.jj;  (c)  AlO1.a3N0.27I  (d) 

ALOi  ooNo.46;  (*)  AlOo.wNo.co;  (0  AlNj.23. 

backfill,  respectively.  The  aluminum  nitride  film 
showed  superstoichiometry1  (AIN  1.2.1)  with  a  small 
amount  of  oxygen  an  the  air-film  interface  which 
would  be  sputter  etched  by  the  Ar  ions. 

The  Alip  peaks  of  A10x N,.  fil  ms  are  shown  in  Fig.  3  as 
x  and  y  vary.  Interesting,  the  A10xNy  films  have  a 
single  A1  ■>„  peak  which  shifts  from  the  high  binding 
energy  [Fig.  3(aa)]  of  A10i..y>  to  the  lower  binding  ener¬ 
gy  [Fig.  3 (/)]  of  AlNi.-v,  as  x  decreases  and  y  increases. 
When  x  is  larger  than  1.3  and  y  is  smaller  than  0.3,  we 
see  no  shift  in  the  Al;„  peak.  In  other  words,  it  seems 
that  a  small  amount  of  nitrogen  (y  <  0.3)  does  not 
appreciably  affect  the  chemical  bonding  states  of  Al, 
and  the  dominant  oxygen  (.t  >  1.3)  bonded  to  Al  be¬ 
haves  as  if  it  were  bonded  as  aluminum  oxide.  As  x 
decreases  below  1.3  and  y  increases  above  0.3,  the  Al;,, 
peak  of  AlOjNv  films  shifts  toward  that  of  AlNi.-a 
because  Al  is  shared  by  more  N  and  less  O.  Since  the 
binding  energy  of  a  single  Al;,  peak  of  A10XNV  film 
shifts  and  the  shape  of  the  Al_>p  peak  does  not  change 
significantly  as  x  and  y  vary,  A10tN\  is  a  compound. 

The  Ou  peaks  shown  in  Fig.  4  behave  quite  similarly 
to  the  AF,,  peaks  in  Fig.  3.  When  x  >  1.3  and  y  <  0.3. 
the  binding  energy  of  Ou  of  A10XNV  film  does  not  move 
(as  in  the  AI2,,  case).  This  is  further  evidence  that 
when  oxygen  is  above  a  threshold  value  (x  s  1.3  in  this 
case),  it  behaves  as  aluminum  oxide.  When  x  <  1.3 
and  y  >  0.3,  the  0_>s  peak  shifts  to  the  lower  binding 
energy. 

The  Nu  peaks  are  plotted  in  Fig.  5.  As  x  increases 
and  y  decreases  [from  (e)  to  (a)],  the  main  N\,  peak 
evolves  from  the  low  binding  energy  (400  e V)  of  A1N|  a 
to  the  high  binding  energy  at  which  both  N  and  O  are 
bonded  to  Al.  When  .r  -  1.42  and  y  =  0.13,  where  the 
Al  and  O  peaks  are  not  separated  from  those  of  AlOi 
a  new  peak  is  observed  al  the  G.2-eV  higher  binding 
energy  than  the  main  N(.  peak,  as  shown  in  Fig.  6. 
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Fig.  4.  Shift  of  the  Oj,  peaks  of  AIO,Ny:  (a)  AIOj.sj;  (b) 
AIOj.42N0.j8i  (c)  AlO1.33No.27!  (d)  A10j.ooNu.4g!  (fi)  AlOo.wNoxo- 


Fig.  5.  Shift  of  the  main  Ntl  peaks  of  A10iNv  and  appearance  of 
new  peaks  at  higher  binding  energy  for  (a)  and  (b):  (a)  AlOi  42N0.18.' 
(b)  AJO1.33N02?;  (c)  AIOj.ooNq  40!  (d)  AlOo84No.60!  (c)  AlNi  23* 

When  x  =  1.33  and  y  =  0.27,  where  the  A1  and  0  peaks 
are  the  same  as  those  of  AIO1.52,  the  new  peak  is  very 
small  and  the  main  peak  moves  0.1  eV  to  a  lower 
binding  energy  (Fig.  6).  When  x  <  1.3  and  y  >  0.3,  the 
new  peak  disappears  and  the  main  Njs  peak  shifts. 

D.  Infrared  Transmittance  and  X-Ray  Oiffractometry 
Analysis 

Aluminum  oxynitride  films  are  transparent  up  to  ~9 
nm  in  the  IR  spectral  range  without  much  absorption, 
as  measured  by  an  Analect  Instruments  FX-6200  Fou¬ 
rier  transform  infrared  spectrometer.  Our  aluminum 
oxynitride  films  exhibited  no  absorption  lines  of  alu- 
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Fig.  6.  Extended  view  of  (a),  (b),  and  (e)  in  Fig.  5. 


minum  oxide  or  aluminum  nitride  which  were  ob¬ 
served  by  Demiryont  e£  a!.2  in  composite  aluminum 
oxynitride  films.  This  agrees  with  the  results  from  our 
XPS  analysis  that  the  aluminum  oxynitride  film  pre¬ 
pared  by  reactive  IAD  is  not  a  mixture  of  aluminum 
oxide  and  aluminum  nitride  but  a  compound.  Also  we 
have  not  seen  any  special  absorption  peaks  which  can 
be  considered  as  aluminum  oxynitride  absorption 
peaks. 

From  Bragg  0-20  x-ray  diffractometer  analysis  it  was 
found  that  our  aluminum  oxide,  aluminum  oxynitride, 
and  aluminum  nitride  films  are  amorphous. 

IV.  Applications 

We  made  a  rugate  filter  by  using  the  step  index 
profile  of  Fig.  7.  Of  the  two  ways  we  investigated  to 
control  the  optical  constants,  we  chose  Case  1  where  we 
varied  the  oxygen  backfill  pressure  while  a  nitrogen  ion 
beam  bombarded  the  growing  A1  film.  We  chose  this 
for  the  practical  reason  that  the  lifetime  of  the  tung¬ 
sten  filaments  in  the  ion  gun  decreases  drastically 
when  oxygen  gas  is  used  as  a  bombarding  species.  The 
backfilled  oxygen  gas  pressure  was  varied  from  0.5  to 
1.4  X  10"4  mbar.  We  avoided  0.4  X  10-4  mbar  of 
oxygen  backfill  at  which  the  highest  refractive  index 
was  obtained  because  downward  fluctuations  in  oxy¬ 
gen  pressure  resulted  in  a  large  undesirable  absorption 
in  the  near  UV.  Thus  n  was  varied  from  1.S1  to  1.63. 
Each  cycle  has  ten  homogeneous  step  index  layers  and 
a  168.6-nm  physical  thickness.  Fifteen  cycles  were 
deposited  on  a  fused  silica  substrate.  Measured  re¬ 
flectance  R  and  transmittance  T and  theoretically  sim¬ 
ulated  R  and  T  from  the  assumed  index  profile  (Fig.  7) 
(neglecting  the  dispersion  of  aluminum  oxynitride 
fiims)  are  shown  in  Fig.  8.  Absorption  (A  =  1  -  R  -  T) 
of  the  filter  (2.5-^m  tolul  thickness)  at  the  peak  wave¬ 
length  is  found  to  be  <0.5%. 

The  discrepancy  in  T  and  R  between  real  and  simu¬ 
lated  filters  can  be  attributed  to  a  fluctuatiqn  of  back¬ 
filled  oxygen  pressure  during  each  layer's  deposition 
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Fig.  7.  Step-index  profile  of  A10,Ny  films  obtained  from  Fig.  1  for  a 
rugate  filter.  One  cycle  consists  of  ten  homogeneous  step  index 
layers  and  has  a  168-6-nm  physical  thickness.  A  sinusoidal  curve  is 
fitted  to  the  step  index  profile. 

and  the  nondispersive  calculation  used  for  the  simula¬ 
tion.  Also,  the  changes  in  refractive  index  were  not 
sharp  steps  but  had  some  slope,  because  the  pressure  of 
oxygen  gas  admitted  into  the  chamber  could  not  be 
increased  or  decreased  instantaneously.  Another  pos¬ 
sible  reason  for  the  discrepancy  is  that  the  nitrogen  ion 
current  density  at  the  substrate  on  each  layer  varied 
slightly  due  to  ion  scattering  from  the  different  oxygen 
backfill  pressure.  It  appears  that  these  slight  varia¬ 
tions  in  each  layer  reduce  the  peak-to-peak  variation  of 
refractive  index  profile  and  cause  a  decrease  in  the 
peak  reflectance  in  the  real  filter  when  compared  with 
the  simulated  one.  Better  control  of  optical  constants 
can  be  accomplished  by  computer  automated  control 
of  gas  pressure,  accounting  for  fluctuations  in  the  opti¬ 
cal  constants  and  thickness  correction. 

Moisture  adsorption  of  this  filter  from  the  atmo¬ 
sphere  was  investigated  by  measuring  the  peak  wave¬ 
lengths  of  the  filter  under  vacuum  (10  mTorr)  and  at 
atmospheric  pressure.  No  shift  of  the  peak  wave¬ 
length  was  observed.  This  indicates  that  the  packing 
density  of  aluminum  oxynitride  films  prepared  in  this 
fashion  is  close  to  unity. 

V.  Discussion 

For  the  first  time,  we  believe,  the  binding  energy 
shift  of  the  Al2p  peaks  of  A10,NV  films  from  aluminum 
oxide  to  aluminum  nitride  is  observed  as  the  stoichio¬ 
metry  changes.  The  differences  in  the  binding  ener¬ 
gies  of  the  M2n  peaks  of  AIOi..-,2  and  A1N|.2s  films  is  2  9 
e V.  In  contrast,  a  2-eV  difference  for  the  A!  >p  peaks  of 
Al  jO:i  and  AIN  films  was  observed  in  an  ionized  cluster 
beam  deposited  film.8 

Similarly,  new  peaks  at  6.5-eV  higher  binding  energy 
than  the  main  Ni,  peak  were  observed  when  AIN  pow¬ 
der  was  heated  to  between  700  and  900°C.  These  new 
peaks  were  attributed  to  oxynitrides.0  In  our  films  the 
new  peaks  (Fig.  6)  show  up  only  when  the  amount  of 


Fig.  8.  Reflectance/?  and  transmittance  T  of  the  real  and  simulated 
rugate  filters.  Fifteen  cycles  were  employed. 


nitrogen  is  small  and  the  amount  of  oxygen  is  close  to 
that  of  AlOt.52-  Since  A1  and  0  peaks  for  x  >  1.3  and  y 
<  0.3  are  at  the  same  binding  energies  and  have  the 
same  shape  as  those  of  A10L52  (Figs.  3  and  4),  possible 
explanations  for  the  new  peak  are  that  (1 )  this  peak  has 
to  do  with  a  reaction  between  nitrogen  and  oxygen  (not 
necessarily  NO  but  a  general  nitrogen-oxygen  intera- 
tion)  and  (2)  A1  needs  the  minimum  amount  of  oxygen 
(threshold  value  x  s  1.3  in  this  case)  to  act  like  AIO1.52. 
When*  =  1.42  and  y  =  0.18,  the  excess  oxygen  (x  >  1.3) 
reacts  with  nitrogen,  and  this  results  in  the  new  peak  at 
the  high  binding  energy  (408  eV),  while  the  rest  of  the 
nitrogen  reacts  with  A1  and  this  contributes  to  the 
main  Nu  peak  (at  401.8  eV)  (Fig.  6).  When  x  =  1.33 
and  y  =  0.27,  the  excess  oxygen  available  for  the  nitro¬ 
gen-oxygen  reaction  is  small.  Thus  the  new  peak 
shifts  to  low  binding  energy  (407  eV)  and  the  intensity 
is  low.  Since  more  nitrogen  (y  =  0.27)  reacts  with  A1 
than  in  the  above  case  (x  =  1.42  and  y  =  0.18),  the  main 
Nij  peak  moves  0.1  eV  to  lower  binding  energy.  As  x 
decreases  below  1.3  and  y  increases  above  0.3,  the  new 
peak  disappears  because  there  is  no  excess  oxygen  to 
react  with  the  nitrogen.  Also  since  both  oxygen  and 
nitrogen  readonly  with  A1  (fort  <  1.3  and  y  >  0.3),  the 
main  Ni,  peak  moves  in  proportion  to  y. 

Those  films  (x  >  1.3  and  y  <  0.3)  in  which  Abp,  Ots, 
and  the  main  N(J  peaks  are  at  the  same  binding  ener¬ 
gies  show  different  transmittances  and  thus  different 
optical  constants.  In  XPS  the  only  difference  be¬ 
tween  these  films  is  the  different  binding  energy  and 
heights  between  new  peaks.  Thus  this  new  peak  must 
be  significant  in  determining  the  different  optical  con¬ 
stants  for  these  films. 

It  is  known  that  when  an  energetic  N?  ion  bombards 
an  A1  surface,  it  is  neutralized  just  before  it  strikes  the 
surface  and  dissociated  into  two  nitrogen  atoms  as  it 
impacts.1011  These  two  nitrogen  atoms  are  implant¬ 
ed,  react  with  two  A1  atoms,  and  result  in  2A1N.  The 
kinetic  energy  of  each  nitrogen  atom  is  less  than  or 
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close  to  half  ot'  kinetic  energy  of  Ni  ion  because  the  N2 
ion  loses  its  kinetic  energy  when  it  collides  on  the 
surface  and  it  dissociates  into  two  nitrogen  atoms.12-14 
On  the  other  hand,  a  N  ion  in  a  nitrogen  ion  beam  is 
neutralized  and  reacts  with  an  A1  atom  with  the  same 
energy  as  it  is  given  by  the  ion  gun. 

In  our  study  we  made  an  A1  film  in  a  nitrogen  backfill 
(2  X  10"4  mbar)  to  check  whether  nitrogen  molecules 
react  with  A1  atoms.  For  this  particular  film  the  sub¬ 
strate  temperature  was  200°C  to  facilitate  the  reac¬ 
tion.  It  was  found  from  R8S  that  the  A1  film  consists 
of  A1  and  0  the  atomic  concentrations  of  which  are  98 
±  2%  and  2  ±  2%,  respectively.  Even  at  the  air-film 
and  film-substrate  interfaces  only  a  small  amount  of 
oxygen  was  found.  In  other  words,  the  nitrogen  in  the 
film  and  at  the  interfaces  is  under  the  detectable  limit 
of  RBS  (- — 0.4  at.  %).  Thus  A!  does  not  chemisorb 
molecular  nitrogen.  This  result  agrees  with  Winters 
and  Kay’s  classification14  of  A1  as  one  of  the  Class  2 
materials  which  do  not  chemisorb  molecular  nitrogen 
but  form  nitrides. 

Based  on  observations  of  Taylor  e£  a(.l0-1!  and  Win¬ 
ters  et  a/.,l2-ir’  we  can  speculate  what  the  nitrogen  ion 
beam  does  in  the  formation  of  aluminum  oxynitride 
films.  Many  molecular  nitrogen  ions  (1250  eV)  are 
neutralized  and  dissociated  into  nitrogen  atoms  on 
collision  with  the  film  surface.  The  resultant  atomic 
nitrogen  atoms  (<625  eV)  penetrate  the  growing  film, 
transfer  their  momenta  to  aluminum  and  oxygen  at¬ 
oms  by  collision  cascades,  and  densify  the  film,10  slow 
down,  and  finally  stop  and  react  with  aluminum  and 
oxygen  atoms.  On  the  other  hand,  the  atomic  nitrogen 
ions  (1250  eV)  from  the  ion  gun  are  neutralized,  and 
the  nitrogen  atoms  with  1250  eV  follow  a  similar  series 
of  steps. 

VI.  Conclusions 

Optical  constant  measurements,  RBS,  and  XPS 
were  combined  to  achieve  an  understanding  of  alumi¬ 
num  oxynitride  films  prepared  by  reactive  ion  assisted 
deposition.  The  refractive  index  at  550  nm  with  a  very 
small  extinction  coefficient  could  be  varied  from  1.65 
to  1.83,  and  a  relationship  between  refractive  index 
and  stoichiometry  was  obtained.  Aluminum  oxyni¬ 
tride  prepared  by  reactive  ion  assisted  deposition  is 
not  simply  a  mixture  of  aluminum  oxide  and  alumi¬ 
num  nitride  but  a  compound.  With  a  step  index  de¬ 
sign,  an  aluminum  oxynitride  rugate  filter  has  been 
successfully  fabricated.  Absorption  of  the  filter  was 
found  to  be  <0.5%.  No  vacuum-to-air  spectral  shift  of 
the  rugate  filter  was  observed,  suggesting  that  the 
packing  density  of  aluminum  oxynitride  films  is  close 
to  unity,  and  the  energetic  nitrogen  ion  bombardment 
densities  the  film  as  well  as  forming  the  compound. 
An  amorphous  aluminum  oxynitride  is  a  good  candi¬ 
date  material  for  an  inhomogeneous  film,  a  rugate 
filter,  and  a  medium  index  film,  because  the  absorp¬ 
tion  is  small  and  the  optical  constants  can  be  readily 
varied  from  those  of  aluminum  oxide  to  those  of  alumi¬ 
num  nitride  by  simply  varying  the  amount  of  reactive 
gases  introduced  in  the  ion  gun  or  alternatively  the 


backfilled  gas  while  the  ion  beam  bombards  a  growing 
Al  film. 
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Abstract 

It  is  necessary  to  take  into  account  scattering  phenomena  in  optical  systems  since  one  is 
interested  in  a  very  detailed  energy  balance.  We  recall  the  principles  of  the  vector  theories 
that  have  been  developed  to  predict  the  spatial  and  spectral  distribution  of  light  scattered  from 
an  optical  surface  or  a  multilayer  coating.  With  the  help  of  a  scattering  apparatus,  we 
emphasize  the  key  scattering  parameters  (roughness,  autocorrelation  length,  isotropy  degree, 

crosscorrelation  laws)  that  are  necessary  to  characterize  surface  defects  and  materials 
microstructure  in  thin  film  form.  In  spite  of  numerous  parameters  involved  in  the 

calculation,  particular  techniques  enable  us  to  extract  all  these  parameters  and  lead  to  good 
agreement  between  theory  and  experiment.  When  the  coating  is  made  of  a  high  number  of 

layers,  investigation  of  experimental  results  is  quite  more  difficult.  However,  with  the  help  of 

a  correctly  chosen  model,  we  can  strongly  reduce  the  number  of  scattering  parameters,  and 

this  allows  an  easy  comparison  between  theory  and  experiment. 

I  -  INTRODUCTION 

Drawing  up  a  very  detailed  energy  balance  (typically  at  better  than  10'3)  in  optical 

multilayers  remains  until  now  a  difficult  problem  because  of  numerous  parameters  that  are 

more  or  less  connected  to  the  microstructure  [1,  2]  of  materials  in  thin  film  form,  such  as 
absorption,  surface  and  volume  scattering,  adsorption,  index  inhomogeneity  and  anisotropy,  ... 

Even  in  the  simplest  case  of  a  single  layer  (figure  1)  deposited  on  a  transparent  substrate, 

we  must  write: 

1  =  R  +T  +  A  De  +  Dj 

where: 

-  R  and  T  are  the  reflection  and  transmission  coefficients  of  the  ensemble  (layer- 

substrate); 

-  A  characterizes  absorption  losses  [3,  4]; 

-  De  characterizes  external  scattering  losses  [S,  6],  which  represent  the  quantity  of 

scattered  light  that  can  merge  out  of  the  stack  and  that  is  easily  measurable; 

-  Di  characterizes  internal  scattering  losses  (fig,  I),  which  represent  the  quantity  of 

scattered  light  embedded  inside  the  stack,  and  that  can  be  absorbed  or  not. 

Fig.  1  -  Scattering  from  the  two  interfaces  (0)  and  (1)  of  a  single  layer 
with  thickness  e.  External  scattering  Dt  characterizes  scattered  waves 
that  can  merge  out  of  the  layer  and  that  can  be  easily  measurable. 

Internal  scattering  Di  characterizes  the  quantity  of  scattered  light  which 
is  imbedded  inside  the  layer  (scattering  at  large  angles  9  in  the  layer). 

In  the  same  way,  to  explain  the  differences  observed  between  calculated  and  measured 
optical  properties  of  these  layers,  we  are  often  constrained  to  take  account  of  a  gradient  index 
inside  the  layer  (7],  which  is  characteristic  of  the  material  growth  in  a  columnar  form.  Let  us 
notice  that  this  columnar  growth  can  also  induce  a  slight  anisotropy  of  the  refractive  index, 
that  can  be  measurable  with  guided  wave  techniques  [8]. 
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Our  laboratory  has  developed  experimental  and  theoretical  tools  to  investigate  each  of 
these  phenomena  [7,  8,  9,  10,  11,  12,  13],  which  permitted  us  to  demonstrate  the 
preponderence  of  surfacic  scattering  phenomena  in  the  limitation  of  optical  coating 
performances  [14]  (in  particular  for  high  quality  mirrors  and  multiple  cavity  Fabry-Perot 
filters  used  in  multi-demultiplexing  systems).  At  this  time,  we  have  at  our  disposal  a  scattering 
apparatus  [IS]  capable  of  measuring  the  scattered  light  in  whole  space  and  a  vector  theory 
[16]  to  investigate  experimental  results.  We  have  noted  an  excellent  agreement  between 
theory  and  experiment  in  many  cases,  especially  for  bare  substrates  [15],  antireflection 
layers  [17,  18]  and  multilayer  mirrors  [19].  However,  comparison  of  theory  and  experiment 
[20]  remains  quasi  impossible  in  the  case  of  complex  coating  designs,  because  of  the  high 
number  of  parameters  involved  in  the  calculation  (5p+4  parameters  for  a  p-layers  stack). 

On  this  subject,  we  will  show  here  how  a  correctly  chosen  model  enables  us  to  strongly 
reduce  the  number  of  scattering  parameters.  This  allows  computation  of  systematic  research 
of  the  scattering  parameters  that  lead  to  the  best  agreement  between  theory  and  experiment. 

In  a  more  general  way,  we  will  try  here  to  give  some  elements  of  response  to  the 
following  questions:  what  are  the  key  parameters  that  characterize  scattering  from  surfaces 
and  multilayer  coatings?  What  techniques  can  be  used  to  point  out  these  parameters?  How  can 
we  improve  the  theoretical  model  to  allow  an  easy  comparison  between  theory  and 
experiment,  and  consequently  extract  the  maximum  information  from  the  measurements?  So 
we  first  propose  an  up  to  date  review  which  will  permit  us,  without  entering  the  detail  of 
calculation,  to  present  some  recent  results  obtained  in  our  laboratory. 

II  -  VECTOR  SCATTERING  THEORY  FOR  OPTICAL  SURFACES  AND  MULTILAYER 
OPTICS:  RECALLS  AND  COMPLEMENTS 

Among  the  vector  theories  that  can  predict  the  spatial  distribution  of  light  scattered 
from  an  optical  surface,  Elson's  [5,  21]  and  Bousquet's  [16]  theories  take  a  particular  place 
because  they  present  the  advantage  that  they  can  also  be  applied  to  the  case  of  multilayer 
coatings.  -  It  is  certainly  useful  to  begin  by  some  recalls  concerning  the  fields  of  application  of 
these  two  theories,  that  are: 

-  the  height  fluctuations  h(r* )  of  the  surface  (figure  2)  are  assumed  to  be  very  small 
with  respect  to  the  wavelength  of  incident  light. 

-  the  slope  of  the  surface  defects  is  very  small  with  respect  to  unity. 

Fig.  2  -  Profile  h(Y )  of  a  rough  surface 
separating  2  medic  of  permittivities  e; 
and  ez  ■  (xoy)  is  the  mean  plane  of  the 
rough  surface. 

Both  mediums  (ex'  and  (e2)  (fig.  2)  are 
assumed  to  be  linear,  homogeneous  and 
isotropic.  We  want  also  to  emphasize  the 
fact  that  these  are  perturbative  theories, 
for  which  the  surface  is  considered  as  a 
deformation  of  a  -rerfectly  fiat  surface 
(with  no  roughness). 

Hence  it  is  obvious  that  these  theories  can  only  be  applied  to  the  case  of  very  good  quality 
optical  surfaces  (mic  -opolishes),  with  no  pits,  no  scratches,  no  dust  ... 

1)  Equivalence  between  Elson's  and  Bousquet's  theories 

Without  entering  the  detail  of  calculation  that  can  be  found  elsewhere  [22],  we  briefly 
indicate  here  the  reason  why  the  two  theories  [16,  21]  are  equivalent  for  the  calculation  of 
external  scattering.  In  fact,  resolution  of  Maxwell's  equations  show  that  the  field  scattered  by 
a  rough  surface  separating  two  media  is  equal  to  the  field  radiated  by  surfacic  currents 
distributed  on  a  perfectly  plane  (z  a  0)  surface  ("ideal  surface")  separating  the  same  media. 

The  scattered  field  li  d  being  defined  as  the  difference  between  the  total  field  £:  existing  in 

space  when  the  rough  surface  is  illuminated  and  the  "ideal"  field  1:°  which  would  exist  in  the 
same  conditions  if  the  surface  was  perfectly  flat,  it  follows  the  two  equations: 


rot  lld  a  j  <B  |i  i?d  +  \t  5(z) 
rot  i?da  -j  (0  81  +■  ~f  5(z) 
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(1-4) 

(1-b) 

where: 

- 1: d  and  it d  are  electric  and  magnetic  scattered  fields 

-  Kt  and  are  magnetic  and  electric  surfacic  currents 
•  S  is  the  Dirac  distribution 

-  ej  m  ei  +  («2  -  el)  H(z)  with  H  being  the  Heaviside  function 

-  e  and  ft  are  permittivity  and  permeability  of  the  media 
■  e  a  2a /X  with  X  being  the  wavelength  of  incident  light. 

According  to  the  calculation  method  that  is  employed  to  obtain  equations  (1),  the  source 

couple  (Si,  J* )  can  be  written  following  Elson: 

=  and  T*=J*t  (where  the  index  t  indicates  that  the  currents  are 

tmgmtial  to  the  ideal  surface  z  =  0) 
or  following  Bousquet: 

5t  =  I?  and  ?  =  +  Ty  (where  the  index  N  indicates  a  normal  component 

of  the  electric  current).  One  can  then  verify  that  the  tangential  components  of  the  fields 
scattered  by  the  source  couples  (hit .  J*  t)  and  ((?,  +  J*y)  are  identical  because  of  the  relation 

j  oa  82  Sit  ■  ^  A  (giud  JNl)t  * 

Consequently  these  two  theories  lead  to  identical  results  as  far  as  calculation  of  external 
scattering  is  concerned. 

2)  Results  for  a  single  surface 

Resolution  of  equation  (1)  is  made  in  the  Fourier  plane  and  leads  to  the  following 
expression  for  the  intensity  scattered  in  a  particular  direction  (9,  p)  of  space  (fig.  3): 

1(9,  $)  *  a(9,  <b)  y(9.  $) 

Fig.  3  Scattering  angles  9  and  p 

that  characterize  a  scattering 
direction,  (xoy)  is  the  mean  plane  of 
the  rough  surface.  (!}  J  is  the  spatial 
frequency  of  the  grating  responsible 
for  scattering  at  (9,  p). 


The  coefficient  a(9,  p)  is  called  "ideal  coefficient"  because  it  does  not  depend  on  the  surface 
defects,  but  only  on  the  illumination  and  observation  conditions.  The  roughness  spectrum  r(9, 
p)  is  the  Fourier  transform  of  the  autocorrelation  function  of  the  profile  h(r ),  and  contains 
all  accessible  information  on  the  interface  defects. 

To  permit  a  better  understanding,  it  is  convenient  to  consider  the  rough  surface  h(7*)  as 
the  sum  of  an  infinity  of  sinusoidal  gratings  (Fourier  decomposition).  In  the  case  of  normal 
incidence,  one  can  easily  show  that  the  fields  scattered  at  symmetrical  directions  (d.p)  and  (0, 
$  +  a)  take  their  origin  in  the  presence,  on  the  rough  surface,  of  one  sinusoidal  grating 
s<j(r*)  =»  2 1  h(#)l  cos  [■?.  -  ¥(?)] 

with: 

fi(^)»FT(h(7*)J=  |i(^)|  ei*^) 


j  COS  d 

1  sin  p 


Through  this  remark,  we  can  see  that  measurement  of  scattering  at  a  given  direction  9  will 
give  us  information  on  a  particular  defect  (fig.  3)  which  spatial  period  is  C  *  jiQ  0  .  This  will 
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be  essential  for  the  calculation  of  roughness  which  will  depend  on  the  different  kinds  of 
measurable  defects  (see  III-2). 

When  the  surface  is  illuminated  at  normal  incidence  with  natural  light,  the  ideal 
coefficient  only  depends  on  8  and  the  scattered  intensity  can  be  written  as:  1(8,  $)  a  a(8)  y(8,  $) 
or  also: 

I  ("?)  a  a(c)  yC?)  where  0=1"?  I .  If  we  rotate  now  the  surface  of  an  angle  o  around  its  normal, 
its  profile  must  then  be  written  in  the  (o  x  y  z)  coordinates  system  as  h'O?)  =  h(  R^Cr*)],  where 
Ra  is  a  rotation  of  angle  a  in  the  plane  oxy.  As  rotation  is  conserved  by  Fourier  Transform,  we 
can  then  measure  the  quantity 
Im(«b  3  Mo)  irtRci  O?)] 

that  is  Im(8,  6)  »  a(8)  7(8,  4  +  a)  *  1(8,  +  +  a) 

This  shows  that  it  is  equivalent  in  these  conditions  to  measure  scattering  in  the  whole  space 
or  to  limit  the  measurements  in  one  particular  plane  $  for  different  angular  positions  a  of  the 
sample  in  its  own  plane.  Our  scattering  apparatus  [15,  16]  is  based  on  this  principle  and  this 
is  also  quite  useful  for  investigation  of  roughness  anisotropy  (see  111-4).  In  the  case  of  an 
isotropic  surface,  the  scattering  distribution  will  be  obviously  symmetrical  with  respect  to  the 
normal  of  the  sample. 


3)  Case  of  a  multilayer  coating 

Ir  the  case  of  a  multilayer  coating  (fig.  4),  each  rough  interface  is  a  source  of 

scattered  light.  The  scattered  intensity  that  can  be  measured  is  due  to  the  superposition  of  the 
waves  emitted  by  these  surfaces  and  we  have  to  take  into  account  both  the  multiple 

reflections  of  each  scattered  wave  inside  the  multilayer  and  the  interference  phenomena 

between  all  these  scattering  sources.  The  preceding  theories  [16,  21]  can  be  easily 

generalized  to  the  case  of  multilayer  coatings  and  lead  to  the  following  expression  for  the 

scattered  intensity: 


1(8,  3  zLr  I  Ci  I  Yii  +  Ci  Cj  Yi  j 

i=0  i*j 

(Cj  »  conjugate  complex  number  of  Cj) 
where  p  is  the  number  of  layers  of  the  stack 
and  Yii  the  roughness  spectrum  of  interface 

(i)  (it  contains  all  information  relative  to  the 
defects  of  this  interface).  The  C,  coefficients 
are  called  "ideal  coefficients"  because  they 
only  depend  on  the  formula  of  the  coating  and 
on  the  illumination  and  observation 
conditions  (incidence.  wavelength  and 
polarization).  They  take  into  account  the 
multiple  reflections  of  the  scattered  waves 
inside  the  stack.  The  Yij  term  is  the  Fourier 
Transform  of  the  crosscorrelation  function 
between  the  profiles  of  surfaces  (i)  and  (j), 
and  takes  into  account  the  interference 
phenomena  between  the  scattering  sources  at 
interfaces  (i)  and  (j).  It  is  commonly  used  to 
define  a  crosscorrelation  coefficient  a;j  by  the 
relation  Yij  “  ctij  Yjj-  This  coefficient  can  be 
seen  as  a  degree  of  coherence  between  the 
scattering  sources  (i)  and  (j).  or  as  a  degree  of 
similitude  between  the  surface  profiles  (i)  and 

(j) . 


Fig.  4  •  Scattering  from  a  multilayer  stack 
with  p  layers  of  thicknesses  «i  and 
indices  ni.  Each  interface  is  a  source  of 
scattered  light. 


In  order  to  simplify,  it  is  generally  assumed  that  a;j  does  not  depend  on  spatial 
frequency,  so  that  we  are  interested  only  in  the  two  extreme  cases  of  uncorrelated  surfaces 
(a,j  *  0)  and  perfectly  correlated  surfaces  (a;j  •  1).  One  also  usually  consider  that  such 
coefficient  ajj  is  a  free  parameter,  independent  of  roughness.  We  will  see  in  section  V  how  a 
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more  realistic  model  can  be  developed.  In  a  general  way,  these  crosscorrelation  laws  are  the 
key  parameters  for  understanding  and  controlling  scattering  phenomena. 

Ill  -  SIMULTANEOUS  USE  OF  THEORY  AND  EXPERIMENT  IN  THE  CASE  OF  A  SINGLE 

SURFACE 

Let  us  first  specify  that  we  can  use  calculation  to  investigate  experiment  only  in  the 
case  of  surfaces  that  are  conformable  to  the  theoretical  model  (see  section  II),  that  are  of  good 
quality,  with  no  defects  such  as  dust  particles,  pits,  scratches  ...  The  experimental  results  that 
are  presented  here  are  always  relative  to  such  surfaces  that  have  been  previously  cleaned  and 
then  observed  with  a  Nomarsky  microscope.  We  intend  here  to  emphasize  the  problems  that 
must  be  solved  when  characterizing  an  optical  surface  [6,  23,  24]  with  scattering 
measurements,  and  to  bring  some  solutions.  We  will  limit  ourselves  to  the  case  of  sufficiently 
opaque  substrates  so  that  we  have  no  difficulties  with  volume  or  back  surface  substrate 
scattering.  In  the  case  of  transparent  sample,  we  usually  make  use  of  the  aluminum  technique 
[13,  25]. 

1)  Extraction  of  roughness  spectrum 

It  is  obvious  that  measurement  of  scattered  intensity  1(8,  <b)  together  with  calculation  of 
ideal  coefficient  a(9,  0)  leads  to  the  roughness  spectrum  y(S,  $).  We  have  verified  that  this 
roughness  spectrum  remains  invariable  when  the  incidence  and  polarization  of  incident  light 
are  modified  [15],  which  shows  that  y  is  a  real  characteristic  of  the  surface  defects.  Such 
results  are  therefore  a  good  test  for  both  the  validity  of  theory  and  the  accuracy  of  the 
measurements. 

In  many  cases,  we  verified  that  the  spectrum  y  of  isotropic  surfaces  could  be  fitted  in 
the  range  of  measurable  spatial  frequencies  with  the  Hankel  Transform  ya  of  the  sum  ra  of  an 
exponential  and  a  gaussian  function  [6]: 


r  (t)  =  Sg  e 


-2 

+  5.  e 


Fig.  5  gives  an  idea  of  the  quality  of  agreement  between  theory  and  experiment  such 
obtained  on  angular  scattering  curves. 


FigJ  -  Comparison 

of  theory  and  experiment 
for  four  black  glasses. 
BRDF.cosQ  [16]  is  the 
scattered  flux  per  unit  of 
surface  and  solid  angle, 
normalized  to  the 
incident  flux.  For  each 
sample  (a,  b,  c,  d),  we 
have  superimposed  the 
calculated  and  measured 
angular  scattering  curves 
at  quasi-normal 

incidence  and  at  the 
wavelength  2  ■  632.8  nm. 
We  have  also  plotted  the 
anisotropy  maps  [15] 
[2d]  for  3  samples.  The 
angular  range  (0  to  90°) 
corresponds  to  scattering 
by  reflection. 
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2)  Calculation  of  roughness 

The  roughness  5  (mean  quadratic  height)  can  immediately  be  obtained  by  integration  of 
the  roughness  spectrum  in  the  range  of  measurable  spatial  frequencies 


S2-  f  y(?)  d? 

For  an  isotropic  surface, 
with  k  =  and  Qm  = 


at  normal  incidence 

k 

we  obtain;  S2  =  In  f  <s  y  (<r)  da 
k  sin  9m 

scattering  measurement  minimum  angle  (with  respect  to  the  sample 


normal). 

However,  we  must  notice  here  (as  shown  in  section  11-2)  that  this  roughness  S  is 
relative  to  particular  defects  responsible  for  external  scattering,  that  have  spatial  periods  C  in 
the  range 

X  £  C  £  —  .  Consequently,  it  is  not  possible  to  definitely  associate  a  roughness  8  to  a  surface, 

since  the  value  of  5  will  depend  on  the  wavelength  under  study,  on  the  angular  range  where 
scattering  measurements  are  performed  and  on  the  illumination  incidence,  according  to  the 
different  kinds  of  defects  that  are  concerned.  Exactly  speaking,  a  roughness  measured  at 
wavelength  Xi  cannot  be  used  for  calculation  at  X2.  Therefore,  one  must  be  careful  as  far  as 
the  statistical  parameters  of  an  optical  surface  are  concerned. 


3)  Autocorrelation  length 

Usually,  the  autocorrelation  length  of  surface  defects  is  taken  as  the  value  Lg  of  the 

preceding  gaussian  function  (Le  >  L  g).  But  this  impiicitely  assumes  that  the  analytic  function 

ya  can  fit  the  roughness  spectrum  7  in  a  range  of  spatial  frequencies  from  0  to  infinity. 
However,  according  to  the  type  of  surface  under  study,  the  errors  such  committed  are 
practically  always  acceptable.  For  that  matter,  we  verified  that  the  exact  value  8  of  the 

roughness  (obtained  by  integration  of  the  roughness  spectrum  7)  is  very  close  to  the  value 

•\J  5e2  +  5.1 .  Anyway,  we  can  take  into  account  the  limited  domain  of  measurable  spatial 

frequencies8  and  obtain  the  autocorrelation  function  by  using  a  convolution  product  [241.  We 

have  also  studied  the  effects  due  to  the  apparatus  function  of  the  scatteromcter  [24], 

4)  Problems  connected  with  anisotropy  of  surface  defects 

Another  problem  which  appears  when  associating  a  roughness  to  an  optical 

micropolish  is  due  to  the  fact  that  anisotropy  defects  are  not  taken  into  account.  In  fact,  some 
surfaces  give  rise  to  scattering  levels  that  are  very  different  according  to  the  scattering  plane 
$  where  measurements  are  performed  (fig.  5).  Such  roughness  anisotropy  leads  to  some  error 
on  the  measurement  of  roughness  and  it  is  interesting  to  be  able  to  quantify  it.  For  that,  we 

calculate  an  anisotropy  function  F(a)  defined  by: 

F(a)  *  Jr  (t>)r[Rs(t>)ldT> 

* 


■  4«2  J  7  (c? )  7  (Ra  C#  )  !  d  # 
t 

where  T  is  the  autocorrelation  function  of  the  surface  defects.  While  r  ft* )  measures  the 
similitude  between  profiles  hfr* )  and  h  (r*  -  t*),  the  function  F(a)  measures  the  similitude 
between  r  (~? )  and  T  (Ra  )  1  .  It  is  the  angular  autocorrelation  of  the  autocorrelation 

function  r  (?)  of  the  surface  profile  h(r*).  In  the  case  of  isotropic  defects,  the  r  function  will 
be  radial  and  the  anisotropy  function  F(a)  will  be  constant.  Then  it  is  tonable  to  define  an 

isotropy  degree  from  the  minimum  of  the  curve  F(a)  (Fig.  6),  and  ate  surface  wiil  be  said 
perfectly  isotropic  if  its  isotropy  degree  is  equal  to  unity.  We  will  later  present  complementary 
results,  in  particular  about  the  variations  of  isotropy  degree  versus  spatial  frequency.  These 
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studies  concerning  roughness  anisotropy  are  quite  promiseful  since  anisotropy  maps  are 
real  signature  of  surface  defects. 


Fig.  6  - 
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IV  -  SIMULTANEOUS  USE  OF  THEORY  AND  EXPERIMENT  IN  THE  CASE  OF 

MULTILAYER  COATINGS 

The  case  of  multilayer  coatings  is  quite  more  complex  because  of  the  high  number  of 
parameters  involved  in  the  calculation.  We  must  take  into  account  the  roughnesses  5e.  5g  and 
autocorrelation  lengths  Le,  Lg  at  each  interface,  as  well  as  crosscorrelation  coefficients 
between  surfaces,  so  that  (Sp  +  4)  parameters  are  required  to  describe  scattering  from  a  p- 
layer  coating.  Though  it  is  always  possible,  in  these  conditions,  to  find  a  family  of  parameters 
leading  to  good  agreement  between  calculation  and  experiment,  this  family  is  not  necessarily 
unique  and  we  must  use  complementary  techniques  to  verify  the  obtained  conclusions.  For 
that,  the  use  of  specific  stacks  judiciously  chosen  [26]  or  particular  effects  such  as 
antiscattering  [17]  allows  us  to  eliminate  any  ambiguity  on  the  choice  of  these  parameters.  It 
is  also  possible  to  use  variations  of  scattering  versus  wavelength  [14,  25,  27]  or  illumination 
incidence,  or  the  aluminum  layer  technique  [13,  25],  to  confirm  the  results.  We  will  try  here 
to  rapidly  describe  these  different  methods  capable  of  pointing  out  the  key  parameters  of 
scattering.  We  will  also  recall  how  these  parameters  give  us  access  to  information  concerning 
the  microstructure  of  materials  in  thin  film  form. 

We  insist  again  on  the  fact  that  all  measurements  here  are  relative  to  thin  films  of  good 
quality,  so  that  calculation  can  be  meaningful.  Observations  with  a  Nomarski  microscope  allow 
us  to  verify  that  the  top  interface  of  the  coating  has  a  quality  very  similar  to  that  of  the 
substrate  before  coating  (no  dust,  pits  or  scratches  ...).  The  coatings  are  produced  with 

multiple  precautions  and  the  conclusions  drawn  here  obviously  depend  on  the  deposition 
techniques  that  are  employed.  In  order  to  obtain  a  detailed  confrontation  between  theory  and 
experiment,  the  coatings  are  generally  produced  on  black  glasses  and  scattering  from  back 
surface  substrate  is  therefore  eliminated. 

I)  Aluminum  technique 

Contrary  to  the  case  of  only  one  surface,  it  is  quasi  impossible  to  isolate  the  flux 
scattered  from  one  particular  surface  of  a  multilayer  stack.  However,  we  have  access  to  the 

roughness  of  the  top  interface  of  the  coating  by  covering  it  with  a  thin  (150  nm)  opaque  layer 
of  Aluminum.  Such  method  obviously  assumes  that  the  deposition  conditions  of  the  metal  are 
good  enough  so  that  all  top  interface  defects  are  reproduced  by  the  aluminum  layer.  On  this 
matter,  our  experimental  results  confirm  these  predictions,  at  least  for  glass  surfaces  with 

roughnesses  greater  than  0.5  nm.  The  mean  roughness  spectrum  has  been  found  to  be 

inchanged  before  and  after  deposition  of  aluminum,  in  the  whole  range  of  measurable  spatial 
frequencies  [25].  As  for  anisotropy  defects,  fig.  7  shows  that  the  results  are  extremely 
satisfying. 
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However,  one  should  take  care  when  using  this  technique  for  substrates  with 
roughnesses  less  than  5  A,  because  the  aluminum  microstructure  can  limit  the  reproduction  of 
substrate  defects  [28]. 


Fig.  7  -  Anisotropy 

functions  F(a)  of  a  black 
glass  before  and  after 
deposition  of  a  thin  opaque 
layer  of  Aluminum.  The  two 
curves  are  practically 
identical. 


Such  metallic  layer  technique  is  quite  useful  because  it  considerably  reduces  the 

number  of  scattering  parameters.  In  the  case  of  a  single  layer  for  instance  (9  parameters),  we 
can  measure  the  roughness  spectra  of  the  substrate  before  coating  and  of  the  top  interface 
after  coating  and  A1  deposition,  so  that  only  the  crosscorrelation  coefficient  remains  to  be 

determined.  However,  as  it  was  shown  by  few  experimental  results,  such  coefficient  a  can 
depend  on  spatial  frequency  [17]  and  can  lie  in  the  range  0  £  a  £  1  [29]. 

In  the  case  of  a  multilayer  mirror  (15  layers),  the  A1  layer  technique  led  to  excellent 

agreement  between  theory  and  experiment  [19,  25], 


2)  Use  of  antiscattering  effect 

Antiscattering  effect  has  already  been  largely  described  [17,  18],  and  is  related  to  the 
fact  that  scattering  can  strongly  be  reduced  after  deposition  of  one  particular  layer  on  a 
substrate  (provided  that  we  have  good  correlation  between  interfaces).  In  the  case  of  the  T1O2 
layer  in  fig.  8,  we  have  shown  [19]  that  reduction  of  scattering  cannot  be  explained 
assuming  uncorrelated  surfaces  (a  *  0):  the  two  scattering  interfaces  show  very  good 
correlation  (as  1). 


Fig.  8  -  Angular 
scattering  curves  measured  at 
quasi-normal  incidence  for  a 
black  glass  before  and  after 
coating  by  a  single  T1O2  layer 
with  optical  thickness  Xo  * 
632.8  nm.  Scattering  losses 
measured  at  Xo  are  reduced  by 
58%  after  coating. 


Antiscattering  effect  also  brings  information  on  the  microstructure  (at  the  scale  of 
optical  spatial  frequencies)  of  materials  in  thin  film  form.  In  the  case  of  the  Si02  layer  shown 
in  fig.  9  for  example,  the  two  angular  scattering  curves  before  and  after  coating  are 
identical,  which  indicates  that  the  two  interfaces  are  perfectly  correlated  (a  =  1)  and  that  they 
have  identical  roughnesses:  the  material  microstructure  is  sufficiently  fine  to  reproduce  ail 
substrate  defects  without  increasing  or  decreasing  their  magnitude.  Furthermore,  we  can  use 
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substrates  of  increasing  quality  until  we  obtain  departure  of  the  two  angular  scattering 
curves  before  and  after  coating,  which  will  be  connected  to  the  limit  of  reproduction  of 
substrate  defects.  Complementary  results  will  be  found  elsewhere  [17]. 


Fig.  9  -  Angular  scattering 
curves  measured  at  normal 
incidence  for  a  black  glass, 
before  and  after  coating  by  a 
S  i  O  2  layer  with  optical 
thickness  \q!  2.  The  two 
curves  are  practically 
identical,  which  indicates  a 
quasi-perfect  correlation  (a  s 
1 )  and  identical  roughnesses 
(  So  »  Si )  at  the  two  interfaces 
of  the  layer. 


3)  Variations  versus  wavelength 

Our  calculation  programs  enabled  us  to  show  that  scattering  losses  vary  in  phase  or  in 
phase  cancellation  with  the  reflection  coefficient  of  the  stack,  according  to  the  value  (  a  =  1  or 
a  *  0)  of  the  crosscorrelation  coefficient  between  interfaces  [27], 

The  example  shown  in  fig.  10  for  one  single  layer  obviously  shows  that  measurements 
of  scattering  versus  wavelength  will  rapidly  give  access  to  the  crosscorrelation  coefficients. 

Let  us  notice  also  the  stability  of  amiscattering  effect  (fig.  11)  over  a  wide  range  of 
wavelengths  in  the  case  of  a  quarterwave  low  index  layer. 


Fig  JO  -  Variations  of  scattering  losses 
(integrated  in  whole  space)  or  TIS  versus 
wavelength  in  the  case  of  ZnS  layer  with 
optical  thickness  Xp/2,  for  two  extreme 
values  of  the  crosscorrelation  coefficient  a. 
The  two  roughnesses  are  assumed  identical 
for  this  calculation. 


Fig  11  Variations  of  scattering  losses  (or 
TIS)  versus  wavelength  for  a  cryolite  layer 
with  optical  thickness  A0/4.  in  the  case  where 
the  ratio  of  the  two  roughnesses  is  equal  to 
(5q,/$1)  -  0.66  (So  it  the  roughness  of  the  top 
interface).  Amiscattering  effect  is  stable  over 
a  wide  range  of  wavelengths. 
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In  the  case  of  coatings  with  a  high  number  of  layers  (as  multipie  cavity  Fabry-Perot 
filters),  the  spectral  variations  of  scattering  [14]  are  however  quite  sensitive  versus 
scattering  parameters,  and  investigation  of  experimental  results  may  be  more  complex. 

4)  Use  of  specific  stacks 

Another  way  to  study  scattering  parameters  consists  in  using  particular  stacks  which 
scattering  losses  show  a  high  sensitivity  versus  scattering  parameters.  Analytical  calculations 
led  us  to  the  choice  of  a  design  made  of  an  odd  number  of  alternated  half-wave  high  index  and 
low  index  layers  [26].  As  shown  in  Table  I,  the  case  of  perfectly  correlated  surfaces  (ct  =  1) 
shows  no  variation  of  scattering  before  and  after  coating,  whatever  the  number  N  of  layers  of 
the  stack.  On  the  other  hand,  we  can  notice  a  high  increase  of  scattering  with  the  number  of 
layers  of  the  coating  (by  a  factor  100  with  11  layers)  in  the  case  of  uncorrelated  surfaces  (a  = 
0). 


Table  l  -  Scattering  losses 
D  by  r  efl  ection 


calculated  for  a  glass 
substrate  before  (N  =•  0) 

N-0 

N-1 

N-3 

N-5 

N-1 1 

and  after  coating  by  an 
odd  number  N  of 

106 

D  (a-0) 

23. 

649. 

1153. 

1660. 

3175. 

alternative  high  index 
(TiO  2 )  and  low  index 

0  (a-1) 

23. 

23. 

24. 

26. 

28. 

(SIO 2)  layers  with  optical 
thickness  XqI2. 

Calculation  is  performed  at  wavelength  Xo  for  two  extreme  values  (a  »  0  and  a  =  1)  of  the 
crosscorrelation  coefficient.  We  can  notice  that  for  correlated  surfaces  fa  *  1),  there  is 
practically  no  change  of  scattering  before  and  after  coating,  whatever  the  number  ,V  of 
alternative  layers.  On  the  other  hand,  for  uncorrelated  surfaces  (a  =»  0),  scattering  highly 
increases  with  the  number  N  of  layers  (a  factor  100  with  N  =  11). 

Some  experimental  results  are  presented  in  Table  II  and  show  that  scattering  is  always 
reduced  after  coating.  We  have  verified  [26]  that  this  can  only  be  explained  by  a  very  high 
correlation  (a  =  1)  between  interfaces  together  with  a  slight  reduction  of  roughness  at 

interfaces  (near  0.7). 


Table  II  -  Scattering  losses  of  3  black 
glasses  measured  before  (N  =  0)  and 
after  deposition  of  an  odd  number  of 
alternative  high  index  (7702)  an d 
low  index  (SiO 2)  layers  with  optical 
thicknesses  Xq!2.  Scattering  is 
always  reduced  after  coating. 


|  N-0 

N-1 

N-3 

N-5 

10®.  D 

41 . 

9.6 

10.0 

13.4 

V  -  BETTER  USE  OF  THE  MODEL 

We  have  seen  how  we  could  have  access  to  the  key  parameters  (roughnesses  and 

crosscorrelation  laws)  of  scattering  in  the  case  of  classical  coatings  and  consequently  obtain 
information  about  each  interface  inside  the  multilayer  stack.  The  results  obviously  depend  on 
the  quality  of  the  substrates  used  for  the  coating  as  well  as  on  the  deposition  techniques  that 

were  employed,  and  therefore  on  the  raicrostructure  of  the  materials  obtained  with  these 
techniques.  This  enables  us  to  modify  the  deposition  conditions  in  order  to  control  the 
scattering  parameters  and  obtain  minimal  scattering  losses  for  a  particular  design.  For 

example,  it  is  preferable  for  a  mirror  to  have  uncorrelated  roughnesses  [30]  while  it  is  the 
contrary  for  a  Fabry-Perot  filter  (in  this  last  case,  scattering  can  vary  by  a  factor  of  100 
according  to  the  value  of  the  crosscorrelation  coefficient).  Unfortunately,  as  far  as  designs 
with  a  high  number  of  layers  are  concerned  (multiple  cavity  Fabry-Perot  filters  for  example), 
the  preceding  specific  techniques  of  comparison  between  experiment  and  theory  are  not 
sufficient  enough  to  investigate  experimental  results  since  the  number  of  scattering 
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parameters  is  too  high.  It  is  difficult  for  instance  to  point  out  the  origin  of  roughness  in 
multilayer  stacks.  We  have  shown  [13,  30]  in  some  cases  that  such  roughness  could  be  due 
both  to  reproduction  of  substrate  defects  (for  low  spatial  frequencies)  and  to  a  residual 
roughness  (for  high  spatial  frequencies)  brought  by  the  material  microstructure. 

We  intend  here  to  briefly  show  how  the  theoretical  model  can  be  improved  to  reduce  the 
number  of  parameters  and  to  permit  an  easy  investigation  of  scattering  phenomena  in 
complex  coating  designs. 


1)  Influence  of  substrate  defects 

Our  experimental  results  show  that  in  many  cases,  the  anisotropy  maps  of  the  substrates 
are  very  similar  before  and  after  coating  (and  after  deposition  of  Al). 

In  the  easy  case  of  a  single  layer,  we  even  show  (see  section  IV-2)  that  the  two  interface 

roughness  spectra  are  identical  in  the  whole  range  of  measurable  spatial  frequencies.  This 
leads  us  to  consider  that  'here  is  a  relation  between  cause  and  effect  [31]  for  the  surface 

profiles  inside  the  multilayer,  relation  which  will  of  course  depend  on  the  quality  of  the 
substrate  [29],  on  the  microstructure  [171  and  thickness  [29]  of  material  layer,  and  also  on 

the  spatial  frequency  [17]  of  the  concerned  defect. 

For  instance,  we  can  write  the  following  relation: 

h>  =  ^2  a*j  *  hi 

where  h,  and  hj  are  profiles  of  interfaces  (i)  and  (j),  and  (*)  indicates  a  convolution  product. 
We  then  obtain  after  Fourier  Transform  :  hiC<?)  =  aij(^)  hj(c?)  •  l  a;j  =  FT  [a;j  (7*)]. 
Consequently,  the  crosscorrclation  coefficient  otjj  represents  the  ratio  of  the  magnitudes  of 

the  gratings  of  same  period  C  *  —  at  interfaces  (i)  and  (j)  [19].  It  varies  with  spatial 

frequency  a  so  that  we  can  consider  the  fact  that  the  material  will  act  differently  according  to 
the  kinds  of  defects  that  are  concerned. 

With  this  "causal"  model,  the  number  of  free  parameters  is  considerably  reduced 
because  roughness  spectra  and  crosscorrelation  laws  are  henceforth  linked  by  the  relation: 


aijC?) 


VI 


iC?)  J^iiCg) 

TjjC?)  6 


Moreover,  the  ajj  coefficients  can  be  calculated  with  the  ajj+i  coefficients.  As  for  the 

argument  'Fij('<?),  one  can  consider  that  it  is  equal  to  zero  provided  that  deposition  of  material 
is  made  near  normal  incidence  and  on  rotating  substrates  (gratings  are  reproduced  with  no 
lateral  shift).  In  the  case  of  a  quarterwave  coating  made  of  only  two  materials,  we  can  also 
assume  that  each  material  will  have  the  same  behavior  whatever  the  layer  of  the  stack  which 
is  made  of  this  material.  Therefore,  we  have  to  consider  only  two  functions  a, j  that 

characterize  the  two  materials:  cti,i+i  =  cti+2,i+3- 

In  that  way,  provided  that  the  substrate  was  measured  before  coating,  the  unknown 
parameters  are  the  crosscorrelation  functions  aoi  and  ai2(c?).  As  it  is  probable  that  the 

greater  the  defect  periods  the  better  they  are  reproduced,  we  can  choose  for  otjjC^)  a  slow 
decreasing  function  analogous  to  a  low-pass  filter  (or  other)  with  cut-off  frequency  <jjjc.  We 

will  not  enter  here  the  detail  of  calculation,  but  we  will  emphasize  the  fact  that  such  causal 
model  allows  computation  for  systematic  research  of  the  parameters  that  show  good  agreement 
between  theory  and  experiment.  A  good  agreement  will  be  obtained  in  that  way  provided  that 
the  materials  do  not  bring  an  important  residual  roughness  (due  to  their  microstructure). 

Figure  12  illustrates  these  results  in  the  case  of  a  Fabry-Perot  filter.  The  two  extreme 
cases  (a  a  0  and  a  a  1)  that  are  generally  used  are  also  represented.  With  our  causal  model  and  a 

cut-off  frequency  oc  ■  =  12.6  4m’1,  the  angular  scattering  curve  is  very  similar  to  one  and 

*-c 

other  classical  cases  of  uncorrelated  (a  a  0)  and  correlated  (a  al)  surfaces,  according  to  the 

angular  range  that  is  concerned.  This  is  explained  by  the  fact  that  scattering  at  low  angles  is 

due  to  large  period  defects  (<t  «  <rc)  that  are  reproduced  (as  1),  while  scattering  at  large  angles 

is  due  to  short  period  defects  that  are  lost  (a(o)  <  l).  Let  us  notice  here  that  the  case  of 
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homothetic  surfaces  (see  II-3)  would  lead  to  the  easy  relation  ay  =  which  recalls  that 

roughnesses  and  crosscorrelation  laws  are  not  independent.  On  the  other  hand,  the  classical 
case  of  independent  roughnesses  is  not  taken  into  account  with  this  model  (see  section  V-3). 

Fig  12  -  Angular  scattering  curve 

calculated  at  the  wavelength  Xq  and  BRSF'octftUtnj _ 

at  normal  incidence  for  a  Fabry-  j  j  j  •  ■  '  •  '  ’ 

Perot  filter  which  design  is:  H  L  H  \  I  ;  ;  1  1  1  !  ;  ; 

L  H  (6L)  H  L  H  L  H  where  H  and  L  \\  j  ,J  \  j  j  j!  [  \  \  /  1 

are  quarterwave  Xq!4  high  index  ;  ;  ;  i, 

(ZnS)  and  low  index  ( cryolite )  \j  ;  /  J  i\  i  ;  St";  \  ;  'J  j 

layers.  We  can  notice  the  presence  fj  yj  \  j\  L  j  ;  \  j  y<  ,1 

of  scattering  peaks  [14,  25].  The  "  W_j  Wf  /  \  \\:  :  JVVS.  / ! 

two  extreme  cases  a  »  0  and  a  =■  1  y  '-y  j  Y-.l  '//j  ! 

that  are  represented  correspond  to  f f  *  :  :  \T*.  ■;/  1  :  i  I  ' 

the  classical  cases  of  correlated  and  j  j  j  vjT  ;f.l  i  i  |i  ’  - 

uncorrelated  surfaces.  The  curve  in  ;  ;  ;  f  i  /  ;  i  ;  wa->. 

full  line  has  been  calculated  t j  0J _  j . ■ .  j .  ...  ■  \  t  •/  j  '•  ■  —  aiphi-o. 

according  to  the  causal  model  and  1  1  ;  tj  f  1  ;  ;  _  ^ootm 

assuming  that  the  crosscorrelation  j  ;  ;  i  y/i  i  I  ;  * 

coefficient  a(  a)  is  a  slow  decreasing  ;  !  j  '■  <>  ■  '  '■  ■  ; 

low-pass  filter  with  a  cut-off  a  to  to  so  so  too  iza  tio  tso  iso 

,  2 it  ,  ,  Tluta 

frequency  ac  m  T~  =*  1 2.0  pm'1 . 

Lc 

The  angular  range  (90°  -*  180°)  corresponds  to  scattering  by  transmission. 


—  sip/a*(. 

—  aip/ut-0. 

—  Ic-  SGQxn  I 


too  120  140  ISO  ISO 

nuto 


2)  "Pseudo-anisotropy”  effects 

In  spite  of  experimental  and  theoretical  difficulties,  we  can  go  further  to  know  until 
where  the  preceding  causal  model  will  be  realistic.  Let  us  consider  for  instance  a  deposition 
performed  at  oblique  incidence  p  on  a  non-rotating  substrate.  As  the  material  flux  reaches  the 
substrate  at  incidence  P,  the  layer  will  grow  under  certain  conditions  in  a  form  of  columns 
with  an  angle  a  with  respect  to  the  sample  normal  (tg  p  =  2tg  a  [32]). 

Therefore  we  can  consider  that  the  different  gratings  are  reproduced  with  a  lateral 
shift  (which  modulus  is  I  a*  I )  depending  on  the  incidence  p  and  on  the  layer  thickness.  In  this 
case,  the  crosscorrelation  coefficient  a  is  not  a  real  number  and  we  must  take  into  account  its 
argument  *P  jj(of)  = .  As  shown  in  fig.  13.  we  then  observe  a  " pseudo- anisotropy"  effect 
which  means  that  scattering  varies  with  the  scattering  plane  0  that  is  concerned  (though  the 
surfaces  are  isotropic). 

Experiments  are  under  study  but  are 
quite  difficult  because  of  the  high 
thickness  of  the  layer;  indeed  we 
must  not  lost  the  relation  between 
cause  and  effect  for  the  surface 
profiles.  Some  results  will  be  given 
later. 

Fig.  13  -  Angular  scattering  curves 
calculated  at  waveleng •‘t  Xq  and  at 
normal  incidence  for  a  cryolite 
layer  with  optical  thickness  13 Xq!  4 
deposited  at  incidence  p  *  25°  on  an 
isotropic  and  non  rotating 
substrate.  The  crosscorrelation 
coefficient  a  is  therefore  a  complex 
number  ( V  *  0)  and  it  appears 
" pseudo-anisotropy "  effects,  which 

means  that  scattering  varies  with  0  a  40  SO  SO  too  m  uo  ISO  iso 

the  scattering  plane  0. 
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3)  Material  grain  size 

Until  now,  we  considered  only  the  reproduction  of  substrate  defects  layer  after  layer, 
without  taking  into  account  a  residual  roughness  due  to  the  material  microstructure  [13,  17], 
However  it  is  possible  to  consider  both  these  two  effects  by  the  relation: 

&*i+l  ai+l.i  *  h*i  +  Si-*- 1 

where  the  first  term  is  relative  to  reproduction  of  defects  and  the  second  term  g;+i  to  the 
residual  roughness  brought  y  the  material.  The  notation  h*;  is  used  here  to  indicate  the 
difference  with  the  preceding  causal  model. 

We  then  obtain: 

i 

Tii*  =Tii+  X  ®ik«iqakqs  i  Sq ! 2 

(k.q) 


r.j*  =  nj  + 


X  «ik  c£ jq  OkqS 
k=l  to  i 
q=i  l  to  j  ' 


where  otjjS  indicates  a  possible  correlation  between  residual  roughnesses  at  interfaces  (i)  and 
(j).  In  the  case  of  a  coating  made  of  only  two  materials,  we  can  assume  that  the  residual 

roughnesses  of  the  two  materials  are  independent,  so  that  a  q,q+2k+i  =  0,  and  therefore  easy 
computation  can  be  performed. 

VI  -  MINIMAL  SCATTERING  LOSSES 

The  techniques  presented  in  the  preceding  sections  will  permit  us,  owing  to  a  better 

understanding  of  scattering  phenomena,  to  adjust  the  deposition  conditions  in  order  -  to  obtain 
minimal  scattering  losses.  However,  we  could  also  try  to  modify  the  coating  design  so  that 

scattering  losses  are  reduced,  without  altering  the  required  optical  properties.  This  is  quite  a 
difficult  problem,  and  we  will  limit  ourselves  here  to  the  case  of  a  simple  example.  For  a 

classical  mirror,  it  has  been  shown  that  it  was  possible  to  shift  the  ideal  electric  field  in  the 

volume  of  the  stack  by  addition  of  particular  layers  on  the  mirror,  without  decreasing  the 

reflection  coefficient.  In  this  way,  absorption  [33]  and  scattering  [30]  can  be  reduced. 
However,  reduction  of  scattering  is  only  50%  and  this  only  occurs  in  the  case  of  unconelated 
surfaces.  This  method  can  also  be  applied  to  the  case  of  Fabry-Perot  filters  (Fig.  14)  and 
scattering  is  then  reduced  by  a  factor  5,  for  uncorrelated  surfaces  (let  us  notice  that  the  case  a 

=  0  can  occur  with  substrates  of  high  quality).  However,  these  results  are  somewhat  limited, 

and  it  is  rather  useful  to  search  for  destructive  interferences  between  the  scattered  waves,  as 
shown  in  section  IV-2  and  IV -4. 


Fig.  14  -  Distribution  of  the 
square  of  the  electric  field  in  a 
Fabry-Perot  filter  which  design 

is: 

HL  HL  HU  H'U  H'  (2 L\  H’  L'H' 
L'H  LH  LH  The  design  is  slightly 
different  with  respect  to  that  of 
classical  quurterwave  Fabry- 
Perot  filters  so  that  the  maxima 
of  electric  field  [30,  33 J  inside 
the  stack  are  shifted  in  the 
volume  of  the  layers.  Surfacic 
scattering  is  therefore  reduced 
by  a  factor  of  5  in  the  case  of 
uncorrelated  surfaces. 
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On  the  other  hand,  in 
particular  for  multi-demultiplexing 
in  optical  telecommunications  [14], 
we  often  have  the  choice  of  several 
designs  that  can  reach  the  required 
isolation  and  transmission 
specifications.  Our  theoretical 
calculations  then  enable  us  to 
choose  the  multiple  cavity  Fabry- 
Perot  filter  which  leads  to  minimal 
scattering  at  critical  directions  of 
transmission  information  (Fig.  15) . 


Fig.  15 


-  Scattering  distribution  from  a  Fabry-Perot 
filter. 


CONCLUSION 

We  presented  a  rapid  overview  of  the  works  that  have  been  performed  in  our 

Laboratory  for  the  study  of  scattering  from  multilayer  optics.  We  expect  promiseful  results 

from  other  experimental  and  theoretical  works  under  study. 

The  case  of  an  uncoated  surface  appears  to  be  well  solved  since  we  can  extract  the 
roughness,  autocorrelation  length  and  isotropy  degree  in  the  range  of  measurable  spatial 
frequencies.  For  multilayer  classical  coatings,  we  have  shown  how  it  was  possible  to  extract 
the  key  parameters  of  scattering  with  the  help  of  particular  techniques  such  as  A1 

overcoating,  antiscattering  effect,  variations  of  scattering  versus  wavelength,  and  specific 
stacks.  These  parameters  (roughnesses  and  crosscorrelation  laws)  give  us  precious 
information  concerning  the  materials  microstructure,  and  it  is  therefore  possible  to  modify 
the  deposition  conditions  to  obtain  the  finest  microstructure. 

In  the  case  of  complex  coating  design,  the  high  number  of  parameters  causes  the 
comparison  between  theory  and  experiment  to  be  extremely  difficult.  However,  we 
demonstrated  that  a  more  realistic  model  for  crosscorTelation  laws  permitted  us  to 

considerably  reduce  the  number  of  parameters  involved  in  the  calculation;  one  can  also 
search  for  pseudo-anisotropy  effects.  In  a  more  general  way,  these  results  allow  computation 
for  systematic  research  of  the  scattering  parameters  that  lead  to  the  best  agreement  between 
theory  and  experiment,  and  therefore  permit  an  easy  understanding  of  scattering 
mechanisms.  At  last,  such  calculation  enables  us  to  choose,  at  the  conception  of  a  multi¬ 
demultiplexing  system,  the  coating  designs  that  lead  to  minimal  scattering  at  the  critical 
directions  where  information  is  transmitted. 
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Comparison  of  the  properties  of  titanium  dioxide  films 
prepared  by  various  techniques 


Jean  M.  Bennett,  Emile  Pelletier,  G.  Albrand,  J.  P.  Borgogno,  B.  Lazarides,  Charles  K.  Camiglia,  R.  A. 
Schmell,  Thomas  H.  Allen,  Trudy  Tuttle-Hart,  Karl  H.  Guenther,  and  Andreas  Saxer 


Fourteen  university,  government,  and  industrial  laboratories  prepared  a  total  of  twenty  pairs  of  single-layer 
titanium  dioxide  films.  Several  laboratories  analysed  the  coatings  to  determine  their  optical  properties, 
thickness,  surface  roughness,  absorption,  wetting  contact  angle,  and  crystalline  structure.  Wide  variations 
were  found  in  the  optical  and  physical  properties  of  the  films,  even  among  films  produced  by  nominally  the 
same  deposition  techniques. 


I.  Introduction 

This  paper  began  as  a  project  of  the  Optical  Society’s 
Optical  Materials  and  Thin  Films  (OM&TF)  Techni¬ 
cal  Group.  The  question  was  asked  “What  are  the 
optical  constants  of  a  titanium  dioxide  thin  film?”  To 
answer  it,  fourteen  university,  government,  and  indus¬ 
trial  laboratories  made  single  layer  TiOn  films  using 
seven  different  processes:  electron-beam  (E-B) 
evaporation,  ion-assisted  deposition  (IAD),  ion-beam 
sputter  deposition  (IBSD),  activated-reactive  evapo¬ 
ration  (ARE),  ion  or  plasma  plating  (IP),  rf  diode 
sputtering  (RFS),  and  dip  coating  from  a  solution  us¬ 
ing  the  Schott  procesi  Dip).1  Commercial  grade  pol¬ 
ished  fused  silica  substrates  were  provided  to  each 
group.  Two  coated  samples  were  made  of  each  type  of 
film  and  the  complete  set  of  samples  was  then  divided 
into  two  groups  for  measurement  and  analysis. 

Initially,  the  optical  properties  of  the  films  were 
analyzed:  the  index  of  refraction,  absorption  and  in¬ 
homogeneity  as  a  function  of  wavelength,  along  with 
the  thickness  as  determined  by  spectrophotometry. 
One  set  of  samples  was  measured  and  analyzed  at 
Ecn'e  Nationale  Superieure  de  Physique  de  Marseille. 
This  set  of  samples  and  the  group  at  Marseille  that 
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performed  the  analysis  are  both  designated  as  Group 
A.  Measurements  were  made  on  the  second  set  of 
samples  by  the  Research  Department  at  Optical  Coat¬ 
ing  Laboratory,  Inc.  (OCLI),  and  the  data  were  ana¬ 
lyzed  at  the  Developmental  Optics  Facility  of  the  U.S. 
Air  Force  operated  by  Martin  Marietta  Astronautics 
Group  in  Albuquerque,  N.M.  This  second  set  of  sam¬ 
ples  and  the  groups  at  OCLI  and  Martin  Marietta  are 
both  called  Group  B. 

The  results  of  the  optical  properties  measurements 
of  twenty  different  titania  films  were  reported  by  each 
group  at  the  OM&TF  Technical  Group  Meeting  at  the 
1986  Annual  Meeting  of  the  OSA.  The  optical  con¬ 
stants  of  different  films  showed  some  correlation  with 
deposition  method,  but  also  indicated  that  there  was 
considerable  influence  of  the  coating  system,  coating 
conditions,  and  operator. 

During  the  following  year  the  films  in  Group  B  were 
subjected  to  a  number  of  additional  measurements. 
These  measurements  included  surface  roughness,  ab¬ 
sorption  (using  photothermal  deflection  spectrosco¬ 
py),  wetting  contact  angle,  and  film  microstructure  as 
determined  by  Raman  spectroscopy.  The  results  were 
reported  at  the  1987  Annual  Meeting  of  the  OSA. 

This  paper  reports  on  the  entire  project,  including 
the  film  preparation,  the  measurements  and  analysis 
techniques,  and  the  results.  We  found  that  the  optical 
and  physical  properties  of  various  films  of  the  same 
material  were,  in  fact,  quite  different.  Thus,  the  ini¬ 
tial  question  about  the  optical  constants  of  a  titanium 
dioxide  film  does  not  have  a  unique  answer. 

A  related  question  might  inquire  about  which  depo¬ 
sition  process  produced  the  best  films.  The  results  of 
the  study  indicate  that  the  answer  to  this  question 
depends  on  the  -riteria  that  are  used  to  define  the 
term,  best:  high  refractive  index,  low  absorption,  low 
scatter,  etc.  We  conclude  that  there  is  enough  varia¬ 
tion  in  the  observed  properties  that  titania  films  can  be 
optimized  for  a  wide  range  of  applications. 
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D.  Film  Preparation 

The  substrates  provided  to  all  participants  were  25.4 
mm  (1  in)  diam  by  1.6  mm  (0.063  in.)  thick  polished 
fused  silica  substrates  of  the  type  frequently  used  for 
commercial  coatings.  Highly  polished  surfaces  (bowl 
feed  polished,  for  example)  would  have  been  prefera¬ 
ble  for  the  evaluation  of  some  properties  such  as  sur¬ 
face  roughness  and  scatter.  However,  the  initial  in¬ 
tent  of  the  project  was  to  look  at  the  optical  properties 
of  typical  coatings  rather  than  the  very  best  ones  that 
could  be  made  under  ideal  laboratory  conditions. 
(Also  the  cost  of  super  smooth  substrates  was  about  ten 
times  higher.)  Even  though  the  substrate  roughness 
was  non-negligible,  it  was  possible  in  most  cases  to 
measure  the  added  roughness  from  the  films.  The 
starting  material  for  the  films  was  supplied  by  each 
participating  group,  since  each  process  required  differ¬ 
ent  forms  of  titanium  oxide  or  titanium  metal. 

Each  participant  provided  two  samples  coated  in  the 
same  run  with  a  nominal  optical  thickness  of  5  quarter 
waves  at  550  nm.  This  thickness  was  chosen  so  that  a 
spectral  transmittance  scan  over  the  region  of  trans¬ 
parency  of  the  titanium  film-silica  substrate  combina¬ 
tion  would  contain  a  sufficient  number  of  maxima  and 
minima  to  allow  the  refractive  index  to  be  determined 
accurately.  The  participants  had  been  told,  however, 
that  the  actual  thickness  was  relatively  unimportant. 
As  is  shown  in  Sec.  IV,  the  average  wavelength  for 
which  the  films  had  an  optical  thickness  of  five  quarter 
waves  was  close  to  570  nm. 

The  participating  groups  and  the  types  of  films  they 
made  are  listed  in  Table  I.  They  had  all  been  asked  to 
provide  the  parameters  of  the  processes  they  were 
using  if  the  information  was  not  company  proprietary. 
These  are  listed  in  Table  II. 

Most  of  the  coating  processes  are  well  known.  How¬ 
ever,  a  brief  summary  is  given  here.  Excellent  reviews 
of  most  of  the  ion  assisted  processes  have  been  given  by 
Martin2  and  Matthews.3  There  are  additional  excel¬ 
lent  review  articles  on  ion  beam  sputter  deposition,4 
reactive  ion  plating,5-6  and  dip  coating.7  Optical  prop¬ 
erties  of  titania  films  prepared  by  different  techniques 
are  summarized  by  Martin;2  the  improvement  in  the 
optical  properties  by  intermittent  ion  bombardment 
has  also  been  discussed.8-9 

Electron-Beam  (E-B)  Deposition:  The  substrates 
are  placed  in  a  vacuum  chamber  with  a  base  pressure  in 
the  10-6  mbar  range.  The  source  is  an  electron  beam 
gun  into  which  titanium  oxide  i3  placed.  A  focused 
beam  of  electrons  in  the  gun  melts  the  material;  it 
evaporates  onto  the  substrates  which  have  been  heated 
to  about  300  °C.  n'  improve  the  stoichiometry,  the 
deposition  is  ger  .ly  done  in  a  partial  pressure  of 
oxygen  of  0.6-3  X  10~4  mbar.10  Typical  evaporation 
rates  are  0.2-0.8  nm/s  Films  deposited  using  an  elec¬ 
tron  beam  gun  or  other  thermal  sources  have  a  colum¬ 
nar  raicrostructure11  which  can  be  porous.  The  result¬ 
ing  surface  roughness  and  grain  boundaries  produce 
light  scattering.12  For  this  reason,  other  techniques 
have  been  developed  to  make  denser  films  whose  prop¬ 
erties  are  closer  to  that  of  the  bulk  material. 


Table  (.  Suppliers  of  Titania  Films 
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Group 

IP 

163.166 

Baleen  AG.  FL-9496  3a leers.  Liechtenstein.  W  Thncni. 

M.  (Caspar 

RFS 

093.145 

fianclte.  Pacific  Northwest  Division.  OailclIe  B!vd..  PO  Box  ‘799. 
Richland,  WA  99352.  W.  T.  Pawlcw.cz 

IBSD 

078.079 

Colorado  State  University.  Department  of  Physics.  Fon  Collins.  CO 
S0523.  James  R-  Sites,  Rons  Rujkorakam 

Dip 

21 3a  21 3b 

Denton  Vacuum.  Inc..  Cherry  Hill  Industrial  Center.  Cherry  1  [ill. 

S)  08003.  Peter  R.  Denton.  Anthony  Musset 

E-B.ARE 

052.057 

Indian  Institute  of  Science.  Instrumentation  and  Services  Unit. 
Bangalore  560  012.  India.  S  Mohan 

E-B 

025.023 

Leybold  AG.  Dept.  UC2,  Sicmenssrr.  100,  8755  Alienau.  West 

IAD 

024,022 

Germany.  R.  Herrmann.  W  Klug 

E-B 

094.097 

Martin  Maneita  Astronautics  Group.  Later  Systems  Technology. 
P.O.  Boa  9316.  International  Airport.  Albuquerque.  NM  87 1 19. 

R.  A.  Schmelf 

E-B 

201.202 

National  institute  of  Standards  and  Technology  (formerly  National 
Bureau  of  Standards).  Bldg.  223.  Room  A329.  Gaithersburg.  MD 
20899.  Albert  Feldman.  Ed  Farabaugh 

IBSD 

120.122 

Optical  Coating  Laboratory.  Inc..  2789  Korthpomi  PVwy.  Santa 

IP 

114.115 

Rosa.  CA  95407.  T  H.  Allen.  T  Tuttle-Han 

E-B 

138.140 

Optical  Sciences  Center.  Umv  of  Arizona.  Tucson.  AZ  85721. 

IAD 

137.139 

H.  A.  Macleod.  John  Lehan 

lAD 

143.144 

E-B 

036.038 

Opto  Mechandc.  fnc..  425  Noah  Dave.  PO  Bok  J6I907, 

Melbourne.  FL  32935.  Ron  R  Willey 

IBSD 

135.136 

Rockwell  International.  Electronics  Operations.  Auionettcs  Marine 
Systems  Division.  3370  Miraloma  Ave..  Mail  Sta.  GE  14. 

Anaheim.  CA  92803.  Austin  N  Kalb.*  Mark  Mildebrath 

E-B 

003.002 

Spectra  Phvt.cs,  (nr.  Optics  Division.  1250  West  MiddlefieJd  Rd . 

E-B 

006.005 

Mountain  View,  CA  94Q42.  Gary  DeBeli 

IAD 

013.014 

University  of  New  Mexico.  Dept  of  Electrical  and  Computer  Engr.. 

IAD 

011.012 

Albuquerque.  N.  Mex.  87131.  J.  R.  McNeil.  Forrejt  Williams 

E-B  Electron  Bom  Deposition  RFS  RF  Dtotie  Sputtering 

IAD  Ion  A iti tied  Deposition  IP  Ion  (or  Plasma)  Plating 

I0SD  ion  Beam  Sputter  Deposition  Dip  Schott  Dip  Process 

ARE  Activated  Reactive  Evaporation 
*  Present  address:  4711  Reeves  Road.  Ojat.  CA  93023 


Table  II.  Coating  Parameters  tor  Titania  Films 


Type 

Sample  » 

Base 

Temp 

Oxvgen 

Rate 

E  gun 

Ion  gun 

Notes 

Pressure 

Pressure 

|10*!mbar 

■c 

[tO'lmbar 

"TVS 

kV 

nvA 

V 

mA 

E-B 

003.002 

006.005 

201.202 

1 

300 

05 

02 

10  5 

a 

b 

c 

025.023 

13 

300 

2.6 

02 

84 

138.140 

26 

136 

26 

04 

052 

1.3 

75 

06 

02 

95 

006 

d 

094.097 

4 

225 

1.3 

0  25 

65 

0  28 

036.038 

26 

300 

26 

03 

70 

004 

IAD 

024.022 

13 

50 

2  6 

016 

84 

750 

100 

e 

137.139 

2.6 

136 

26 

04 

300 

32 

f 

U3. 144 

2.6 

137 

i.6 

04 

300 

30 

f 

013.014 

008 

175 

2 

0  15 

64 

0.16 

300 

0JJ.0J2 

008 

275 

1,3 

015 

64 

0  16 

500 

h 

IBSD 

135.136 

00) 

26 

0  093 

1300 

MO 

, 

078.079 

06 

50 

06 

0025 

1 100 

30 

120.122 

2.6 

100 

12 

1400 

175 

k 

ARE 

057 

1.3 

75 

06 

0.1 

95 

006 

d 

RFS 

093.145 

0  1 

250 

40 

01 

U50 

I 

IP 

165.166 

114.115 

Dip 

2l3a,b 

Notes 

j  Smalt  chamber 
b  Large  chamber 
c  Ion  cleaned  (t/bimres 
J  Crvopumned  system 
e  50uA/crn, 
f  Intermittent 
g  JO  wA/cm* 


h  ■lOuA/em* 

.  VUin:  beam 

;  Argon  hi  lO'm.Mf  ihnWnetS 
»  jnjiii'n  noted  m  films 
i.  IAD  50  V.  200  mA 
I  JW/cmv  Argon  1  )  i  10 :  mbar. 
substrate  bus  voltage  i(C  V 


3304  APPLIED  OPTICS  /  Vol.  28.  No.  15  /  15  August  1969 


Ion-Assisted  Deposition  (I AD):  A  beam  of  ions, 
generally  argon,  is  directed  onto  the  substrate  along 
with  the  arriving  titania  molecules.  The  ion  beam  is 
usually  produced  by  a  Kaufmann-type  ion  source. 
The  ions  break  up  the  normal  columnar  growth  and 
densify  the  film  by  impulse  transfer.13  Films  pro¬ 
duced  by  IAD  have  smoother  surfaces  and  show  less 
grain  structure.  Hence  they  have  lower  scattering 
levels  than  films  produced  by  electron  beam  deposi¬ 
tion. 

Ion-Beam  Sputter  Deposition  (IBSD):  A  beam  of 
energetic  argon  ions  is  used  to  sputter  material  from  a 
titanium  metal  or  titanium  oxide  target  onto  the  sub¬ 
strates.  The  argon  beam  is  produced  by  admitting 
argon  gas  into  an  ion  gun  at  a  pressure  of  about  10~2 
mbar.  Electrons  from  a  thermionic  emitter  ionize  the 
gas.  Argon  ions  from  the  source  (having  energies  in 
the  1000  to  1400  eV  range)  are  extracted  into  the 
coating  chamber  using  dual  grids  to  form  a  well-de¬ 
fined  ion  beam  which  is  directed  at  the  sputtering 
target.  The  coating  chamber  is  differentially  pumped 
to  a  back  pressure  in  the  range  of  10-4  to  5  X  10~5  mbar. 
Oxygen  is  usually  also  admitted  at  a  partial  pressure  of 
1-2  X  10"4  mbar  to  react  with  the  sputtered  material. 
In  dual  ion  beam  sputtering,  the  second  ion  source 
bombards  the  growing  titania  film  with  low-energy 
oxygen  ions  to  control  its  stoichiometry  or  for  sub¬ 
strate  precleaning.  The  low  energy  of  the  oxygen  ions 
(50-100  eV)  minimizes  coating  damage.  Films  pro¬ 
duced  by  IBSD  are  dense  and  homogeneous;  film  den¬ 
sity  can  be  close  to  that  of  the  bulk  material. 

Activated  Reactive  Evaporation  (ARE):  Titania 
films  are  made  using  an  electron  beam  gun  to  evapo¬ 
rate  titanium  oxide  onto  the  substrates  in  a  reactive 
gas  atmosphere.  An  additional  discharge  source  is 
used  through  which  the  reactive  gas  passes.  The  reac¬ 
tive  gas  is  ionized  or  excited  and  emerges  from  a  nozzle. 
It  then  impinges  on  the  substrate  along  with  the  coat¬ 
ing  material  in  vapor  form.  The  deposition  in  the 
presence  of  ionized  and  excited  gas  enhances  the  reac¬ 
tivity,  thereby  ensuring  complete  oxidation  and  im¬ 
proving  the  stoichiometry  of  the  films.  As  a  result,  the 
films  are  nearly  absorption  free.  The  energy  of  the 
ions  is  so  low  that  it  causes  minimum  disturbance 
(dissociation)  to  the  growing  film. 

Ion  Plating  or  Plasma  Plating  (IP):  In  this  process 
titanium  oxide  or  titanium  metal  is  melted  in  an  elec¬ 
tron  beam  evaporator.  The  E-B  crucible  serves  as  the 
anode  of  a  low-voltage,  high-current  plasma  arc  which 
is  maintained  with  argon  gas  at  about  2  mbar  pressure 
in  the  source  compartment.  A  hot  cathode  ionizes  the 
argon  gas.  The  pressure  differential  at  a  1-2  mm 
diam.  hole  leading  to  the  main  vacuum  chamber  (at  a 
background  pressure  below  10'8  mbar)  extracts  the 
argon  plasma.  Oxygen  is  admitted  at  a  pressure  of 
~10"3rnbar  to  react  with  the  titanium  metal  ions  at  the 
substrate  surface  and  produce  titanium  dioxide.  The 
large  electron  current  formed  in  the  plasma  discharge 
not  only  ionizes  the  vapor  above  the  crucible  (M+, 
MO*)  but  also  the  reactive  gas  in  the  main  chamber 
(Oj).  The  plasma  sheath  thus  produced  in  the  main 


chamber  is  in  contact  with  the  dielectric  substrates 
which  sit  on  isolated  holders.  Electrons  leaving  the 
plasma  more  easily  than  the  heavier  ions  charge  the 
substrates  negatively  with  a  self  bias  of  -5  V  to  -60  V. 
This  bias  attracts  and  accelerates  the  positive  ions 
(titanium  vapor,  oxygen  and  argon)  from  the  plasma 
environment.  The  impact  of  ions  in  this  energy  range 
densities  the  growing  film  by  a  momentum  transfer 
process.14  Resputtering  does  not  seem  to  occur  except 
possibly  as  preferential  sputtering  from  asperities  pro¬ 
truding  from  the  film  surface  and  of  coadsorbed  impu¬ 
rities.  This  process  helps  to  make  the  film  smooth  and 
eliminate  impurities  such  as  water  vapor. 

RF  Diode  Sputtering  (RFS):  An  rf  sputtering  tech¬ 
nique  has  been  used  at  Battelle  Pacific  Northwest 
Laboratories18  for  depositing  metal  and  dielectric 
films.  Titanium  dioxide  is  produced  by  rf  sputtering  a 
titanium  disk  in  a  mixture  of  argon  and  oxygen  at  a 
total  pressure  of  ~0.03  mbar.  The  substrates  are  posi¬ 
tioned  on  a  second  disk,  located  parallel  to  the  first 
disk,  and  ~30  mm  away.  The  substrate  holder  may 
also  be  an  rf  electrode.  Rf  power  applied  to  the  source 
electrode  induces  a  negative  bias  of  ~1500  volts,  caus¬ 
ing  Ar+  ion  bombardment  and  sputtering  of  titanium 
and  oxygen  atoms  from  the  oxidized  titanium  surface 
by  momentum  transfer.  The  atoms  travel  to  the  sub¬ 
strate  with  energies  of  about  10  eV,  and  condense  on 
the  substrate  to  form  a  dense  film.  Activated  oxygen 
produced  in  the  interelectrode  discharge  assures  com¬ 
plete  reaction  and  stoichiometry.  Any  desired  mixed 
phase  composition  (x  anatase,  1  -  x  rutile)  can  be 
made,  from  pure  anatase  to  pure  rutile,  by  selection  of 
the  oxygen  percentage  in  the  sputtering  gas.  Higher 
oxygen  percentages  and  rf  substrate  bias  (increased 
oxygen  activation)  produce  more  rutile  phase.  Grain 
size  is  adjustable  over  the  range  from  amorphous  (<50 
A)  to  600  A  by  selection  of  the  rf  source  power  (sub¬ 
strate  temperature). 

Dip  Coating  (Dip):  Dip  coating  is  entirely  different 
from  the  vacuum  processes  described  previously.  An 
organometallic  compound  (metal  alkoxide)  is  made  in 
solution;  it  contains  titanium  and  an  organic  radical. 
The  sample  is  dipped  into  the  solution  and  the  excess 
liquid  is  drained  off.  The  sample  is  then  exposed  to 
moist  air  to  hydrolyze  the  alkoxide  to  an  oxide.  Next 
the  film  is  heated  to  700-800°  C  to  drive  off  the  remain¬ 
ing  organic  materials  and  to  densify  the  film.  The 
organometallic  solution  can  be  made  to  favor  a  high  or 
low  refractive  index,  dense  or  porous  films.  The  film 
made  for  this  study  was  dip  coated  from  a  solution 
designed  to  produce  low  index  anti-reflection  coatings. 

III.  Measurements  and  Analysis 

A.  Optical  Constants 

1.  Wideband  Spectrophotometric  Method 

The  method  used  by  Group  A  at  Marseille  to  obtain 
n,  k,  and  the  inhomogeneity  A n/n  has  been  described 
previously18"19  and  is  summarized  here.  Reflectance 
and  transmittance  were  measured  at  151  different  wa- 
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velengths  in  the  spectral  range  400-1000  nm  using  a 
double-grating  spectrometer  designed  and  construct¬ 
ed  specifically  for  measurements  of  optical  constants. 
The  average  angle  of  incidence  was  2°,  and  the  optical 
beam  was  close  to  being  parallel  (// 20).  The  illuminat¬ 
ed  area  on  the  sample  was  a  rectangle  4  X  6.8  mm2. 
The  R  and  T  measurements  were  equally  spaced  in 
wavenumbers  (cm-1)  so  the  wavelength  intervals  var¬ 
ied  from  1.6  nm  at  the  short  wavelength  end  to  9.9  nm 
at  the  long  wavelength  end. 

Since  the  substrates  were  not  wedged  as  they  were 
for  the  measurements  in  Ref.  16,  a  correction  had  to  be 
made  for  the  reflection  from  the  back  surface  of  the 
substrate.  An  identical  uncoated  substrate  was  used 
for  this  purpose.  The  measurements  and  analysis  are 
as  follows:  The  transmittance  r0  of  the  bare  nonab¬ 
sorbing  substrate  is 

to-Tyr/d-RA).  (U 

where  T0  and  Rq  are  the  transmittance  and  reflectance, 
respectively,  at  the  front  surface,  and  T,  and  R,  are  the 
transmittance  and  reflectance  at  the  rear  surface  (see 
Fig.  1).  The  transmittance  tc  of  the  coated  substrate  is 

r,  -  TCTJ( l  -  RJ),),  (2) 

where  Tc  and  Rc  are  the  transmittance  and  reflectance, 
respectively,  at  the  surface  coated  with  the  titanium 
dioxide  film. 

The  reflectance  for  the  bare  substrate  is 


\  SUBSTRATE 

VlLM 


Fig.  1.  Diagrams  showing  the  quantities  used  in  the  calculations  of 
the  optical  constants  of  the  titania  films.  The  Ts  and  Rs  are  intensi¬ 
ty  transmittances  and  reflectances,  respectively. 


Fig.  2.  Derived  transmittance  T,  and  reflectance  ft,  for  IBSD  film 
#  120  measured  by  Group  A. 


Pa  -  R0  +  7*iRy(  1  -  RJi,).  (3) 

and  for  the  film-coated  surface  is20 

p,  -R,  +  TlRJ(\  -r^r,),  (4) 

For  a  nonabsorbing  substrate 

R0+r0-i,  (5i 

from  which  it  follows  that 

T5-RS-(T0-Rg)ir0  +  fl0)  -  T„-R0.  (6) 

At  each  wavelength,  the  quantities  measured  are  rc,  r0, 
pc  and  po.  Then  with  the  ratios  Qt  =  rc/r0  [Eqs.  (1)  and 
(2)),  and  Qr  =  pc/po  [Eqs.  (3)  and  (4)],  and  Eq.  (5),  Rc 
and  Tc  can  be  determined  without  any  approxima¬ 
tions.  The  relation  f'  T,  is 

T,  -  3  +  OR,.  (7) 

where  0  =  ToQtHI  —  RoR,)  and  a  = 

The  expression  for  Rc  is 

R,  -(y$*-|JJR,)/(l (81 

where  y  =  [Ro  +  Rr(To  -  /?o)]/(l  —  RoRr)- 

For  fused  silica  substrates,  the  reflectance  is  given 
by 

R0  -  R,  -  (1  -  n.jVCl  +  n,)3,  (9) 

where  the  refractive  index  ns  of  the  fused  silica  is 
assumed  to  be  known  (from  a  handbook  or  glass  cata¬ 
log).  The  value  for  To  is  obtained  from  Eq.  (5)  if  the 
substrate  is  nonabsorbing.  It  is  also  possible  to  take 
into  account  absorption  in  the  substrate,  but  the  calcu¬ 
lations  and  calibration  are  more  tedious. 

Since  Rc  and  Tc  are  the  reflectance  and  transmit¬ 
tance  of  the  Ti02  film  between  two  semiinfinite  me¬ 
dia — fused  silica  and  air — the  values  obtained  from 
the  above  analysis  can  be  used  for  theoretical  calcula¬ 
tions  of  the  optical  constants  of  the  film.  (The  angle  of 
incidence  in  the  spectrometer  was  2°.  In  the  calcula¬ 
tions  we  assume  the  case  of  normal  incidence.) 

The  accuracy  of  the  measured  reflectance  and  trans¬ 
mittance  of  the  Ti02  films  was  lower  than  it  would 
have  been  had  the  films  been  deposited  on  wedged 
substrates.  On  plane  parallel  substrates  the  uncer¬ 
tainties  in  R  and  T  were  about  ±  0.003. 

An  example  of  reflectance  and  transmittance  data 
for  IBSD  Ti02  f.Im  #  120  is  shown  in  Fig.  2;  digitized 
reflectance  and  transmittance  data  were  recorded  to 
four  decimal  places,  and  were  then  used  to  calculate 
the  optical  constants. 

For  the  calculations,  a  linear  variation  of  the  refrac¬ 
tive  index  is  assumed  throughout  the  thickness  of  the 
layer,  and  the  three  following  parameters  are  deter¬ 
mined  simultaneously: 

(1)  The  mean  complex  index  of  refraction  n*  =  n  - 
ik,  which  is  a  function  of  the  wavelength  \. 

(2)  The  variation  An/n  of  the  refractive  index  in  the 
thickness  of  the  layer,  with  the  sign  convention  An  =  ng 
-  n„  where  n0  is  the  refractive  index  at  the  air-film 
interface  and  n,  at  the  film-substrate  interface. 

(3)  The  geometrical  thickness  d. 
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Fig.  3.  Calculated  refractive  index  versus  wavelength  for  IBSD  film 
#  120  measured  by  Group  A.  The  solid  curve  is  the  average  refrac* 
tive  index,  and  upper  and  lower  dotted  curves  are  no  at  the  air- film 
interface  and  n,  at  the  film-substrate  interface,  respectively.  The 
error  bars  were  calculated  assuming  uncertainties  of  ±  0.003  in  both 
Te  and  Rf. 


For  the  initial  calculations,  approximate  values  for  n 
can  be  obtained  from  Rc  or  Tc.  Values  of  k  were 
calculated  using  the  method  described  in  Ref.  13.  In 
the  inhomogeneity  calculation,  the  film  is  split  into  ten 
sublayers  of  equal  (physical)  thickness  with  different 
refractive  index  n.  The  extinction  coefficient  k  is 
assumed  to  remain  constant  for  each  sublayer.  Scat¬ 
tering  losses  are  not  taken  into  account  in  the  calcula¬ 
tions;  the  value  obtained  for  the  extinction  coefficient 
assumes  that  all  significant  losses  are  caused  by  ab¬ 
sorption.  The  final  values  for  n,  k,  and  A n/n  for  each 
wavelength,  and  the  thickness  d  of  the  film,  were  auto¬ 
matically  calculated  using  the  method  described  in 
Refs.  18  and  19. 

An  example  of  calculated  values  of  n  vs  \  for  TiO? 
film  #  120  is  shown  in  Fig.  3.  The  error  bars  are  based 
on  the  assumption  that  the  errors  in  Rc  and  Tc  were 
each  ±  0.003.  The  upper  dotted  curve  is  the  calculat¬ 
ed  refractive  index  no  (at  the  air-film  interface),  while 
the  lower  dotted  curve  is  the  calculated  refractive  in¬ 
dex  n i  (at  the  film-substrate  interface).  Further  de¬ 
tails  of  the  measurements  and  the  fully  automated 
calculations  are  given  in  Refs.  16-19. 

Finally,  it  is  convenient  to  express  n  as  a  function  of 
wavelength  in  the  form  of  a  Cauchy  function  such  as 

n  -  A  +  B/XJ  +  c/x\  (10) 

where  the  coefficients  A,  B,  and  C  are  obtained  during 
the  calculations.  The  refractive  indices  for  most  films 
measured  by  Group  A  are  given  in  this  form  in  Table 
IV  in  Sec.  IV. 

2.  Modified  Valeev  Turning  Point  Method 

The  method  of  obtaining  n  and  k  used  by  Group  B 
(OCLI  and  Martin  Marietta)  has  also  been  described 
previously.16  However,  the  main  points  are  given 
here.  Each  sample  was  scanned  in  transmission  and 


reflection  over  the  spectral  range  350-900  nm  using  a 
Cary  17-DX  spectrophotometer  interfaced  to  a  Hew¬ 
lett  Packard  9825T  desk-top  computer.  The  sample 
was  in  a  moderately  converging  // 5.6  beam,  and  the 
illuminated  area  on  the  sample  was  a  rectangle  ~4-5- 
mm  long  and  1-mm  wide.  The  wavenumber  spacing 
between  data  points  was  81.6  cm'1.  The  transmit¬ 
tance  was  measured  using  a  procedure  similar  to  that 
used  by  Group  A:  The  transmittance  rc  of  the  coated 
part  and  ro  of  an  uncoated  fused  silica  reference  piece 
(sample  #126)  were  measured,  and  the  ratio  Qr  = 
tc/to  was  used  in  the  analysis.  By  using  this  ratio,  the 
effect  of  substrate  absorption  was  minimized.  The 
reflection-losses  for  the  uncoated  substrate  and  for  the 
coated  sample  were  calculated  using  the  known  refrac¬ 
tive  index  vs  wavelength  for  used  silica.  These  were 
taken  into  account  in  the  analysis. 

The  film  reflectance  was  measured  by  using  a  spe¬ 
cially  designed  single  bounce  reflectance  attachment 
on  the  spectrophotometer.  In  order  to  balance  the 
sample  and  the  reference  beams,  a  similar  single 
bounce  reflectance  attachment  with  an  uncoated  sap¬ 
phire  substrate  was  used  in  the  reference  beam.  The 
light  beam  was  incident  on  the  sample  at  10°  although 
the  deviation  from  normal  incidence  was  ignored. 
Only  the  reflectance  from  the  first  surface  of  the  sam¬ 
ple  was  measured.  The  second  surface  reflection  was 
eliminated  by  attaching  a  fused  silica  wedge  onto  the 
back  of  the  sample  with  index  matching  fluid.  To 
eliminate  the  effects  of  variations  in  detector  sensitiv¬ 
ity,  source  intensity,  etc,  with  wavelength,  an  uncoated 
BK-7  wedge  was  scanned  for  the  background.  A  new 
background  scan  was  run  for  every  four  samples. 

The  reflectance  Rc  of  the  titania  film  was  deter¬ 
mined  from  the  relation 

R, »  eSJe*.  (11) 

where  pr  is  the  measured  reflectance  of  the  titania- 
coated  substrate  and  pb  is  the  measured  reflectance  of 
the  BK-7  reference  wedge.  The  primes  indicate  that 
Pc  and  pb  are  measurements  of  the  reflectance  of  a 
single  surface  and  do  not  include  the  second  surface 
reflection.  The  term  Rb  is  the  theoretical  reflectance 
of  a  BK-7  surface  calculated  from  an  expression  simi¬ 
lar  to  Eq.  (9)  using  the  index  of  the  BK-7  glass. 

Values  of  n,  k  and  the  inhomogeneity  A  n/n  were 
determined  using  a  modified  turning  point  method 
described  in  detail  in  Ref.  16,  as  adapted  from  the  work 
of  Valeev.25-22  This  method  uses  the  transmittance 
values  at  the  turning  points  (maxima  and  minima  in 
the  transmittance  vs  wavelength  curve)  and  gives  the 
optical  constants  only  at  the  half  wave  points.  The 
method  was  extended  to  take  into  account  inhomoge¬ 
neity  in  the  films.  Turning  points  were  determined 
after  the  transmittance  ratio  data  Qt  had  been 
smoothed  using  eleven  points  for  the  smoothing  pro¬ 
cess.23  The  steps  in  the  analysis  were  as  follows: 

(1)  Transmittance  ratios  at  the  minima  were  inter¬ 
polated  to  the  wavelengths  of  the  maxima,  and  were 
designated  Q71nl. 

(2)  The  average  value  of  n  and  an  approximate  val¬ 
ue  of  k  were  determined  at  the  half  wave  points  from 
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QTmax>  Qrmti  X,  and  m  using  the  method  of  Valeev.21’22 
Here  m  is  the  order  of  interference  of  the  maximum  at 
which  the  calculations  are  being  made. 

(3)  A  more  accurate  value  of  k  was  determined  from 
the  relation 


h  ”  nA/lxm,  (12) 

where  the  absorptance  A  is  given  by 

U3) 

/?Cmin  is  the  reflectance  of  the  coated  surface  measured 
at  the  minimum  corresponding  to  the  mth  order  maxi¬ 
mum  of  Qt  and  Tcmu  is  the  maximum  transmittance  of 
the  coating  only,  determined  in  a  manner  analogous  to 
Eq.  (7). 

(4)  The  degree  of  inhomogeneity  A n/n  was  deter¬ 
mined  from  the  reflectance  minima  using  the  relation 

Cm  -  nlR'n,  -  RJ/UARJ,  (14) 

where  A  n  =  no  —  n,-  as  before  and  it  is  assumed  that  the 
average  value  of  n  determined  in  step  (2)  is  given  by  n 
=  (no  +  n,)/2.  The  factor  4.4  has  been  arrived  at 
empirically.16 

(5)  The  physical  thickness  d  of  the  film  was  calcu¬ 
lated  at  each  half  point  from  the  relation 

mX/2  -  nd,  (15) 

where  nd  is  the  optical  thickness  of  the  film. 

Graphs  of  n  and  k  vs  X  are  given  in  Figs.  4  and  5  for 
two  films  prepared  by  ion  beam  sputter  deposition. 
Note  that  the  data  points  for  n  and  k  are  obtained  at 
the  wavelengths  corresponding  to  integral  multiples  of 
a  half  wave  optical  thickness  of  the  films  (maxima  in 
the  transmittance  vs  wavelength  curves).  The  inho- 
mogeneity  for  n  and  uncertainty  in  h  are  determined  at 
the  wavelengths  corresponding  to  odd  integral  multi¬ 
ples  of  a  quarter  wave  optical  thickness  of  the  films 
(minima  in  the  transmittance  vs  wavelength  curves). 
In  Fig.  4  the  film  has  a  small  but  nonzero  k  value  which 
decreases  with  wavelength,  while  the  k  value  for  the 
film  in  Fig.  5  is  large  at  the  shortest  wavelength  but 
smaller  elsewhere.  The  bars  associated  on  the  refrac- 
ti  .e  index  curves  extend  from  n,  at  the  lower  end  to  no 
at  the  upper  end,  indicating  the  degree  of  in  homogene¬ 
ity.  The  bars  on  the  k  curves  indicate  the  uncertainty 
in  the  k  values. 

B.  Rim  Absorption 

In  addition  to  obtaining  values  of  the  extinction 
coefficient  k  from  measurements  of  the  film  reflec¬ 
tance  and  transmittance,  film  absorption  was  mea¬ 
sured  using  photothermal  deflection  spectroscopy.24 
A  schematic  diagram  of  the  apparatus  is  shown  in  Fig. 
6.  The  film-coated  substrate  is  placed  in  a  cell  con¬ 
taining  CC1<  which  has  a  large  temperature  coefficient 
of  the  refractive  index.  An  argon  ion  laser  pump  beam 
(X  »  514  nm)  or  a  focused  Xe  light  source  at  normal 
incidence  heats  the  film.  The  He-Ne  laser  probe 
beam  (X  =  633  nm)  is  directed  through  the  CCU  paral¬ 
lel  to  the  sample  and  very  slightly  displaced  from  it. 
As  the  film  absorbs  the  pump  light  it  conducts  heat  to 


Fig.  4.  Calculated  refractive  index  n  (upper  curve)  and  extinction 
coefficient  k  (lower  curve)  vs  wavelength  for  IBSD  film  #  136  mea¬ 
sured  by  Group  B.  The  vertical  bars  on  the  n  curve  indicate  the 
wavelengths  where  the  inhomogeneity  was  measured;  the  length  of 
the  bars  indicates  the  difference  no  -  n,-.  The  length  of  the  vertical 
bars  on  the  k  curve  indicates  the  uncertainty  in  k. 


Fig.  5.  Calculated  refractive  index  n  (upper  curve)  and  extinction 
coefficient  k  (lower  curve)  vs  wavelength  for  IBSD  film  #079  mea¬ 
sured  by  Group  B.  The  vertical  bars  have  the  same  meanings  as  in 
Fig.  4. 


Fig.  6-  Schematic  diagram  of  the  photothermal  deflection  spec¬ 
troscopy  apparatus. 


the  CCI4,  changing  its  refractive  index.  This  causes 
the  He-Ne  probe  beam  to  deviate  slightly.  The 
change  of  position  is  sensed  by  a  position  detector. 

The  relationship  between  the  deflection  angle 
shown  in  Fig.  6,  and  the  absorptance  of  the  sample  is 
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where  n0  is  the  refractive  index  of  CC1<  at  room  tem¬ 
perature,  (dn/dT)  is  the  change  of  /to  with  tempera¬ 
ture,  L  is  the  interaction  length  between  the  probe  and 
pump  beams,  i.e.,  the  pump  beam  diameter  on  the 
sample,  or  *,  is  the  thermal  conductivity  of  the 
deflecting  medium  or  the  substrate  material,  kd  or  k, 
the  thermal  diffusion  coefficient  of  the  deflecting  me¬ 
dium  or  the  substrate  material,  /o  the  intensity  of  the 
incident  pump  beam  in  W/cm2,  zq  the  distance  be¬ 
tween  the  sample  surface  and  probe  beam,  and  A  the 
absorptance  of  the  sample  (A  =  1  —  R  —  T,  neglecting 
scattering).  There  are  only  two  unknown  quantities  z o 
and  A.  Therefore,  a  measurement  is  first  necessary 
with  a  sample  of  known  absorptance  in  order  to  deter¬ 
mine  zq.  Once  z o  is  known,  the  absorptance  of  other 
samples  can  be  determined. 

To  calculate  k  from  the  measured  absorptance  A  of 
the  titania  film,  the  expressions  given  by  Bubenzer  and 
Koidl25  were  used.  However,  the  film  refractive  index 
and  thickness  are  also  needed.  These  values  were 
obtained  from  the  reflectance  and  transmittance 
curves  measured  by  Group  B  for  the  same  films.  Simi¬ 
lar  results  for  the  k,  values  were  obtained  with  a  pro¬ 
gram26  which  calculates  the  k  value  from  the  measured 
absorptance  and  refractive  index.  The  calculated  val¬ 
ues  of  n  and  t  agreed  with  those  of  Group  B  within  the 
limits  of  the  experimental  error  (somewhat  larger  than 
the  errors  of  Group  B).  In  the  conversion  of  A  to  k, 
extreme  accuracy  is  not  needed  in  n  and  t. 

C.  Surface  Roughness 

Surface  roughness  was  measured  at  two  places  on  all 
films  of  Group  B  using  the  Talystep  Surface  Profiling 
Instrument27  at  the  University  of  Alabama  in  Hunts¬ 
ville.  The  instrument  had  a  diamond  stylu-  with  1  -|im 
radius  of  curvature;  a  2-mg  loading  was  used.  Two 
different  profile  lengths  were  measured  at  each  site:  1 
mm  (1000  Mm)  and  100  Mm.  The  longer  profile  includ¬ 
ed  roughness  from  both  the  titanium  dioxide  film  and 
the  commercially  polished  fused  silica  substrate,  while 
the  roughness  on  the  shorter  profile  was  primarily 
caused  by  the  film.  In  both  cases  there  were  2628  data 
points,  with  sampling  distances  of  0.38  Mm  and  0.038 
Mm,  respectively.  The  lateral  resolution  of  the  1-Mm 
radius  stylus  was  ~0.1-0.2  Mm  for  this  type  of  surface. 

Autocovariance  functions28  were  calculated  from 
the  profile  data,  and  the  autocovariance  lengths  deter¬ 
mined.  The  autocovariance  function  is  a  measure  of 
the  correlation  properties  of  the  surface  roughness.  It 
is  the  product  of  two  copies  of  the  same  surface  profile 
as  one  is  shifted  relative  to  the  other.  The  amount  of 
lateral  shift  between  the  two  profiles  is  the  lag  length. 
The  autocovuricnce  length  is  the  value  of  the  lag  length 
at  which  the  autocovariance  function  drops  to  1/e 
(0.368)  of  its  value  at  zero  lag  length.  Very  short 
autocovariance  lengths,  less  than  0.5  Mm,  indicate  that 
the  granularity  of  the  titanium  dioxide  film  is  the 


dominating  feature.  Longer  autocovariance  lengths 
show  that  the  substrate  roughness  is  more  important. 

Power  spectral  density  functions  were  calculated 
from  the  autovariance  functions.28  These  give  the 
frequency  spectrum  of  the  surface  roughness.  Since 
the  roughness  appeared  to  be  uniform  over  the  film 
surface,  the  surface  statistics  were  averages  oi  mea¬ 
surements  made  at  only  two  places  on  each  film. 

D.  Microstructure 

Film  properties  can  be  measured  in  a  variety  of  ways 
to  show  film  microstructure,  degree  of  crystallinity, 
stress,  adhesion  to  substrates,  etc.  In  the  titania  film 
study,  only  a  limited  analysis  of  film  structure  was 
attempted,  partly  to  avoid  destructive  measurement 
techniques,  and  partly  because  of  limited  time  and 
availability  of  people  to  perform  the  specialized  mea¬ 
surements. 

An  attempt  was  made  to  investigate  film  stress  by 
placing  the  coated  substrates  between  crossed  polariz¬ 
ers,  illuminated  by  an  extended  white  light  source,  and 
looking  for  the  characteristic  cross-shaped  pattern. 
Selected  samples  of  Group  B  were  tested  as  well  as 
uncoated  substrates.  No  evidence  of  stress  could  be 
observed  on  either  the  coated  or  uncoated  substrates. 
Flatness  of  selected  coated  and  uncoated  substrates 
was  measured  in  a  Zygo  interferometer  at  the  Universi¬ 
ty  of  Alabama  in  Huntsville  to  see  if  there  was  evidenc.. 
of  substrate  curvature  produced  by  stress  in  the  *ania 
films.  The  optical  figure  of  the  commercial  grade 
uncoated  substrates  was  several  fringes  from  being  fiat 
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A  Measured  j"d  analvred  at  Ecole  Nationals  Sunene  ■  ,»e  Pvvxi<)iie  de  Marseille  Trance 
B  Measured  at  OCLI,  Santa  Rosa.  CA.  artalvred  at  V  *  i  Marietta  Astronaunc »  Group 
Alt  muerque.  NM 

C.  Measured  and  anal y ted  by  A  SaxrratLi  *rv,  of  Alabama  in  Ituniw  :  ALandfAf. \, 
Rome.  Italy 

C*  n  was  calculated  from  R  and  T  measurements  and  interpolated  toMJ  nr?,  k  measured 
si  514  nm  by  phooihemul  deflection  spectroscopy 
*  Rectangular  substrate  did  not  fit  m  sample  holder 
••  Inhomngeneity  presented  satisfactory  n<aturenscnt 
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(X  =  633  nm),  and  varied  considerably  from  one  sub¬ 
strate  to  another.  Thus  the  surface  flatness  was  not 
good  enough  to  detect  small  changes  produced  by 
stress  in  the  films. 

Two  measuring  techniques  were  used  with  some  suc¬ 
cess  to  investigate  film  microstructure.  They  are  de¬ 
scribed  below: 

1.  Contact  Angle  for  Wetting 

One  way  to  determine  the  type  of  structure  in  titania 
films  is  to  measure  the  contact  angle  for  wetting. 
There  are  two  possible  structures  for  crystalline  tita¬ 
nia — anatase  and  rutile.  The  wetting  properties  of 
these  materials  are  different,  and  are  also  different 
from  the  wetting  of  amorphous  fused  silica.29  Wetting 
is  governed  by  an  interaction  force  at  a  solid-liquid 
interface,  which  is  expressed  by  the  Hamaker  con¬ 
stant.29  Anatase  has  a  constant  of  3.5  X  10"20,  while 
rutile  has  a  constant  of  8  X  10-20;  these  differ  by  a 
factor  >2.  Further,  silica  has  an  intermediate  con¬ 
stant,  6  X  lO'20.  In  order  to  see  whether  differences 
could  be  discerned  between  films  prepared  by  differ¬ 
ent  techniques,  water  drops  were  placed  on  the  sur¬ 
faces  of  Group  B  films  and  photographs  were  taken 
aiming  the  camera  parallel  to  the  film  surface  to  show 
the  drops  in  cross  section.  Contact  angles  were  mea¬ 
sured  on  the  photographs.  Differences  were  seen  be¬ 
tween  films,  and  are  reported  in  Sec.  IV. 

2.  Raman  Spectroscopy 

A  Raman  microprobe  was  used  to  investigate  the 
crystalline  structure  in  the  titania  films.  The  two 
forms,  anatase  and  rutile,  have  Raman  scattering 
peaks  at  different  frequencies,  and  thus  are  easily  dis- 
tingishable.  The  primary  anatase  peak  is  at  143  cm-1, 
while  the  rutile  structure  has  peaks  in  the  vicinity  of 
422  cm-1  and  595  cm-1.30  All  the  films  in  Group  B 
were  studied  with  the  Raman  microprobe  at  the  Uni¬ 
versity  of  Alabama  in  Huntsville.  The  instrument  was 


T.bl«  IV.  Cauchy  C  ontt.nl,  (Et).  <10)1  (ortho  Relractlyu  Indie  o)  Film* 
In  (ho  Wavelength  Rang#  400-1000  tun  Meant? od  by  Group  A;  the 
Wavalanyth  It  fn  4m 
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2.0898 

1289 

3-534 

2.1222 

1353 

11970 

3.780 

3492 

IAD 

024 
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1356 

a  Jobin  Yvon  Instruments  S.A.  MOLE  U1000,  with 
excitation  by  a  Spectra-Physics  Series  2020  argon  ion 
laser.  The  focused  spot  on  the  surface  was  ~5  ^m. 
Stokes  transitions  were  measured  on  the  long  wave¬ 
length  side  of  the  exciting  514.5-nm  argon  line.  No 
rutile  was  seen  in  any  of  the  films,  but  anatase  was 
observed,  as  reported  in  Sec.  IV.  Rutile  has  previously 
been  found  in  laser  annealed  titania  films31  and  films 
subjected  to  laser  damage.30 

IV.  Results 

A.  Optical  Constants 

The  values  of  the  optical  constants  measured  by 
Groups  A  and  B  using  the  techniques  described  in  Sec. 
III.A  are  given  in  Tables  III  and  IV.  Further  informa¬ 
tion  is  shown  in  Figs.  7-11.  Several  interesting  obser¬ 
vations  can  be  made  about  the  values  of  the  optical 
constants,  and  are  discussed  in  this  section. 

First,  in  Table  III  are  listed  values  of  n,  k,  and  A n/n 
for  all  films  at  a  wavelength  of  550  nm,  as  measured  by 
Groups  A  and  B,  respectively.  There  is  excellent 


Fig.  7.  Refractive  index  n  vs  extinction  coefficient  k  for  the  films  in 
Group  B. 


2  39  for  all  films. 
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Fig.  9.  Calculated  refractive  index  vs  wavelength  of  E-B  deposited 
films  measured  by  Group  A. 


Fig.  10.  Caicuiated  refractive  index  vs  wavelength  for  IAD  films 
measured  by  Group  A. 


‘  Fig.  11.  Calculated  refractive  index  vs  wavelength  for  IBSD,  IP, 
ARE  and  dip  coated  films  measured  by  Group  A. 


agreement  in  the  n  values  measured  by  each  group  on 
films  prepared  at  the  same  time.  This  shows  consis¬ 
tency  in  the  film  preparation,  and  that  the  two  differ¬ 
ent  measurement  techniques  were  free  of  systematic 
errors.  All  the  n  values  agreed  within  0.03  except  for 
two  films  which  were  observed  to  be  nonuniform.  The 
refractive  indices  of  the  films  spanned  a  range  from 
2.18-2.61,  indicating  a  large  variation  among  the  vari¬ 
ous  films.  Most  of  the  films  had  very  low/e  values  at  X 
=  550  nm;  in  general  the  k  values  measured  by  Groups 


A  and  B  were  in  agreement.  Small  film  inhomogenei¬ 
ties  An/n  were  measured  in  most  cases,  showing  that 
the  preparation  techniques  gave  uniform,  homoge¬ 
neous  films.  Only  two  films  were  too  inhomogeneous 
to  yield  a  value  of  An/n,  while  a  third  had  a  large 
negative  inhomogeneity,  i.e.,  the  film  refractive  index 
was  highest  at  the  film-substrate  interface.  The  film 
prepared  by  dip  coating  was  initially  coated  on  both 
sides  of  a  1  in  (25  mm)  X  6  in  (152  mm)  fused  silica 
substrate.  However,  the  film  on  one  side  was  re¬ 
moved  with  a  saturated  solution  of  ammonium  bifluor¬ 
ide.  The  part  was  then  cut  into  two  1  in  (25  mm)  X  3  in 
(76  mm)  pieces;  one  of  the  pieces  was  sent  to  each 
gToup  for  measurement.  The  measurement  of  k  and 
A n/n  may  have  been  complicated  because  the  coating 
was  composed  of  several  layers,  perhaps  resulting  in  a 
nonuniform  index. 

In  order  to  determine  whether  there  was  any  relation 
between  the  optical  constants  and  the  method  of  film 
preparation,  n  is  plotted  vs  k  in  Fig.  7  for  all  films  in 
Group  B.  The  n  and  k  values  for  most  of  the  E-B  films 
cluster  around  the  low  end,  while  other  deposition 
techniques  yield  films  with  higher  n  values  and  some¬ 
what  higher  k  values.  The  ARE  and  dip-coated  films 
are  exceptions;  these  had  the  lowest  refractive  indices 
of  all.  The  dip  coating  was  made  from  standard  solu¬ 
tions  used  in  making  production  multilayer  antireflec¬ 
tion  coatings,  and  was  intended  to  have  a  low  refractive 
index.  With  other  liquid  precursor  solutions  dense, 
high  refractive  index  films  can  also  be  made  by  dip 
coating. 

Figure  8  shows  a  histogram  of  the  refractive  indices 
for  all  films  measured  by  Group  A.  Again  the  E-B 
deposited  films  have  the  lowest  n  values,  the  IAD  films 
have  higher  n  values,  and  the  IBSD,  IP  and  RFS  films 
have  the  highest  n  values. 

Other  quantities  were  also  compared:  n,  k  and  An 
vs  the  temperature  of  the  substrate  during  deposition, 
and  n  and  k  vs  oxygen  pressure.  No  correlation  was 
seen  between  any  of  these  quantities  and  the  method 
of  film  preparation. 

The  wavelength  dependence  of  the  refractive  index 
measured  by  Group  A  is  given  in  the  form  of  Cauchy 
constants  in  Table  IV,  and  as  curves  of  n  vs  X  for  the 
different  types  of  films  in  Figs.  9-11.  The  variation  of 
n  with  X  was  almost  identical  for  all  E-B  films  with  the 
exception  of  film  #036.  All  IAD  films  had  similar 
wavelength  dependences  except  for  film  #011.  The 
IBSD  and  IP  films  not  only  had  the  highest  refractive 
indices  but  also  the  highest  short  wavelength  absorp¬ 
tion.  The  wavelength  dependence  for  ARE  film  #057 
was  quite  similar  to  that  for  the  E-B  films.  Dip  coated 
film  #213a  had  a  unique  wavelength  dependence. 

The  wavelength  dependence  of  k  was  also  measured 
by  both  groups.  In  general  k  was  largest  at  the  short¬ 
est  wavelengths.  However,  since  the  curves  of  k  vs  X 
were  quite  variable,  only  two  examples  are  shown 
(Figs.  4  and  5). 

Table  V  lists  the  optical  and  physical  film  thick¬ 
nesses  obtained  from  the  analysis  of  the  reflectance 
and  transmittance  measurements.  The  optical  thick- 
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T*W«  V.  Phy*lc*l  (d)  and  Optical  (nd)  Thicknaaaaa  of  THanfa  Film*  at  X 
■  550  nm  Obtalnad  from  tha  Measurements  of  Groups  A  and  B;  ail 
TTilcknesse*  Ara  in  nm 
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A:  Measured  tnd  analyzed  31  Ecolc  National*  Supeneut*  de  Physique  de  Marseille.  France 
B  Measured  ai  OCLI.  Sanu  Rosa.  CA:  analysed  11  Mamn  Manetu  Astronauts*  Group. 
Albuquerque.  NM 

*  Reel  insular  substrate  did  not  fit  in  sample  holder 

♦  Samples  052  and  037  measured  by  Groups  A  tnd  B 


nesses  nd  are  for  X  =  550  nm  rather  than  for  the  five 
quarter  .wave  interference  maximum  wavelength, 
which  varied  around  an  average  value  of  570  nm.  The 
excellent  agreement  between  physical  and  optical 
thicknesses  obtained  on  pairs  of  films  prepared  at  the 
same  time  again  validates  the  two  analysis  techniques, 
and  shows  that  the  pairs  of  films  had  nearly  identical 
properties.  In  the  only  case  where  there  is  a  large 
discrepancy  between  the  two  samples  in  a  pair  (IBSD 
films  #078  and  #079),  colors  were  clearly  evident  on 
sample  #078  indicating  a  nonuniform  film  thickness. 

Figure  12  shows  a  histogram  of  the  optical  thick¬ 
nesses  for  the  films  in  Group  A.  Twelve  films  had 
optical  thicknesses  lying  between  675  and  775  nm,  4 
films  had  lower  values  and  3  had  higher  values.  The 
average  optical  thickness  for  all  films  was  711  nm;  the 
wavelengths  for  which  the  films  had  a  5X/4  optical 
thickness  were  not  tabulated. 

Physical  film  thicknesses  could  not  be  measured 
with  a  profiling  instrument  because  the  film  edges 
were  not  sharp.  In  most  cases  colored  rings  showed 
that  the  film  thickness  was  decreasing  over  a  distance 
of  several  millimeters.  It  was  also  not  possible  to 
remove  the  film  nondestructively  from  part  of  the 
substrate. 

B.  Absorption 

The  results  of  the  direct  absorption  measurements 
using  photothermal  deflection  spectroscopy  are  pre¬ 
sented  in  Tables  III  (C*)  and  VI.  Since  an  argon  laser 
X  =*  514  nm  was  used  to  heat  the  film,  the  k  values  in 
Table  III  and  the  absorptance  values  in  Table  VI  are 


Fig.  12.  Histogram  of  the  optical  thicknesses  of  the  films  in  Group 
A,  calculated  at  X  *  550  nm. 


Fig.  13.  Absorptance  as  measured  by  photothermal  deflection 
spectroscopy  vs  rms  roughness  for  selected  films  in  Group  B. 


for  that  wavelength.  However,  the  k  values  in  Table 
III,  column  C*  are  in  reasonable  agreement  with  those 
calculated  by  Groups  A  and  B  from  spectrophotomet- 
ric  measurements.  The  two  films  having  the  highest  k 
values  from  optical  measurements  (IBSD  film  #122 
and  RFS  film  #  145)  also  have  high  k  values  calculated 
from  photohermal  deflection  spectroscopy  measure¬ 
ments.  The  agreement  for  E-B  film  #  023,  which  also 
had  a  high  measured  absorption,  is  not  so  good.  Later 
measurements  on  other  titania  films  made  at  several 
wavelengths  indicate  that  there  is  a  20-50%  decrease  in 
k  between  514-550  nm. 

The  relation  between  the  surface  roughness  of  some 
films  and  their  absorption  at  X  =  514  nm  is  shown  in 
Fig.  13.  As  discussed  in  Sec.  IV. C,  the  short-profile- 
length  roughness  is  characteristic  of  the  film  surface 
only.  The  plot  of  the  film  absorptance  (as  determined 
by  photothermal  deflection  spectroscopy  at  514  nm)  vs 
the  surface  roughness  of  the  100-utn  profiles  (from 
Table  VI)  shows  clearly  that  the  rougher  films  also 
have  more  absorption.  Since  only  two  films  prepared 
by  each  technique  were  measured  by  photothermal 
deflection  spectroscopy,  it  is  not  known  whether  there 
is  a  linear  relation  between  surface  roughness  and  ab- 
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sorption.  However,  it  is  clear  that  the  dependence  of 
absorption  on  surface  roughness  is  different  for  the 
different  types  of  films.  The  ion  plated  films  show  the 
least  increase  in  absorption  with  increasing  roughness. 
A  possible  explanation  for  these  results  is  that  there  is 
an  enhanced  coupling  of  the  pump  light  energy  (at  514 
nm)  for  the  rougher  surfaces  because  of  enhanced  scat¬ 
tering  into  the  Ti02  film.  The  different  slopes  for  the 
different  deposition  techniques  seem  to  be  related  to 
the  microstructure  of  the  films.  The  fairly  porous 
columnar  structure  known  for  E-B  films  provides  for 
an  abundance  of  foreign  adsorbates  in  the  film,  where¬ 
as  the  almost  perfectly  dense  EP  films  prevent  such 
intrinsic  absorption,  and  these  films  stay  clean. 


Ttbla  VI.  Measured  Roughness,  Correlation  Length,  and  Absorptance  of 
Tlfanla  Film*  In  Group  B 
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Fig.  15.  Surface  profile  (1000-moi  profile  length)  for  a  typical 
uncoated  fused  silica  substrate  in  Group  B. 
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Fig.  16.  Surface  profile  (100-jim  profile  length)  for  IBSD  film 
#i36. 


Fig.  14.  Measured  rms  roughnesses  of  films  in  Group  B. 


Fig.  17.  Surface  profile  (100-Mm  profile  length)  for  IBSD  film 
#122. 
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C.  Surface  Roughness 


AUTOCOVARIANCE  FUNCTION 


The  measured  roughnesses  and  autocovariance 
lengths  of  the  films  in  Group  B  are  listed  in  Table  VI. 
The  films  are  in  the  order  of  increasing  roughness. 
Roughness  values  for  the  1000  ^m  profile  length  in¬ 
clude  both  the  substrate  roughness  and  the  added 
short  spatial  wavelength  roughness  of  the  titania  films. 
Film  roughness,  when  present,  generally  dominates  in 
the  lOO-ftm  profiles.  Film  roughnesses  (100-Mm  pro¬ 
file  length)  are  shown  for  the  different  film  types  in 
Fig.  14.  It  is  clear  that  no  one  type  of  film  preparation 
method  gives  the  smoothest  films;  smooth  films  (and 
also  rougher  ones)  were  produced  by  all  methods. 

Figures  15-21  show  selected  surface  profiles,  as  well 
as  autocovariance  functions  and  power  spectral  densi¬ 
ty  functions  obtained  from  the  profiles.  A  1000-^m- 
long  profile  of  one  of  the  uncoated  commercially  pol¬ 
ished  fused  quartz  substrates  is  shown  in  Fig.  15.  The 
roughness,  0.58  nm  rms,  is  typical  is  the  roughness 
values  of  the  bare  substrates  which  ranged  from  0.4- 
0.8  nm  rms;  most  were  about  0.6-0.7  nm  rms.  Rough¬ 
nesses  for  the  100-Mm  profile  lengths  were  in  the  0.2- 
0.4  nm  rms  range,  with  0.3  nm  rms  being  a  typical 
value.  These  values  are  inferred  from  a  careful  exami¬ 
nation  of  all  the  profiles. 

In  Figs.  16  and  17  profiles  are  shown  for  one  of  the 
smoother  and  one  of  the  rougher  films,  both  of  which 
were  made  by  ion  beam  sputter  deposition.  The  pro¬ 
file  of  the  smoother  film  is  indistinguishable  from  that 
of  a  bare  substrate,  while  the  profile  of  the  rougher  film 
clearly  shows  the  grain  structure  of  the  film. 

The  autocovariance  functions  for  the  same  smoother 
and  rougher  films.  Figs.  18  and  19,  have  very  different 
characters,  especially  with  regard  to  the  initial  shapes 
of  the  curves.  Note  that  the  autocovariance  lengths 
differ  by  nearly  10  times.  The  short  autocovariance 
length  for  the  rougher  film  is  typical  of  film  granulari¬ 
ty,  while  the  longer  autocovariance  length  for  the 
smoother  film  is  caused  by  the  substrate  roughness, 
with  the  film  adding  no  additional  structure. 

Power  spectral  density  functions  for  the  same  two 
films  are  shown  in  Figs.  20  and  21.  A  vertical  arrow 
marks  the  spatial  frequency  corresponding  to  the  auto¬ 
covariance  length  of  each  film.  The  power  spectral 
density  function  for  the  smoother  film  is  indistinguish¬ 
able  from  that  of  an  uncoated  substrate.  It  has  a 
maximum  contribution  from  the  lowest  spatial  fre¬ 
quencies  (long  spatial  wavelengths).  The  highest  spa¬ 
tial  frequencies  are  less  than  10-4  of  the  maximum 
lowest  spatial  frequency.  By  contrast,  the  power  spec¬ 
tral  density  function  for  the  rougher  film  has  contribu¬ 
tions  from  a  larger  range  of  spatial  frequencies  with 
much  less  decrease  for  high  spatial  frequencies  (short 
spatial  wavelengths). 

The  range  of  spatial  frequencies  that  would  produce 
scattering  in  the  visible  spectral  region  for  a  He-Ne 
laser  beam,  X  =•  0.633  nm,  at  normal  incidence  and 
scattering  angles  of  0.4-90°  is  shown  by  the  horizontal 
double  ended  arrows  in  Figs.  20  and  21.  It  is  clear  that 
the  rougher  film  should  scatter  much  more  than  the 


Fig.  18.  Autocovariance  function  for  IBSD  film  #  136. 


AUTOCOVARIANCE  FUNCTION 


Fig.  19  Autocovariance  function  for  IBSD  film  #  122. 


SPATIAL  FAf OUlMCY  (micron*  ') 


Fig.  20.  Power  spectral  density  function  for  IBSD  film  #  1 36.  The 
vertical  arrow  marks  the  spatial  frequency  corresponding  to  the 
autocovariance  length.  The  horizontal  double  ended  arrow  marks 
the  range  of  spatial  frequencies  that  produce  scattering  into  angles 
between  0.4*  and  90®  for  a  waveh  igth  of  0.633  nm  at  normal  inci¬ 
dence. 
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smoother  one.  Studies  of  angle  resolved  scattering  of 
the  Group  A  samples  have  been  made  by  Group  A. 
However,  since  their  He-Ne  laser  spot  size  was  ~4  mm 
in  diameter,  scattering  from  substrate  defects  domi¬ 
nated.  The  results  of  these  measurements  will  be 
reported  in  a  forthcoming  paper.32 

0.  Microstructure 

The  measured  contact  angles  for  wetting  and  film 
crystallinity  as  determined  from  Raman  microprobe 
measurements  are  given  in  Table  VII.  The  cosine  of 
the  wetting  angle  is  plotted  vs  rms  roughness  in  Fig.  22. 
(Films  that  contained  anatase  are  indicated  by  the 
solid  circles  in  the  figure.)  The  cosines  of  the  wetting 
angles  divide  approximately  into  two  groups,  those 
which  are  larger  than  that  of  an  uncoated  fused  silica 
substrate  and  those  which  are  smaller.  (Large  cosines 
correspond  to  small  contact  angles.)  It  is  immediately 
obvious  that  there  is  no  consistent  relation  between 
the  method  of  film  preparation  and  the  contact  wet¬ 
ting  angle.  Also,  the  film  roughness  is  not  correlated 
with  the  contact  wetting  angle. 

Film  crystallinity  as  measured  by  the  Raman  micro¬ 
probe  showed  that  some  of  the  films  contained  ana¬ 
tase,  but  none  contained  rutile;  most  were  amorphous. 
In  Fig.  23  are  shown  Raman  microprobe  spectra  for  two 
IAD  titania  films;  one  had  a  large  anatase  peak  and  the 
other  had  none  so  is  presumed  to  be  amorphous.  Only 
eight  of  the  films  in  Group  B  showed  any  evidence  of 
anatase,  ranging  from  large  peaks  for  E-B  Him  #023 
and  IAD  film  #012,  to  medium  peaks  for  E-B  films 
#038,  #202,  and  #005,  and  RFS  film  #145,  and 
small  peaks  for  IAD  film  #014  and  dip-coated  film 
#213b.  As  expected,  the  uncoated  fused  silica  sub¬ 
strate  showed  no  evidence  of  crystallinity. 

There  may,  however,  be  a  relationship  between  the 
contact  wetting  angle  and  the  film  crystallinity,  as 
predicted  by  the  Hamaker  constant  mentioned  earlier. 
Five  of  the  eight  anatase-containing  films  had  large 
contact  wetting  angles  (small  cosines);  these  included 
all  the  films  whose  roughnesses  were  0.5  nm  rms  or 
greater  except  for  IBSD  film  #  122,  which  was  rough 
but  contained  no  trace  of  anatase.  On  the  other  hand, 
four  of  the  seven  films  that  had  small  contact  wetting 
angles  (large  cosines)  contained  no  anatase.  Excep¬ 
tions  were  E-B  films  #023  and  #005  and  dip-coated 
film  #213b.  Interestingly,  half  of  the  anatase  con¬ 
taining  films  were  prepared  by  E-B  evaporation.  The 
relation  between  contact  wetting  angle  and  film  crys¬ 
tallinity  is  by  no  means  conclusive,  but  it  appears  that 
anatase-containing  films  are  likely  to  have  larger  con¬ 
tact  wetting  angles  (smaller  cosines)  than  uncoated 
fused  silica. 

V.  Conclusions 

There  were  large  variations  in  the  properties  of  films 
produced  by  each  technique.  Film  refractive  index, 
absorption,  roughness,  contact  angle  for  wetting,  and 
crystallinity  seemed  to  depend  not  only  on  the  tech¬ 
nique  but  also  on  the  particular  coating  chamber  and 
parameters  used.  However,  some  general  conclusions 


ROWER  SPECTRAL  OENSITT  FUNCTION 


SPATIAL  FREQUENCY  (morons  ■’) 

Fig.  21.  Power  spectral  density  function  for  IBSD  film  #  122.  The 
arrows  have  the  same  meanings  as  in  Fig.  20. 


Fig-  22.  Cosine  of  the  contact  angle  for  wetting  vs  rms  surface 
roughness  for  the  films  in  Group  B  and  a  bare  fused  silica  substrate. 
The  solid  data  points  are  for  films  containing  anatase  (see  Table 
VII). 


Fig.  23.  Raman  microprobe  spectra  for  IAD  films  #012  and  #022 
in  Group  B.  Note  the  anatase  peak  in  the  upper  curve. 
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TaWe  VIL  Roughness,  Contact  Angle  for  Wotting,  and  Crystallinity  of 
TllanJa  Film*  in  Group  B 
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Alabama  m  Huntsville.  AL 

Film  crystallinity  measured  by  Raman  mcroprobe  by  Boon  Loo  and  Oavid  Bums  at  University 
of  Alabama  m  Huntsville,  AL 


can  be  drawn:  The  refractive  indices  of  the  E-B  films 
as  a  group  were  lowest,  IAD  films  had  intermediate 
refractive  indices,  and  the  IBSD  and  IP  films  had  the 
highest  indices  (Figs.  8-10).  The  short  wavelength 
refractive  index  was  only  slightly  higher  than  at  550 
nm  for  the  E-B  films;  there  were  larger  differences  for 
the  IAD,  IBSD,  and  IP  films.  There  was  a  partial 
correlation  between  film  crystallinity  and  contact  wet¬ 
ting  angle  in  that  nearly  all  of  the  films  with  large 
contact  wetting  angles  also  contained  crystalline  ana- 
tase  (Fig.  22).  Also,  all  but  two  of  the  films  rougher 
than  0.6  nm  rms  contained  anatase  (Fig.  22). 

If  one  defines  the  best  films  by  some  criterion,  for 
example,  those  with  the  highest  refractive  index,  low¬ 
est  absorption,  smallest  inhomogeneity,  or  smoothest 
surface,  there  is  no  one  technique  that  will  consistently 
give  this  type  of  film.  All  techniques  produced  films 
with  some  excellent  qualities,  but  best  and  worst  films 
were  produced  by  each  technique. 

Excellent  agreement  was  obtained  for  values  of  opti¬ 
cal  constants  and  film  thickness  for  the  corresponding 
films  measured  by  Groups  A  and  B,  indicating  that 
there  were  minimum  errors  in  the  measurement  tech¬ 
niques  and  a  high  degree  of  consistency  between  pairs 
of  films  prepared  at  the  same  time.  In  the  few  cases 
where  agreement  was  not  obtained  between  two  films 
in  a  pair,  the  films  had  noticeable  color  variations 
across  their  surfaces,  indicating  nonuniform  optical 
properties. 

Finally,  it  should  be  emphasized  that  all  films  were 
provided  at  the  request  of  the  organizers  of  this  study, 


and  were  not  company  products  or  in  any  way  repre¬ 
sentative  of  the  best  that  a  group  could  do.  In  some 
cases  the  films  were  the  first  of  their  kind  produced  by 
a  group,  and  better  ones  have  been  made  since.  Many 
of  the  techniques  used  to  produce  the  titania  films  are 
new,  and  much  better  films  are  expected  to  be  pro¬ 
duced  in  the  future. 
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